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Chapter I 

GENERAL CONSIDERATIOre 


1.1 Introduction 

In Loth external and Internal aerodynamics, most of the problems 
involve boundary layer flows over curved surfaces. Some of these flows 
are fully developed, while others are not. In the past, investigators 
have been content to treat these boundary layer flows by assuming 
that the curvature of the mean flow streamlines in planes normal to 
the surface, and the associated static pressure variation across the 
boundary layers do not influence the flow significantly, if the radius 
of curvature of the surface is much larger than the boundary layer 
thickness. With such simplification, the many problems in viscous 
aerodynamics can be solved by making use of the vast amount of data 
collected in the study of fully developed pipe /channel flows and two- 
dimensional boundary layers. This assumption is known to be correct 
approximately in laminar flow, since the effects of the additional 
curvature of the mean flow streamlines are known to be of second order 
smallness (see van Dyke 1962). These small effects arise from the 
kinematics of curved flow, and the curvature of mean flow streamlines 
has little influence on the magnitude and distribution of the viscous 
stresses . 

On the other hand, turbulent flows are very sensitive to the 
curvature of the mean flow streamlines. The additional curvature 
influences the flow in a manner which is not known at present, but 
is believed to be significant enough to change the detail flow mechanism 
completely. This is evident from an examination in some detail of the 
observed phenomena in a turbulent boundary layer on a flat surface . 

Very close to the wall, there appears to be a region where the 
fluid motion is still predominantly viscous, and the velocity rises 
steeply. Further away, the flow becomes unstable, and finally, a 
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region is reached where the entire flow is involved in turbulent motion. 
Recent detailed experimental studies by Schraub and Kline (1965) 
have revealed that the viscous region is not truly undisturbed. 

Rather, relatively large elements of low velocity fluids adjacent to 
the surface of the viscous sublayer periodically lift off the surface 
and move towards the fully turbulent region. There they join a pattern 
of decaying turbulence. The mechanism responsible for this phenonfenon 
is not fully understood. However, it can be attributed to tte insta- 
bility of the flow near the outer edge of the viscous sublayer, and also 
due to the action of the outer flow which operates mainly by pressure 
fluctuations transmitted to the sublayer. Since continuity dictates 
that fluid must replace the elements which lift off the surface, elements 
of highly energetic fluid will have to move in from further out. In 
turn, these highly energetic fluid may supply the energy to eject the 
elements from the surface. At any rate, turbulence in the f ull y 
turbulent region is generated and maintained by the elements originating 
at the surface . 

With this picture in mind, one can see why turbulent flows are 
very sensitive to curvature of the mean flow streamlines. For flows 
over convex surfaces, the centrifugal force cn a fluid particle must 
be balanced by an inward pressure gradient. If a particle is moving too 
slowly, its centrifugal force is too small, and it moves inwards. 

As a result, the fluid elements, which lift off from the surface carrying 
with them the velocities at the points where they come from, will have 
a smaller centrifugal force at their new locations. Therefore, they 
will be pushed inwards by the pressure gradient and the interchange of 
momentum and energy between the faster and the more slowly moving 
fluids are hindered. As a result, the boundary layer is thinner at 
the wall. Observations of flows between rotating cylinders by Taylor 
( 1936 ) confirmed this and showed a very large reduction in turbulent 
shear stress. In contrast with this, the destablizing effect of 
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centrifugal forces on concave walls leads to the formation of Taylor - 
Gortler type vortices which promote transition to turbulence. Therefore , 
it is quite possible that turbulent boundary layers on curved surfaces 
may differ appreciably from that on a flat plate with the same distri- 
bution of external pressure. 

Before proceeding any further , it is felt that a discussion of 
previous work is most essential in that it not only points the way 
for further work in curved turbulent flows, but that it also indicates 
the logical line of attack for the present investigation. For purpose 
of clarity, the discussion is divided into two sections. The first 
section deals mainly with fully developed flows, while the second touches 
on turbulent boundary layers along curved surfaces. In anticipation of 
the fact that the flow phenomenon in turbulent boundary layers along 
concave surfaces are different from that along convex surfaces, distinc- 
tion will be made between these two types of flow in the discussion 
below. 

1.2 A Brief Discussion of Previous Work 

1. 2.1 Fully Developed Curved Turbulent Flows 

Since the flow near a surface is determined to a large extent by 
local conditions and to a lesser extent by the action of the outer flow 
which operates mainly through pressure fluctuation transmitted to the 
region near the surface, one can expect the flow in the vicinity of the 
wall to be the same, be it a boundary layer flow or a fully developed 
flow. This is true for straight flows and data obtained from pipe 
flow measurements lead to the establishment of the Law of the Wall. 

With only the additional effects of the curvature of the mean flow 
streamlines which are known to be of second order smallness only if the 
flow is viscosity dominant, one would expect the curvature to have very 
little effect on the flow close to the wall*. As a result, one can expect 

*This fact is at least partially supported by the vast amount of flat plate 
pressure gradient data in which the mean flow streamlines are also curved 
due to rapid boundary layer growth. 



the Lav of the Wall to hold also in a curved flow. However , the 
measurements of Watte ndorf (1935) in two curved channels of constant 
curvature and a ratio a of half channel width to mean radius of 
curvature of l/l9 and 1/9 do not lend evidence to such an argument. 

Rather, Wattendorf found that the deviation from the Law of the Wall was 
in opposite directions for opposite signs of curvature, and that the 
deviation increased with curvature. 

Eskinazi and Yeh ( 1956 ) did similar measurements in a curved 
channel of constant curvature and with a = l/l9 . They found that the 
deviation from the Law of the Wall does not start until about ^t/v - 200. 
Indeed, the deviation was in opposite directions for opposite signs of 
curvature. This did seem to indicate the existence of the Law of the Wall 
region. From their data, Eskinazi and Yeh concluded that u/u T is not 
a function of ^*t/v , but also of y/r or some combination of the 

two. 

For fully developed turbulent flow in a curved channel of constant 
curvature, the turbulent shearing stress would, in the momentum transfer 
theory, be given by 


- uv = v {— — (ur)) 
e r &r 


where v is the eddy viscosity. This requires the shearing stress to 
6 c* 

vanish at Hur) = 0 . However, both Wattendorf and Eskinazi and Yeh 
or 

found that the point of zero uv did not coincide with the point of 
d , V 

zero c-KurJ . Hence, Eskinazi and Yeh concluded that there was a 
dr 

region between these two points in which the product uv r—(ur) was 

dr' — 

positive. This meant that in this region, the turbulent energy in u 2 
was being suppressed by Reynolds stress working on the mean momentum 
gradient . 

Besides mean flow measurements Eskinazi and Yeh also made detailed 
turbulence measurements and found that the turbulence intensities were 
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greater on the outer (concave) vail and smaller on the inner (convex) 
wall compared to corresponding points in a straight channel. In an 
in viscid, irrotational curved flow , the radial equilibrium of a fluid 
element is stable if the radial gradient of angular momentum is negative. 
Eskinazi and Yeh were the first to show that these considerations also 
apply to curved turbulent flows . 

1 .2.2 Turbulent Boundary Layers on Curved Surfaces 

The effect of curvature on turbulent boundary layer was first 
investigated by Wilcken (1930) on both the convex and concave surfaces. 

It was then followed by Schmidbauer (1936) "who studied the flow over 
convex surfaces only. In both these experiments, s/r was at most .025 . 
From his measurements, Wilcken concluded that the eddy viscosity was 
much greater for the flow near the outer (concave) wall than that near 
the inner (convex) wall, and that the mixing length at the outer wall 
was found to be considerably larger than for a flat surface. Although 
no turbulence measurements were made, these results seem to lend evidence 
to the fact that with concave surface, the turbulence is strongly increased 
while with convex walls, it is reduced. Wilcken explained this by arguing 
that the centrifugal forces at the outer wall promote, while it diminishes 
near the inner wall, the turbulent exchange between adjacent fluid layers. 
Hence, the rapid thickening of the boundary layer on the outer wall is 
due also to the increased turbulence activity. The opposite is true on 
the convex wall. In connection with this, Schmidbauer 's data also 
showed a decrease in the boundary layer thickness. However, Wilcken 's 
results were somewhat obscured by disturbing side influences, which 
were mainly caused by pressure conditions at the beginning and end of 
the channel. In addition, the disturbing effect of the secondary flow 
may also have influenced the boundary layer developnent on both the inner 
and outer walls of the channel. 

The work of Schneider and Wade (1967) did nothing to eliminate the 
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secondary flow influences. In fact, it is quite a bit more pronounced 
because the aspect ratio of the test section varies from 1 at the 
entrance to 2 at the exit. Also, the flow in the test section was one 
of constant acceleration. As a result, it was little wcnder why their 
data do not correlate well with the Law of the Wall. 

Through detailed measurements in a l 8 (f bend, Patel ( 1968 b) 
demonstrated that the boundary layer on the convex wall was relatively 
two-dimensional near the center line for the first 80 ° of the bend. 

Thus, he managed to establish that the boundary layer development on 
the convex wall near the center line of an 80 ° curved duct was relatively 
from the influence of secondary flow. Xn contrast to previous work, 
the boundary layer investigated by Patel ( 1968 a) has a max imum value of 
6 /r ~ 0.1 . However, Patel did nothing to change the pressure distri- 
bution on the convex wall of the curved duct, and in the opinion of the 
present author, his data is subjected to the same shortcomings as that 
of Wilcken in that it is hard to separate the curvature effects from 
the pressure gradient effects. 

The flow between two concentric cylinders with the inner one in 
motion while the outer one remains at rest affords a good example of 
the unstable stratification causes by the additional centrifugal forces 
acting on the fluid particles. Since the velocity is higher near the 
inner cylinder, the fluid particles experience a higher centrifugal 
force, and as a result, the fluid particles will have a tendency to 
move outward. Taylor (1938) was the first to observe that when a certain 
Reynolds number is exceeded, longitudinal vortices with axes located 
along the circumference begins to appear. A similar vortex system was 
being observed by Gortler (1940) for flows along concave walls. The 
analogy can also be drawn between the flow along a concavely curved 
surface and that along a heated horizontal plate. Gortler (1959) 
showed that the bouyant force in the thermally stratified layer also 
gives rise to a system of longitudinal vortices. Noting that the 



mechanism responsible for the existence of the longitudinal vortex system 
is the same for these three different flows, one will naturally ask 
whether such a system would also exist in a fully turbulent flow along 
a concave surface , and if so, would the system of vortex be stationary. 

In addition, one would like to find out the effect of such a vortex 
system on the development of the boundary layer along the concavely 
curved surface . 

Tani (1962) was the first to demonstrate the existence of the 
longitudinal vortices in the turbulent boundary layer along a concave 
wall. According to Tani, the resultant wave amplitude of the vortices 
decreased as the radius of curvature of the surface was increased. More 
recently, Patel (1968b) found that the longitudinal vortices also 
appeared on the concave wall of his l8Cf channel, and suggested that 
before studying the influence of concave curvature on a two-dimensional 
turbulent boundary layer , one should examine these vortices in detail 
so as to determine the curvature parameter that would govern the appearance 
and strength of these vortices. 

1.3 Present Objectives 

A great many mathematical models of physical phenomena are very 
complicated and are highly nonlinear in nature. With the present 
knowledge of nonlinear equations, most of these models cannot be 
solved in full. Usually, simplifying assumptions are necessary to 
reduce these equations to a more manageable form. In so doing, the 
researcher is faced with the difficulty of deciding which effects or 
variables are more important. Experimental studies are designed to 
provide such information. Not only will the experimental results 
reveal the most significant features of a physical phenomenon (provided 
the experimental work is designed and carried out properly), they will 
also provide empirical correlations for the theoretical researcher, 
thus, enabling the researcher to make further simplifications of the 


7 



mathematical model and eventually obtain a solution to the particular 
problem . 

In connection with this, Coles (1962) has pointed out that the 
most useful data are those in which one physical effect or parameter is 
varied at a time. Well designed experiments should be able to separate 
the various physical effects and isolate them for Investigation, thus 
allowing the researcher to comprehend the significance of these various 
effects easily. Fully developed turbulent flows in pipes and channels 
and flat plate turbulent boundary layer with the pressure gradient 
carefully controlled, have led Clauser (1954, 1956) to the formulation 
of the concept of the equilibrium boundary layer. 

In keeping with the idea that one variable should be varied at 
a time, the present experiments are so designed that curvature alone is 
isolated for investigation. Due to the presence of the pressure gradient 
across the boundary layer in flows along curved walls, it is hard to 
avoid the streamwise pressure gradient in the flow along a flat surface 
joining smoothly to a curved surface. However, the pressure change in 
the streamwise direction should be minimized, so that the flow can 
recover from the pressure effect in a relatively short distance. Means 
of doing this are incorporated in the boundary layer tunnel that is 
specially designed for the present investigation. (This will be discussed 
further in Chapter II.) Provisions are also made for the variation of 
streamwise pressure distribution along the curved wall, so that both 
constant pressure and pressure gradient data can be obtained. Because 
of the need to be able to establish constant pressure on both convex 
and concave surfaces, it is apparent that a single test section cannot 
provide such versatility. Therefore, separate test sections, one with 
convex curvature and one with concave curvature, are required. Again, 
previsions are made in the tunnel so that the test section can be 
changed without too much trouble. No attempt has been made to obtain 
data for several values of the curvature parameter, instead it has 
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been decided to study the flow phenomenon in more detail in addition to 
the usual velocity profile measurements. 

Thus, the objective of the present investigation is to provide 
both empirical knowledge and physical understanding in the following 
areas. These areas are: (a) constant pressure flow along convex surface, 

(b) separating flow along convex surface, and (c) constant pressure flow 
along concave surface. In case (a), the curvature parameter is kept 
constant, but is varied along the flow in case (b). For both cases (a) 
and (b), velocity profiles were obtained from pitot measurements across 
the boundary layers while a rotating wire technique proposed by Fujita 
and Kovasznay (1968) was used to measure the various components of the 
Reynolds stress tensor. 

In view of the fact that the flow along a concave surface is 
different from that along a convex surface, the system of longitudinal 
vortices that exists was mapped out first using fixed hot-wires. Then, 
the rotating -wire technique was used to measure the various components 
of the Reynolds stress tensor at both the positions of the crest and 
trough of the wave system. 

Following a discussion of the various test equipment and their 
qualification in Chapter II, the results are examined and critically 
analyzed in Chapter III, and the Law of the Wall is established for 
flows along convex surfaces. In Chapter IV, a self-consistent set of 
equations for curved turbulent boundary layer flows is derived, together 
with an eddy viscosity which is modified to include curvature effects. 

The predictions by this model are then compared with the experimental 
results. Chapter V summarizes all the results and recommendations for 
future work are also given . 
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Chapter II 

EXPERIMElfEAL PROGRAM 


2.1 Curved Wall Tunnel 

2.1.1 Description of Tunnel 

The construction of the curved vail tunnel used in the present 

investigation is described in detail in Appendix A. The diagrammatic 

layout of the tunnel is shovn in Figure 1 and a picture of the curved 

wall tunnel is given in Figure 2. Briefly, the wind tunnel is of the 

open-return suction type and is powered by a two speed, 10 H.P. fan. 

The Reynolds numbers per unit length corresponding to these two speeds 

5 5 

are 4.37 x 10 and 7*56 x 10 respectively. The wind tunnel is less 
quiet and the vibration due to the fan is more severe at the higher 
speed. As a result, all runs on both the convex and concave test vails 
were conducted at the lever Reynolds numbers. Originally the entry to 
the tunnel was made up of three layers of honey-comb flow straightener 
separated by screens. However, it was found that the unit cannot provide 
sufficient straightening effect to prevent unsteadiness in the flow. 
Later, through observation of tufts, the unsteadiness was shown to be 
due to the unsteady ingestion of vortices which foimed on the floor of 
the room. To remedy the situation, two more layers of honey-comb 
straightener separated by screen together with a bell -mouth were added 
to the entrance of the tunnel. The traversing of tufts behind the 
last screen then showed that the entire unit has sufficient straightening 
effect to prevent unsteadiness in the flow, and that the flow is quite 
parallel . 

A contraction section that has a contraction ratio of 6 to 1 is 
installed behind the flow straightener. The straight section that 
follows is 4 feet long and has a cross-section of 6 inches by 48 inches. 
The remainder of the tunnel consists of the curved test section which 
has a fixed test wall of varying curvature and an adjustable opposite 
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vail, the exit diffuser and the fan and motor housing. The adjustable 
wall allows the pressure distribution on the curved test wall to be 
adjusted to give any arbitrary pressure distribution. 

Two curved test sections were made; one has a convex test wall, 
while the other has a concave test wall. The whole tunnel was so constructed 
that the test section could be connected to and disconnected from the straight 
section and the exit diffuser with relative ease. In order to accommodate 
the adjustable wall in the curved test section, the exit diffuser was 
constructed in such a manner that one movable wall was all that was 
required to fit the test section with either the convex or concave test 
walls. End wall jets were installed at the entrance to the curved 
test section to limit the secondary flow in the test section to small 
regions near the end walls. In addition, it was found necessary to have 
a side wall jet installed on the wall opposite the convex test wall. 

However, with a concave test wall, the side wall jet was found unnecessary. 
Further discussions of the significance of the side wall jet and the end 
wall jets will be given in Sections 2 . 3-2 and 2.3.3 respectively. 

Wall static pressure taps made of .032 inches O.D. stainless 
tubing were installed in the walls of the tunnel. For the straight 
section where the wall is 9/16 inches thick, the tubings were cemented 
into slots milled into the walls of the tunnel. The extruded tubings 
were then sanded flush with the surface. For the curved test section 
where the walls are only 1/32 inches thick, it was found necessary to 
solder a square metal plate to the tubings. The tubings were then 
installed into holes drilled into the wall, and then secured in place 
by soldering the plate onto the wall. Again, the tubings were sanded 
flush with the surface. Thus, the resultant pressure taps were of the 
sharp edged, deep hole type. 

2 . 1.2 Geometry of Curved Test Walls 

In designing the curved wall, an original objective was to achieve 
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an equilibrium constant pressure profile characterized, however, by a 
constant value of the curvature parameter, 6 */r . Using simple flat 
plate turbulent boundary layer theory, an estimate of 6*(x) is 
obtained by assuming 


6*(x) 


.037x 



( 2 . 1 . 2 - 1 ) 


The empirical relation provides values of R(x) if the curvature 
parameter is maintained constant. For the convex wall, the constant 
is taken to be .01, while for the concave \jall, the constant is assumed 
to be .007. The results are displayed in Figures 3 and 4 for the curved 
walls respectively. 


2.2 Measurement Techniques and Instrumentations 

2.2.1 Wall Static Pressure 

A row of wall static pressure taps is provided on the tunnel 
center line in both the straight test wall and the curved test wall. 

In addition, two rows of pressure taps, 14 inches above and below the 
tunnel center line, are provided on the curved test walls. With the 
convex test wall, the pressure taps are spaced as follows. The first 
pressure tap in the tunnel center line is situated at x = 4 inches 
(see Figure 5), then every 10 inches until x = 44 inches, then every 
inch until x = 58 inches, then every 2 inches until x = 78 inches where 
the last pressure tap is situated. The off center line pressure taps 
are situated at x = 51, 57, 62, 68, 72 and 76 inches respectively. 

With the concave test wali, the pressure taps in the tunnel center line 
start at x = 6 inches (see Figure 6), then every 10 inches until x = 46 
inches, then every 4 inches until x = 5*+ inches, then every inch until 
x = 70 inches, then every 2 inches until x = 110 inches where the last 
pressure tap is situated. The off center line pressure taps start at 
x = 46 inches and are spaced at 4 inches apart until x = 110 inches. 
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The pressure taps along the tunnel center line are used for the measure- 
ment of the potential velocity at the surface while the off center line 
pressure taps are used only for checking the pressure variation in the 
transverse direction. They also serve to indicate the influence of the 
secondary flov (see Section 2.3*3)* 

The wall static pressure taps are connected to the negative port 
of a Pace Model CP510-.1 PSID pressure transducer. The positive port 
of the pressure transducer is connected to a reference pressure which 
is provided hy a pitot tube situated in the free stream near the tunnel 
entrance. The transducer output is connected to a DISA Model 55 I>30 
Digital DC Voltmeter which can be read to .001 volts. The transducer 
has a maximum output of 5 volts and this corresponds to a pressure of 
3.51 inches of alcohol (sp. gr. = 0.791)* Static calibration of the 
pressure transducer is carried out against a micro -manometer for the 
range 0 — 1 inch of alcohol, and against a manometer which 
reads to within .02 inch of alcohol for the pressure range of 1 - 4 
inches of alcohol. The calibration was checked from time to time and 
it was found that the calibration curve was quite repeatable. 

2.2.2 Yaw Measurement s 

In the wall regions the centrifugal acceleration of the low velocity 
end-wall boundary layer flows cannot balance the pressure gradient 
impressed by the main stream. The resultant end-will cross flows or 
secondary flows influence the flow near the tunnel center line and it 
may no longer be two-dimensional. Theoretically, if the aspect ratio 
of the channel is large enough, the secondary flow will be confined to 
a small region near the end walls, and the flow in the core will be 
essentially two-dimensional. The present tunnel has an aspect ratio of 
8 in the straight section and varies to a minimum of 6.2 in the curved 
test section with the convex test wall (see Figure 3)* I n order to find 
out the extent of the secondary flow, yaw measurements inside the boundary 
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layers were made at the tunnel center line and at planes above and 
below the tunnel center line. 

For the case of constant pressure flow along convex wall , yaw 
measurements were made at Stations 3, 4 , 5, 10, 11 and 12. For the 
case of separating flow along convex wall , yaw measurements were made 
at Stations 3, 4 , 5> 8 , 9 and 10. The various locations of these 
stations are shown in Figure 5* For the case of constant pressure 
flow along concave wall , yaw measurements were made on the convex wall 
opposite the concave test wall. These measurements were made at x = 70 
inches and x = 98 inches and at planes 16 inches above and below tunnel 
center line. The locations of these stations are shown in Figure 6 . 

A Conrad probe (used as a null direction probe) is used to measure 
flow directions inside the boundary layer. The particular probe geometry 
(see Figure 7) is chosen because of its high sensitivity and its zero 
scale effect (see Bryer, et.al. 1958 ). The probe has a stem of l /8 
inches O.D. stainless steel tube and a goose-neck (Figure 8 ) of .063 
inches O.D. stainless steel tube. With this shape, the edge of the 
probe is in line with the axis of the probe. Therefore, by aligning 
the axis of the probe with the radial lines at the point of measurement, 
the traverse across the boundary layer will be perpendicular to the 
surface. The Conrad probe is inserted into the rotary device (Figure 9) 
which in turn is secured onto the probe carrier. The probe carrier is 
slid into a vertical stand and secured tightly in the level where the 
probe is to be introduced into the test section. The bottom of the 
stand is clamped to the floor of the tunnel and the top of the stand 
is also clamped to one of the ribs on the tunnel wall (see Figure 10). 
The probe is then introduced into the test section. 

The introduction is through holes 7/16 inches diameter on the wall 
opposite the test wall. A plexiglass plug is machined to fit into the 
hole and flush with the wall. A hole just large enough to allow the 
probe to pass through is drilled onto the plug. The plug-probe assembly 
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is then taped onto the wall. When this prole is securely in place, it 
is traversed towards the test wall. An electrical circuit is arranged 
so that when the probe is touching the curved wall, a light goes on. 

When the probe is brought back to such a position that the light just 
goes off, the reading on the vernier mounted on the probe carrier is 
noted. This gives the zero reading for y . This way of positioning 
the probe has an accuracy in y of + .001 inches 

The probe is rotated until the pressure is the same in the two 
openings. The pressure difference between the two openings is measured 
by a Pace Model CP51D-.1 PSXD pressure transducer which is connected to 
a DISA 55D30 Digital Voltmeter, and is indicated by a zero reading on 
the voltmeter. However, due to pressure fluctuations, the voltmeter 
can be at best read to + .01 volts. The angle through which the probe 
has been rotated can be read from the dial on the rotary probe carrier 
(Figure 9). Depending on the mean flow velocity, this corresponds to 
an accuracy of + 0.5 3 in the determination of the flow direction. The 
sensitivity of the Conrad probe was determined for two different mean 
flow velocities in the fully developed pipe flow system and the result 
is given in Figure 8. Also, the null angle of the Conrad probe was 
determined in the pipe flow system. A more detailed explanation of the 
positioning of the probe, the accuracy of the yaw measurements, and the 
determination of the null angle is given in Appendix B. 

2.2.3 Velocity Profiles 

In most boundary layer measurements, the mean velocity is measured 
with a small pitot-static tube which is traversed in the direction per- 
pendicular to the wall. The basis for these measurements is the incom- 
pressible, frictionless Bernoulli equation 



The mean velocity thus measured depends on how accurate the pitot-static 
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tube measures the static pressure. The presence of this probe causes 
the streamlines to diverge as the leading edge of the probe is approached. 
This changes the curvature of the streamlines and as a result the local 
static pressure. The effect of the leading edge is felt many probe 
diameters downstream as the streamlines return to parallelism with those 
of the main flow. On the other hand, the disturbances created by the 
stem propagate upstream and produce a local variation from true static 
pressure. However, the two effects cause opposite changes in the stream 
pressure, and therefore a position can be found in the probe to locate 
the static pressure holes such that the effects of leading edge and stem 
cancel each other and again the true static pressure is measured. 

Pitot -static probes of such design are commonly used in boundary layer 
measurements on flat surfaces where the probe is essentially parallel 
to the surface. For boundary layer flows along curved surfaces, the 
probe will no longer be parallel to the surface. As a result, the mean 
flow approaching the probe is at an angle of attack, and the probe 
will no longer register the true static pressure. In view of this, 
the conventional pitot -static probe was not used for velocity profile 
measurements. Rather, a total head probe was used to measure the local 
stagnation pressure, and the velocity was calculated from the incompressible 
frictionless Bernoulli equation together with the y momentum equation 
(see Section 2.4.1). 

The stem of the total head probe is made from l/8 inches O.D. 
stainless steel tube while the goose-neck is formed from .032 inches O.D. 
stainless steel tube flattened at the tip to an outside vertical dimension 
of .008 inches with an opening vertically of .004 inches (Figure 8). 

Rogers and Berry (1950) found that the response of such a flat nose 
probe was quite independent of yaw angle of 15 degrees and less. The 
accuracy of the present probe was checked against a standard pitot -static 
probe in a fbJQy developed turbulent pipe flow system. The wall static 
pressure at the point of measurement was used as the reference for the 
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total-head prole. Assuming the static pressure to be constant across 
the pipe the probe gives the dynamic head directly. The measurements 
obtained using the pitot-static probe and the present probe were identical. 

For the convex test section, velocity profile measurements were 
made at a total of twelve stations situated on both the straight and 
curved sections. Of the twelve stations, eight were situated along the 
tunnel center line, two were situated at a plane 14 inches above tunnel 
center line and the other two were situated at a plane 14 inches below 
tunnel center line. The locations of the center line stations together 
with the locations of the off center line stations are shown in Figure 5 . 
For the case of constant pressure flow, the off center line stations 
were located directly above and below Stations 4 and 11, while for the 
separating flow case, ttey were located directly above and below Stations 
4 and 9 . This was because the flow was near separation at Station 11 for 
the separating flow case. 

The positions of the various stations were measured from the entrance 
to the straight section. Since the constant pressure flow experiment 
was run first, the locations of the various stations on the test (inner) 
wall were selected and radial lines were marked on both the top and 
bottom walls of the tunnel at these locations. The adjustable (outer) 
wall was then installed and the correct pressure distribution was set 
up on the test wall (see Section 2.3.2). The intersection between the 
radial lines and the outer wall would then give the locations of the 
corresponding measuring stations on the outer wall. The total-head 
probe is introduced into the tunnel from the outer wall and is secured 
in place in much the same way as that used to secure the Conrad probe. 

The probe is carried by the probe carrier shown in Figure 10, and is 
positioned radially much the same way as that used to position the 
Conrad probe. The probe carrier is mounted in the same manner (see 
Figure 10). Again, the same electrical circuit is used to indicate 
whether the probe is just touching the wall. Since the surface of the 
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straight section was made of white panelite, it does not provide 
electrical connection. To remedy this, a thin coat of silver paint 
was painted on the tunnel center line of both side walls of the straight 
section. 

The same outer wall was used for the separating flow experiment, 
but it was adjusted to give the correct pressure distribution on the 
test wall. Therefore, if the same holes were to be used to introduce 
the probe, their corresponding positions cn the test wall will be 
different. The locations of these stations are also shown in Figure 5 . 

It is seen that the positions of Stations 1 through ^ are not changed. 

The correct positioning of the probe radially depends on aligning the 
probe with the radial lines inscribed on both the top and bottom walls 
of the tunnel. However, with the separating flow case, there was only 
one set of radial lines inscribed on the lucite top. It was impossible 
to inscribe a corresponding set of radial lines on the bottom wall 
without taking the outer wall apart. This was not advisable once the 
correct pressure distribution has been set up. 

The outlet of the total -head probe is connected to the negative 
port of a Pace Model CP51D-*1 psid pressure transducer. The positive 
port of the transducer is connected to a reference pressure which is 
provided by a pitot tube situated in the free stream near the tunnel 
entrance. The transducer output is connected to a DISA Model 55D30 
Digital DC Voltmeter which can be read to .001 volts. The calibration 
of the transducer is described under Section 2.2.1. 

At the time the velocity profile measurements were made on the 
convex wall with constant pressure distribution, the hot-wire equipment 
was not available. When the hot-wire equipment was ready, the velocity 
profile at Station 11 for the constant pressure flow case vas made with 
the hot-wire equipment. The velocity profile calculated from the total- 
head probe measurement was then compared to the hot-wire measurement and 
the two profiles overlap each other (Figure 37)* Because of its convenience, 
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the total-head probe was also used to measure the velocity profiles 
for the separating flow case. However, the total-head probe was not 
used to measure the velocity profiles in the case of the concave test 
wall. Rather, hot-wire was used. The hot-wire technique is described 
in Section 2.2.5. 

Anticipating the fact that a system of longitudinal vortices 
would exit along the concave wall, the z positions of the locations 
where measurement was to be made could not be selected without a know- 
ledge of the vortex system. The x positions of the stations were 
chosen at 24, TO and 96 inches respectively. In order to determine the 
z positions, a hot-wire technique (see Section 2.2.6) was used to traverse 
the flow at x = 24, TO and 96 inches respectively in the z direction. 

The z positions were then taken to be the tunnel center line or the 
positions of the trough and crest of a wave. The reason for this 
choice of z positions is discussed in Section 3*3 * As a result of the 
these traverses, the positions of the five locations chosen are shown 
in Figure 6. Shown also are the z distance of these stations as 
measured from the tunnel center line. The hot-wire probe is introduced 
through the inner (adjustable) wall. The same method as before was used 
to secure the probe in place and the same technique as that used to 
position the Conrad probe radially was used to position the hot-wire 
probe. As to the y measurement, a description of the method used is 
given in Section 2.2.5. 

2.2.4 Hot-Wire Equipment and Calibration 

The hot-wire equipment wacd in the present experiment is shown in 
Figure 11. It consists of the following standard, commercially available 
units: one TSI Model 1010A Constant Temperature Anemometer, one DISA Model 
55D10 Linearizer, one DISA 55D35 True RMS Voltmeter, two DISA 55D30 Digital 
Voltmeters, one Texas Instrument X-Y Plotter, one Pace Associates X-Y Plotter 
and one Techtronic Twin Beam Oscilloscope. The sensors and the probes were 
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obtained from Thermo System Inc. The sensors are tungsten vires and 
their characteristic dimensions are .00015 inches in diameter and .05 
inches in sensitive length. Both ends of the sensor are copper 
plated. Throughout the whole experiment , factory mounted vires were 
used. The characteristic dimensions of the probe (see Figure 12) are: 
length of stem, 1.5 inches; diameter of stem, .059 inches; length of 
prongs, .25 inches; distance between tips of the prongs, .0 6 inches. 

The prongs are gold-plated. 

A small portable calibration tunnel was built for hot -wire cali- 
bration. The tunnel is of the open-return type and is powered by 
a variable speed fan. A diagram of the calibration tunnel is shown in 
Figure 13. The air enters through a smooth bell -mouth which is followed 
by a layer of honey-comb straightener and a fine mesh screen. This is 
followed by a settling chamber that is three diameters long. A con- 
vergent nozzle accelerates the flow by a 4:1 ratio, and at the same time 
provides an axisymmetric contraction to the flow. The straight calibrating 
section is located behind the nozzle, and the station used for calibration 
is situated one diameter downstream of the nozzle exit. The air velocity 
in the tunnel can be varied continuously from zero to 100 ft/sec. The 
velocity at the calibrating station is monitored by a fixed pitot probe 
located on the center line upstream of the nozzle entrance and a wall 
static hole drilled diametrically opposite the calibrating station. The 
turbulence level was measured to be .2^. when the air velocity in the 
tunnel was 90 ft /sec. 

The relation between the output voltage of a constant temperature 
anemometer and the mean flow velocity is given by the equation 

1 

2 m , 

e = A + BU (2. 2. 4-1) 

where A , B and m are constants to be determined by direct calibration. 
In arriving at this equation, free convection is neglected. Therefore, 
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in general A is not equal to the square of the zero flow voltage, 
and this complicates the calibration. The procedure finally adopted 
was to determine the constants A and B for different values of 
the exponent l/m by least-square fitting the data to equation (2. 2. 4-1). 
The chosen value of m is the one that gives the minimum deviation which 
is defined as 

1 

a 2 (m) - E ( e 2 - A(m) - B(m) U^} 2 (2. 2. 4-2) 

All calibration 
points 

2 

The data obtained showed no clear minimum for a . This means 
that the choice of m was not critical for the type of sensors used 
and the range of velocities considered. A plot of e 2 versus 
for different values of m is given in Figure 14. On the other hand, 
it was found that A/e Q where e Q is the zero flow voltage, corre- 
lates very well with m (Figure 15). As a result, l/m was chosen 
to be .42 and A/e Q was taken to be .89 according to the recommendation 
of DISA. 

A diagrammatic layout of the hot-wire equipment is shown in Figure 
l6. For the present purpose, the linearizer which is connected in series 
with the constant temperature anemometer is just considered to be an analog 
computer having a well defined transfer function. A typical calibration 
curve of the linearizer output is given in Figure 17 . 

2.2.9 Hot-Wire Techniques 

Mean Velocity Measurements 

Once the hot-wire has been calibrated, it is a simple matter to use 
it to measure the mean flow velocity. The mean value of the linearized 
hot-wire signal follows the law 

e = hQf(<t>) (2. 2. 5-1) 

L 
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where Q is the magnitude of the velocity vector, <> is the angle 
between the normal to the sensor and the direction of the mean velocity, 
and h is the calibration constant. If the cross-flow is negligible 
compared to the mean flow U , (this is true for the two cases of 
constant pressure flow along convex and concave walls ) then Q = U . 

The function f is symmetric with respect to its argument and is also 
normalized to 1 when its argument is zero, therefore by aligning the 
hot-wire so that its normal coincides with the flow direction, 

e = hU (2. 2,5-2) 

which gives the mean velocity directly. 

The straight hot-wire probe (figure 12) with a straight wire was 
introduced Into the tunnel in the same manner as the Conrad and total- 
head probes and was carried by a rotating mechanism - probe carrier com- 
bination shown in Figure 10. Due to the fragility of the wire, a 
different method was used to determine the y position of the wire. 

First, a dummy probe (it can be a hot-wire probe with broken wire) is 
used. As before, an electrical circuit is used to indicate contact 
between the probe and the wall. The reading on the vernier is noted, 
and the probe is retracted until it is completely outside of the tunnel. 
Then the dummy probe is replaced by a hot-wire probe. Knowing the exact 
measurements of the hot-wire probe, the dummy probe and an average value 
of the backlash of the probe carrier, the y position of the wire can 
be determined to an accuracy of + .008 inches. A detailed description 
of this method is given in Appendix B« 

Turbulence Measurements 

The standard technique used in the measurements of turbulent 
stresses is the employment of an X -prose or a V-probe. However, this 
method requires accurate alignment of the probe, and this poses serious 
problems in the present investigations. In addition, if the longitudinal 
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vortex system does exist along the concave test vail, the flow will no 
longer be two-dimensional and the application of the X-probe or V-probe 
will be quite doubtful. For these reasons, a rotating -wire technique 
that is a modification of a technique proposed by Fujita and Kovasznay ( 1968 ) 
is selected. The method is essentially the same as the one used by 
Bissonnette (1970), therefore the method will not be discussed in detail 
here. Instead, the reader is referred to the thesis of Bissonnette. 

The hot-wire probes used are the straight and the slanted probes 
shown in Figure 12. With the straight probes two types of wire are used; 
one is straight and the other is a k 5° slanted hot-wire. With the slanted 
probe, only the straight wire is used. 

The turbulent stresses can be determined by two sets of measurements, 
one set is given by the straight probe with the straight wire, while 
the other set is given either by the slanted probe with the straight 
wire or the straight probe with the slanted wire. For reasons that 

Be apparent later, the turbulent stresses in the case of constant 
pressure flow along convex wall were measured using straight wires on 
both the straight and slanted probes. For the other two cases, the 
turbulent stresses were measured using straight and k 1 ? slanted wires 
on the straight probe only. 

With the straight wire on either the straight probe or the slanted 
probe, the wire rotates in the xoz plane. The coordinate systems 
for both the straight wire and the slanted wire are shown in Figure l 8 . 

Both the straight probe and the slanted probe were introduced into the 
tunnel the same way as before, and the roation was chosen to span 
approximately -1*5° < <t> < • In order to minimize the backlash errors 

the recording of the linearized signal was always performed in the same 
direction of rotation. A sample recording of the mean value and the 
mean square value of the linearized signal versus the angle & of 
rotation ( straight -wire on straight probe) by means of X-Y plotters are 
shown in Figures 19 and 20 respectively. The recordings were performed 
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in both directions of rotation to account for the hysteresis due to the 
different time constants of the low pass filter and the RMS neter. The 
true signal was taken to be the average curve drawn between the two 
traces. Recording the signal in this manner 'vs.s very time consuming 
(roughly about 15 minutes per experimental point) and in order to 
reduce this experimentation time, the two time constants were chosen 
to be approximately the same. Assuming comparable hysteresis in both 
cases, recording in only one direction was required, and results obtained 
were practically the same as those obtained by recording two traces. 

With the 1+5° slanted wire on the straight probe, the wire generates 
a conical surface with axis parallel to the oy axis when the probe is 
rotated (Figure 18). In order to avoid the disturbance due to the longer 
prong, the probe is positioned in such a manner that the longer prong 
is always behind or parallel to the shorter prong during rotation. 

For the case of constant pressure flow along convex wall, turbu- 
lence measurements were made at Stations 1, 7, 9 and 11, while for the 
case of separating flow along convex wall, turbulence measurements were 
made at Stations 1, 4, 9; H and 12. For the case of constant pressure 
flow along concave wall , turbulence measurements were made at all 
five stations. 

2.2.6 Mapping of the Longitudinal Vortex System 

In order to make meaningful measurements on the concave wall, the 
question of the existence of the longitudinal vortex system has to be 
resolved. From the data of Tani (1962) it can be seen that the vortex 
system is very likely to be stationary, otherwise Tani could not have 
detected it with a pitot traverse. For the present investigation, a 
hot-wire set at a constant distance away from the concave wall is used 
to traverse the flow in the z direction. The mechanism used to carry 
the hot-wire probe and drive it in the z direction is shown in Figure 21. 

Briefly, the probe adaptor used in conjunction with the TSI miniature 
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hot-vire probe is firmly secured to a movable block which slides within 
two groves on the strut of the probe carrier (Figure 22). The leads 
of the adaptor are connected to a co-axial wire which has a BNC connector 
on its other end. The straight hot-wire probe can then be plugged into 
the probe adaptor. A constant speed motor is used to drive the movable 
block along the strut. The speed of the motor can be varied by regulating 
the supplied voltage. The linear movement of the probe is translated 
into an electrical signal through a ten turn pot, so that its movement 
can be recorded on a X-Y plotter. The probe carrier is so constructed 
that the traverse in the z direction has a span of 21 inches, and 
is centrally located with respect to the tunnel center line. At both 
ends of the traverse, there is a limit-switch. When the movable block 
is located at either end, it triggered the switch and the motor is turned 
off. The pot, the motor and the gear assembly are located on top of the 
strut, and a piece of £ inches thick foam rubber is cemented on the 
bottom of the strut. The length of the strut including the piece of 
rubber is a bit longer than 4 feet. 

phe straight hot-wire probe in position, the strut was intro- 
duced into the tunnel from the top. It was then fastened onto the 
lucite top by screws as shown in Figure 23. The three rectangular 
holes on the top wall of the tunnel at x = 24, JO and 96 inches respec- 
tively were located in such manner that when there were no spacers in 
front of the probe carrier assembly, the hot-wire was exactly .2 
inches away from the wall. By moving the spacers to the front, the 
hot-wire can be set at y = .4, .6, .8 and 1 inches away from the wall. 

The straight hot-wire was positioned to give the maximum voltage 
output. The linearized signal was connected to one arm of the X-Y 
plotter. The X axis of the X-Y plotter was calibrated to indicate 
the z position of the probe. The origin of z is taken to be the 
tunnel center line, and z is positive when measured upward and negative 
downwards. The traverse was always in the same direction, i.e. from top 
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to bottom, and the speed of the motor vas regulated to give a linear 
speed of approximately 6 in/min. 

Since the purpose of this investigation was to detect and map 
out the longitudinal vortex system, no quantitative measurements of 
the vortices were made. Once the vortex system was mapped out, the 
locations of the stations where boundary layer measurements were to be 
made could be determined. 

2.3 Qualification of Curved Wall Tunnel 

2 . 3 .I Two-Dimensional Equilibrium Turbulent Flow in Straight Section 

Uniformity and steadiness of the flow in the straight section of 
the tunnel was checked by observing tufts attached to the walls of the 
tunnel. The equilibrium nature and two -dimensionality of the flow were 
checked by actual measurements and comparison with Klebanof^s (1955) 
data under similar circumstances. 

First, the boundary layer was checked to see if it was turbulent. 

To do this, velocity profile measurements were made at Station 1 and 
x = 40 inches. Results indicated that the flow was still laminar up to 
Station 1 and the boundary layer thickness at x = 40 inches was much 
less than 1 inch. The flow was then tripped a small distance downstream 
of the entrance to the straight section. The tripping was effected by 
a slightly stretched piano wire covered by electrical tape. The 
diameter of the wire was chosen to give a boundary layer thickness at 
x = 40 inches of approximately 1 inch. For the sake of symmetry, the 
flows on the other three walls were also tripped. With this arrangement, 
the flow was found to be turbulent at Station 1, and that the boundary 
layer thickness vas approximately .6 inches. The total pressure profile 
at Station 1 was then measured using the flattened total-head probe and 
the velocity profile was calculated together with the skin friction at 
the wall (see Section 2.4.2). The defect plot of the velocity profile 
was then compared with the zero pressure gradient profile of Mellor 
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and Gibson (1966). The comparison is excellent and the result is 
shown in Figure 24. The calculated profile is given by the solid line 
in the figure. The result indicates that an equilibrium turbulent flow 
exists at Station 1. 

The rotating -wire technique was used to measure the Reynolds 
stresses. First, the straight probe with the straight wire was used 
and this gives the following three Reynolds stresses, namely u 2 , uw 
and w 2 (see Section 2.4.3 for turbulence data reduction). Then the 
slanted probe with the straight wire was used. In order to get uv 
and v 2 with this probe, vw has to be either known or assumed zero. 

In light of the fact that for two-dimensional flow, due to symmetry 
about the center line plane, the off-diagonal elements of the Reynolds 
stress tensor, except uv , are necessary zero. Therefore, the assump- 
tion of vw = 0 was made so that uv and v 2 could be calculated. 

In addition, if the measured uw were indeed very small it would provide 
an independent check for the assumption, while at the same time, it would 
q]_sq serve to show the two-dimensional nature of the flow. To see if 
this is true, the measured uw was plotted against the measured uv 
at Station 1, and is given in Figure 25. It can be seen that nowhere 
inside the boundary layer was (uw/uv) greater than .05 • This shows 
that the flow was relatively two-dimensional. To further substantiate 

this claim, the measured uv/u 2 x , u'/Upv, , v 'Aj pw > w '/ u pw and UY / U,V ' 
were plotted against y/& and compared with the flat plate data of 
Klebanoff (1955). These are shown in Figures 2 6 to 30. In general, the 
measured data compare favorably with Klebanoff' s data except near the 
wall. This is because the rotating -wire technique is not very accurate 
in region of high shear (see Bissonnette 1970). From these measurements, 
it can be concluded that the flow was in equilibrium and was indeed two- 
dimensional . 

After the first set of experiments with the convex wall was finished, 
the pressure distribution on the test section was set up to give a 
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separating flow in the curved tunnel. In order to make sure that the 
same entrance flow exists in this case as in the constant pressure 
flow case , measurements at Station 1 were repeated again. Surprisingly, 
the measured velocity profile did not show any noticeable difference, 
hut the measured uw was found to he much larger than before (at one 
point uw/uv even amounts to ~.2). Also, the measured uv/u^ was 
noticeably smaller than the previous data. After a long and tedious 
check on the equipment and the rotating-wire technique (including the 
test of the technique in the fully developed pipe flow) it was found 
that the flow was not transisting properly. The improper transition of 
the flow was caused by the loosening of the tripping wire. As a result, 
the transition was not uniform along the walls, thus creating a certain 
skewness in the flow, and the off-diagonal elements of the Reynolds 
stress tensor was no longer zero. To correct for this, a tripping 
device proposed by Hama (1957) was used. 

Briefly, the tripping device consists of isoceles tr ian gles made 
of either electrical tap or aluminum sheet. These triangles are then 
cemented to the walls of the tunnel so that the vertices of the triangles 
are facing the flow. According to Hama, the device was very efficient, 
and anchored the transition right at the tripping device. With this 
modification, the measured uw/uv was again found to be less than .08 
across the boundary layer. However, the boundary layer at Station 1 
was found to be approximately 20$ thicker than the previous measurement. 
The measured velocity profile together with the various Reynolds stresses 
at Station 1 for the case of separating flow along convex wall are shown 
in Figures 24, and 26 to 30. Again the flow was in equilibrium and was 
two-dimensional . 

When the test wall was changed to the concave wall, the same 
measurements were performed at Station 1. However, in this case, the 
wall opposite the previous test wall was the test wall. The results 
of these measurements are again plotted in Figures 24, and 26 to 30, and 
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they show that the flow was in equilibrium and two-dimensional . In 
addition, they show that the flow in the straight section was very 
uniform. 

2.3.2 Pressure Distribution on Curved Test Section 

In the course of setting up the wall static pressure distribution 
on the convex test wall it was found that the pressure decreased 
slowly as the flow approached the entrance to the curved section, and 
then steeply as the flow entered the curved section. The pressure drop 
amounts to more than 6 0$. of the reference dynamic head. An effort 
was made to reduce this amount of pressure drop by increasing the 
cross-sectional area right at the entrance, but to no avail. Later 
it was found that the flow was separated on the adjustable wall and 
right at the entrance to the curved test section. Therefore, increasing 
the cross-section area at this section merely increased the separated 
region, and the flow on the test wall was undisturbed. To correct this, 
a side wall jet was installed on the adjustable wall and right at the 
entrance to the curved section. The jet is merely an opneing, four 
feet high, on the side wall of the tunnel. The pressure difference 
between the outside and inside of the tunnel provides sufficient 
momentum to the jet. The opening of the jet is controlled by a flexible 
flap attached to the straight wall. The flap extends 8 inches into the 
straight section, thus providing some control on the flow approaching 
the curved section. The jet and flap control mechanism are shown in 
Figure 31 . With this arrangement, it was found possible to reduce the 
pressure drop to about 30$. of the reference dynamic head. However, 
the pressure drop cannot be eliminated completely. Because of this, 
the final constant pressure distribution set up on the convex wall was 
at a different level from that of the straight wall. Immediately after 
the pressure drop, there was a slight pressure increase. All these 
occurred within a distance of 6 inches, from x = hG inches to 52 inches. 



After that the adjustable vail vas adjusted to give the minimum possible 
variation of vail static pressure on the convex vail. The above adjust- 
ment vas done vith no regard to the secondary flow in the tunnel. 
Therefore, it vas no surprise to find that the secondary flovs measured 
at Stations 3 and 5 vere quite appreciable. Means vere introduced to 
control the secondary flov (this is discussed in Section 2.3.3) and 
this changed the pressure distribution. The secondary flows were 
adjusted to a minimum and the adjustable wall vas again set to give 
the best pressure distribution. The final pressure distribution obtained 
for the case of constant presuure flow along convex vail is shown in 
Figure 32. The final setting of the adjustable wall is given in Figure 3. 
Finally, the wall static pressure at planes 14 inches above and below the 
tunnel center line were also measured. If these measurements were plotted 
on Figure 32, they would overlap the tunnel center line measurements. 
Therefore, for the sake of clarity, the off-center line wall static 
pressure measurements are not shown. 

In setting up the pressure distribution on the convex wall for 
the separating flow case, the objective vas to have a linear decelerating 
potential velocity at the wall (Figure 82). By pulling back the adjustable 
this was relatively easy to set up in the curved section between 
x = 50 inches to 60 inches. After x = 60 inches, the pressure distri- 
bution began to level off. Again, the adjustable wall was pulled back, 
but there was relatively little change in the pressure distribution at 
the test wall. At this point, the flow was suspected to have separated, 
but the line of separation was not known. In an effort to locate the 
separation line, tufts were attached to the convex -wall between x = 60 
inches and 70 inches at an equal spacing of two inches apart. The effort 
proved to be futile, but it did seem to indicate that separation did not 
occur along a straight line from the top to the bottom of the tunnel. 

There vas clear indication that the flow separated first near the top 
and bottom of the tunnel and further downstream along the central 
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portion of the tunnel . 

Attempts were made to move the separation curve further downstream 
by readjusting the whole adjustable wall. The results was a non-linear 
distribution of decelerating potential velocity at the convex wall. 

After much trial and error, the original pressure distribution was set 
up. This is given in Figure 3 2. The final setting of the adjustable 
wall is shown in Figure 3. The linear distribution of decelerating 
potential velocity at the wall is displayed in Figure 8l, together with 
the actual measured potential velocity at the wall. As far as can be 
made out from the tufts, separation occured somewhere between x = 66 
inches and TO inches. Later measurements of the velocity profiles lend 
evidence to this. The secondary flows at Stations 3, 5, 8 and 10 were 
tten measured, and the end-wall jets were adjusted to give the minimum 
secondary flow at these stations. The pressure distribution was mea sirred 
again, but in this case no repeated adjustments were made as in the case 
of the constant pressure flow. The final results are those mentioned 
above. The secondary flow measurements are discussed in detail in the 
next section. 

The wall static pressure at planes 14 inches above and below the 
tunnel center line was measured. The difference between these measure- 
ments and the tunnel center line measurements were very small at the 
leading section of the convex wall, and became noticeable near separation. 
However, for the sake of clarity, they are not shown in Figure 32. 

For the concave test wall, the pressure at the wall increases as 
the flow enters the curved section. In order to reduce this pressure 
increase, the flow in the region upstream and downstream of the entrance 
to the curved section was accelerated by adjusting the adjustable wall. 

As a result, no side wall jet was necessary. The adjustable wall extends 
one foot into the straight section, and this provides sufficient 
adjustment for the region upstream of the entrance to the curved section. 
Again, after much trial and error, the final constant pressure distribution 



set up on the concave vail is at a level about 22 $ higher than at the 
straight section. The pressure distribution is given in Figure 33, 
while the final setting of the adjustable wall is shown in Figure 4. 

In this case, the variation of the wall static pressure in the stream- 
wise direction is greater than the corresponding case with the convex 
test wall. Efforts were made to reduce these variations but to no avail. 
The wall static pressure at planes 14- inches above and below the plane 
of symmetry was also measured. Again, there were noticeable differences 
between these measurements and the tunnel center line measurements. To 
show the extent of these variations, the off center line wall static 
pressure measurements are also shown in Figure 33- 

As in the case of the constant pressure flow along convex wall, 
the secondary flows at planes l 6 inches above and below the tunnel center 
line were measured, and the end wall jets were adjusted to give the 
minimum secondary flow. The pressure distribution was then re-set and 
the secondary flow was again measured and readjusted if necessary. The 
whole procedure was repeated until a satisfactory pressure distribution 
was obtained together with a reasonable secondary flow at the stations 
where measurements were made. The secondary flow measurements are 
discussed in detail in the next section. 

2 . 3.3 Secondary Flows in Curved Test Sections 

As has been explained in Section 2.2.2, secondary flows exist near 
the end walls of the tunnel because of the longitudinal curvature of 
the test section. The influence of this end-wall flow can extend to 
the tunnel center line. 

When the tunnel was first designed, it was thought that if it has 
a large aspect ratio, the secondary flows would be essentially confined 
to small regions near the end walls. However, initial yaw measurements 
at Stations 3, 5, 10 and 12 on the convex wall, with the constant pressure 
distribution along it, revealed that the measured yaw angles at Stations 3 , 
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and 5 were approximately equal and were opposite in sign but they 
reached a maximum of ~ 10° . The maximum reached at Stations 10 and 
12 was about 2 (f . To reduce the secondary flows, end wall jets were 
installed on the top and bottom walls of the tunnel and right at the 
entrance to the curved section (Figure 31). A series of four tangential 
jets were installed on each end wall. The jets were made of l/8 inches 
I.D. tube, and were supplied by a high enough pressure to give critical 
flow at the jet exits. Thus, the amount of mass flow added to the 
tunnel flow is quite small while the amount of momentum added is enough 
to compensate the momentum defect of the end wall boundary layers at the 
position of the jets. This additional momentum is sufficient to reduce 
the secondary flows considerably even at Stations 10 and 12. 

With the end wall jets installed, yaw measurements were again made 
at Station 3, 5, 10 and 12. The maximum yaw angles at Stations 3 and 5 
and 10 and 12 were found to be less than 2° and 4° respectively. By 
monitoring the yaw angles at Stations 3 and 5, repeated measurements of 
the pa’ es sure distribution along the convex wall and continually adjusting 
the adjustable wall to give the minimum variation in the streamwise 
pressure, it was possible to set up the best constant pressure distribu- 
tion on the convex wall. The resultant yaw measurements at Stations 3, 
b j, 10, 11 and 12 are shown in Figure 34. The flow in the tunnel was 
relatively two-dimensional. Later velocity profiles and turbulence 

r .eusurements lend support to this claim. 

Essentially the same end wall jets operating at the same conditions 
were used to control the secondary flows in the case of separating flow 
along convex wall. In this case, because of the added complexity m 
setting up the correct pressure distribution, the procedure followed 
was different from that adopted in the case of constant pressure flow alon 
convex wall. The correct pressure distribution was set up with the end 
wall jets operating at the optimal conditions. Then, yaw measurements 


were made at Stations 3, 4, 5, 8, 9 and 10. The results are shown in 
Figure 35 , and for the sake of clarity, the yaw measurements at Stations 
4 and 9 are not shown. The maximum yaw measured at Station 4 was less 
than 1° and that at Station 9 was less than 2° . Although the flow along 
the tunnel center line is quite symmetrical, the secondary flows at 
Stations 8 and 10 are considered large. Attempts were made to reduce 
this by installing larger jets to the end walls. This managed to 
reduce the secondary flows. However, the pressure distribution along 
the convex wall was disturbed, and it was not possible to obtain a 
linear decelerating potential velocity at the wall. After much trial, 
it was finally decided to go back to the original set up. Later 
velocity profile and off center line wall static pressure measurements 
seem to indicate that the flow was relatively two-dimensional up to 
x ~ 60 inches. Downstream of x = 60 inches, the influence of the 
secondary flow can no longer be discounted. 

For the case with the concave test wall, secondary flow measurements 
were taken on the adjustable wall opposite the concave test wall. The 
reason for doing this is because of the existence of the system of 
longitudinal vortices. Under the influence of the longitudinal vortices, 
the secondary flow measurements would depend to a large extent on where 
the measurements were taken. As a result, the true influence of the 
end-wall secondary flows cannot be estimated. 

From previous experiences, yaw measurements at tunnel center 
line were found to be unnecessary. For this reason yaw measurements 
were made at planes l6 inches above and below tunnel center line at 
x - 70 and 96 inches for the case of constant pressure flow along concave 
wall. The same end wall jets were used to limit the secondary flows to 
small regions near the end walls. The pressure distribution was set up 
on the concave wall, and yaw measurements were made at the off center 
line stations at x = 70 inches. The supplied pressure of the jets 
was then adjusted to give minimum yaw at the two stations where 
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measurements were made. The pressure distribution was measured again, 
and if considerable changes were noticed, the above procedure was repeated 
all over until a satisfactory pressure distribution on the concave test 
wall and a relatively small yaw profile across the boundary layer on 
the adjustable (convex) wall were obtained. The final yaw measurements 
at x = TO and 9 6 inches are given in Figure 36. 


2.4 Data Reduction 

2.4.1 Calculation of Velocity Profiles 

Assuming = pkU 2 , Bernoulli's equation may be written as 
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where the reference pressure P^ is here taken to be the total pressure 
in the potential core so that P^ - P r -* 0 outside of the boundary 
layer. Therefore the potential velocity is given by 
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Within the boundary layer ky < 0.1 and it is possible to simply 
represent the potential velocity according to U p d U^l-ky) (see Section 
4.2). Furthermore, the integrals in equations (2. 4.1-1) and (2.4.1 2), 
which represent the difference between the static pressure calculated 
from the actual velocity and the potential velocity, are small. Therefore 
it is possible to write 

F t ~ ?r . + e ~ 2k y (2. 4. 1-3) 
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In Figure 37 the result of using equation (2. 4. 1-3) to reduce 
the data is compared with that of equation ( 2 . 4 . 1 - 1 ) which is also 
normalized with equation ( 2 . 4. 1 - 5 ) instead of equation ( 2 . 4 . 1 - 2 ) after 
setting y = 0 . The latter would be the consistent normalization if 
equation ( 2 . 4. 1 - 1 ) were used, but then no distinction would be visible 
in Figure 37- 

As discussed in Section 4.2, the present procedure is consistent 

with the boundary layer equations used to compute theoretical profiles. 

It should be noted that throughout this investigation ky < 0.1 and 
-ky 

e ~ 1 - ky . However, it appears exceptionally convenient to retain 
the complete exponential form. 


2.4.2 Skin Friction Deduction 

The skin friction at the wall was not measured, hut rather it was 
obtained from the Clauser plot (1956). For curved flow, the skin 
friction is defined as 
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( 2 . 4 . 2 - 1 ) 


and since the Clauser plot is independent of the reference velocity 

chosen to make the velocity profile non-dimensional, the reference 

velocity is here chosen to be U ; The velocity profile was calculated 

PW yU 

and u/lJp^ was plotted aginst log . On top of this was superimposed 
several plots of the Law of the Wall for different values of * The 

line that passes through the most numbers of experimental points is 
taken to give the correct for the measured profile. The correct 

C f becomes a bit difficult to decide for the near separation profiles. 

Two sample plots of the velocity profiles for the determination of C 
are shown in Figure 38 . 



2. In 3 Calculation of the Reynolds Stresses 

The General equations for analyzing the hot-wire signal verc 
derived by Bicsonnettc (1970). However, they arc recorded here for 
reference. These equations are specialized for the use of a linear! 
w ith cal ibration constant h • 
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e is the angle between oz axis and the normal to the wire, 
<l> , 5 and p are defined in Figure l8. 
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The quantities with an asterisk refer to flow quantities defined 
in a frame of reference attached to the wire (wire coordinate system) 
while those without an asterisk refer to flow quantities defined in a 
frame of reference attached to the tunnel (tunnel coordinate system). 

For the two cases where the test wall was convex, all measurements 
were made in the tunnel center line, therefore the cross-flow was very 
small (except near separation in the separating case) and p ~ o . with 
the straight wire on the straight probe, e = */2 and the wire coordinates 
coincide with the tunnel coordinates. As a result, equation (2. 4. 3 -4) 
is much simpler and ? = u 2 * , 7 = 7* , ^ ^ , etc. With the 

straight wire on the slanted probe, q ~ 46° and the turbulent stresses 
are given by equation ( 2 . 4 . 3 - 4 ). 

For the case of constant pressure flow along concave walls, measure- 
ments were taicen at the position of the crest and trough of the wave 
system and it was argued that the cross-flows at those two positions were 



very small, (for reason why refer to Section 3*3) hence p ~ 0 . This 
again allows the equations for the signal analysis to be further reduced. 

The x coordinate of each plot is calibrated to give the angle of 
rotation 6 . Once 6 and p are known, <t> is given by either 
equations (2. 4.3-1) or (2. 4. 3-5)- The function f(<t>) is normalized by 

the maximum on the trace, i.e. e(0) and is therefore given by 
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For the sake of simplicity, the functional form proposed by 
Champagne et. al. (1967) for f( ) is used, i.e. 
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where k is a constant, and is determined for every run by fitting 

2 

equation (2.4.3-10) to the experimental data. Taking k to be the 
mean value over all data points, it is given by 
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With the function f(4>) defined, the turbulence stresses u , w and 
uw* can be obtained by least square fitting equation (2. 4. 3 -3) to 

the experimental record of e T ^ . These values are then substituted into 

L * *”2* * 

equation (2. 4.3-7) &nd the other Reynolds stresses uv , v and vw 

are computed by a similar technique. 

As explained earlier, the cross -flows are small for all cases of 

measurements and hence no corrections are necessary. Hence, in adapting 

Bissonnette T s (1970) data reduction program for the present purpose 0 

was set equal to zero. The accuracy of the rotating -wire technique was 

tested in a fully developed turbulent pipe flow experiment. From the 

measured axial pressure drop, the shear stress distribution across the 
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he measured shear 


pipe can be calculated and this is used to check t 
stress profile. In general, the accuracy of the rotating-wirc technique 
is very good except near the wall. This experiment is described in 
detail in Appendix C. 
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Chapter III 
DISCUSSION OF RESULTS 


In the following, the discussion is divided into three sections. 

The first deals with the results of the constant pressure flow along 
convex wall, the second concentrates on the results of separating flow 
along convex wall, while the third discusses the results of the constant 
pressure flow along concave wall. The discussion taKcs the following 
format. The mean flow data are analysed first and this is then followed 
hy a discussion of the turbulence data. The results of these experiments 
jipe tabulated and. are given in Tables 1 to 35* Tables 1 to 12 contain 
the results of the moan flow and turbulence? measurements of the constant 
pressure flow along convex vail experiment. Tables 13 to 25 contain the 
data of the mean flow and turbulence measurements of the separating 
flow along convex wall experiment , and the mean flow and turbulence 
data of the constant pressure flow along concave wall experiment are 
given in Tables 26 to 35. For all three experiments, the free stream 
turbulence level in the straight section is less than 
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that the flow was very nearly two-dimensional were the small yaw angles 
measured at Stations 3,> ^ > 5? 10* H and 12 (Figure 30* and the 
apparently identical pressure distributions at tunnel center line and 
planes above and below the tunnel center line. Besides, the velocity 
profiles measured at Stations 3 and 5 and 10 and 12 were no different 
from those measured at Stations 4 and 11 respectively (Figure 39). To 
further chock the two-dimensionality oi the flow, the von e.n man 
momentum integral was integrated according to the method proposed by 
Coles ( 1\68 ) . 
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The von Karman momentum integral is derived in Section 4.2, and 
is given by equation (4.2-22). It can be written as 
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where the subscript o denotes measurements in a reference station. To 
integrate this equation, values of 5 , 0 and u^ at the various stations 

are required. The values of 5 and Q are obtained by direct numerical 
integration of the measured profiles, which are shown in Figure 39, 
while the skin friction is obtained from the Clauser plot. The values 
thus obtained are plotted verses x in Figure 79, together with the 
shape factor (h) development. From these measured values, equation 
(3*1-1) is integrated numerically and the result is shown in Figure 40. 

It can be seen that all along the convex test wall PL ( it denotes the 
quantity on the left hand side of equation (3-1-1)) is not too much 
d iff eren i from PR (it denotes the quantities on the right hand side of 
equation 3*l“l))> hence the two-dimensional momentum integral is satisfied 
and the flow can be said to be relatively two-dimensional. Further evi- 
dence of two-dimensionality is provided by subsequent turbulence measure- 
ments, however, this will be discussed later. 



Since the skin friction is not measured independently, the momentum 
balance calculation depends to a large extent on the validity of the 
Clauser plot to determine C for curved flow. To demonstrate the 
validity of this method it is necessary to show that the Law of the 
Wall holds even for turbulent flow along convex surfaces . The 
measured velocity profiles were plotted to chow the Law of the Wall 
region using the u determined from Clauser plot. If the Law of the 

T 

Wall indeed holds for curved flows as well as for plane flows, then a 
straight line having a slope of 5.6 can be drawn through all the data 
points in the wall law region. The result is shown in Figure 4l, where 
it can be seen that the measured profiles all show the existence of 
the wall law region. Also, all the profiles shown start to deviate 
at about the same point (yu^/v ~ 200) where the flat plate profile 
(Station l) begins to deviate from the Law of the Wall. Thus, it can 
be seen that the Law of the Wall, which is given by 


— = 5.6 log — - + B (3-1- 2 ) 

u v 

T 

also holds for flow along convex surfaces. However, this should come as 
no surprise , because in this region the mean flow streamlines are 
essentially prallel to the surface. 

By comparison with plane flow, the Defect Law for convexly curved 

flow can "be written ac 
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If the defect region overlaps with the wall law region, then the function 



F(y/A) can be shown to be 


F(^) ~ - 5-6 log | + A 


(3-1-4) 


Mellor and Gibson (19 66 ) have pointed out that A is not truly a 

constant. Instead, it Is a function of the equilibrium parameter 

p = — — ; and for p = 0 , which is the flat plate flow, A ~ - 0.6 . 

T w 

Since the Law of the Wall also holds for constant pressure flow along 
convex surface, therefore, the Defect Law as given by equation ( 3 .I- 3 ) 
should also hold. However, A should also be a function of some 
curvature parameter. A semi-log plot of (Up-U)/u T verses y /a will 
indeed show that A depends on curvature . This is displayed in 
Figure 42. Therefore, it can be seen that A = A(&/r ,p).* 

The assumption of an overlap region means that U can be described 
equally well by the Law of the Wall (3.1-2) or by the Defect Law ( 3 .I- 3 ). 
These equations my be added to give the skin friction equation 


(f) 

Of 


i u. 


pw 


— 5*6 log R~ 77 + B + A 


(3-1-5) 


where Aoq = U 5 
pw 

R_. is the Reynolds number based on 6^ 


is substituted to give the final equation and 
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In obtaining equation (3*1-5); it is assumed that u ~ U as 

P P v; 

y/A 0 . The skin friction equation thus deduced is quite general. 
It applies to all types of boundary layer flows over plane or convex 


The reason for suggesting 5 /r is because for the present experiment 
8 /R is constant. Besides, 5 remains relatively constant along 
the convex vail, thus, rendering 5 /r a very natural curvature 
parameter to use. 
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surfaces so long as A(6/R,p) is known. 

For convex curvature, A in positive, therefore, according to 
equation (3.1-5), C f will decrease faster than the corresponding flat 
plate flow. This means that the flow is less " turbulent -lixe , and as 
a result, the flow cannot sustain as high a pressure gradient as the 
the corresponding flat plate flow. In other words, under the action 
of the same adverse pressure gradient, the flow over a convex surface 
will separate first. This is a direct consequence of the reduced mixing 
activities between fluid layers, which can be seen from the following 
explanation . 

For flows over convex surfaces, the centrifugal force on a fluid 

particle must be balanced by an inward pressure gradient. If a particle 

is moving too slowly, its centrifugal force is too small, and it moves 

inward. Hence , it can be seen that the fluid particles that move away 

from the surface carry with them the velocities at the points where they 

* 

come from will have a smaller centrifugal force at their new locations . 
As a result, they will be pushed inwards by the pressure gradient and 
hence the interchange of momentum, energy, etc., between the faster 
and the more slowly moving particles are hindered. Thus mixing of 
momentum in the boundary layer will be reduced. This is evident from 
the measurements. At Station 1, the boundary layer thickness, 5 , is 
about 0.55 inches, and 5 grows to about 1 inch at the entrance to 
the curved section. The boundary layer thickness at Station 2 is about 
1 inch, and it remains constant up to Station 11 (Figure 39)- Thus, it 
can be seen that the boundary layer has not grown at all under the 
influence of convex curvature • 

By comparison with the flat plate equilibrium boundary layer, 


*This is Prandtl ' s mixing length argument which says that the linear 
momentum of the fluid particle is conserved when it is displaced from 
layer y to a new location at layer y ± dy . 
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of the boundary layer. Therefore, a plot of — — verses % would 

Uj A 


show whether the present flow has reached an equilibrium state or not. 
Such a plot is given in Figure 1+3 . From this plot, it appears that the 
flow has reached equilibrium at Station 11. However, this condition 
is not sufficient to demonstrate the equilibrium nature of the flow. 
Further evidence should be obtained from the turbulence measurements. 


Thin point will be discussed again in the following. 

The decrease in mixing activities is also evident from the 
turbulence measurements. The distribution of turbulence intensities 
and turbulence energy are shown in Figures 44, 45 and 46 respectively. 
It can be seen that there are significant decreases in turbulence 


intensities across the boundary layer. Of course, the results are 
influenced by the favourable pressure gradient at the entrance to the 
curved test section. 

After the favourable pressure gradient, which only extends to 
Station 2, the only external force that acts on the fluid particles is 
the centrifugal force created by the convex curvature. In the absence 
of wall curvature, the boundary layer would recover from the 
favourable pressure gradient and flat plate equilibrium would again be 
reached at about Station 9 (later calculations lend evidence to this). 
This means that mixing between fluid layers and turbulent diffusion 
would bring the level of turbulence right back to the flat plate values*. 


However, under the influence of convex curvature, the turbulence 
intensities are prevented from recovering. This is clearly evident from 
the measurements at Stations 7, 9 and 11 (Figures 44, 45 and 46). 

Convex curvature acts to prevent mixing and diffusion of fluid particles 
outward, therefore, after the decrease in turbulence activities caused 


As evidence by the fact that the turbulence measurements at Station 1 
correlate well with Klebanoff 's (1955) data. This shows that the distri- 
butions of u , /u 1 , v'/Uq and v' /u 1 are the same for all equilibrium 
flat plate boundary layers. 



by the favourable pressure gradient, the centrifugal force prevents any 
further increase in turbulence intensities. This is why the turbulence 
measurements at stations 7, 9 and 11 are quite similar to each other. 

It should be pointed out that there are very little turbulence activity 
in the outer part of the boundary layer, i.e. y/s > 0.4 . Instead of 
increasing from Stations 7 to 11, the turbulence intensities decrease . 

On the other hand, the turbulence intensities near the wall increase 
slightly from Stations 7 to 11. This indicates that the effect of 
curvature is much greater in the outer part of the boundary layer than 
near the wall. On the other hand, the similarity of the distribution 
of the turbulence intensities at Stations 7 > 9 and 11 seems to indicate 
that the flow has reached a new equilibrium state. If this is so, then 
the shear stress profiles at Stations 7, 9 and 11 should also be similar 
From examination of the measured shear stress profiles in Figure 47 it 
appears that the flow has not reached an equilibrium state. 

The flow, after recovering from the initial favourable pressure 
gradient, settles into a kind of quasi-equilibrium state where an inner 
region of rotational mean flow and non-zero Reynolds stress (y/& <0.4 in 
Figure 47) is embedded in a larger region (y/b < l) of rotational mean 
flow but zero Reynolds stress. Presumably the inner region would 
slowly grow until it coincides with the mean rotational region. Only 
then will it be an equilibrium flow in the conventional sense. In the 
outer part of the boundary layer, it is clearly evident that the 
"curvature effects" counteract the usual shear stress producing mechanism. 

Since the shear stress is not measured by a conventional method, 
there is always doubt as to the reliability of the data, especially the 
"inviscid like" nature of the flow in the outer part of the boundary 
layer. For the present case, a slanted probe with a straight wire is 
used to measure uv , v^~ and vw ■ It was pointed out in Section 2.4*3 
that in reducing the data for uv and v 2 , vw is assumed zero. This 
assumption is justified because of the relative two -dimensionality of 
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the flow • As an independent check, the measured uw (which can be 

obtained directly by rotating the straight probe with a straight wire 

in the u and w plane) is found to be very nearly zero. A plot of 

uw and uv measured at Stations 1 and 11 is given in Figure 25. It 

can be seen that nowhere is uw greater than % of uv . This not 

only serves to show the true two-dimensional nature of the flow in the 

curved test section, but also the reliability of the measured shear stress 

profiles and the fact that the stress vanishes when y/b > 0.4 . 

The distribution of shear correlation coefficient is given in 
Figure 48, For equilibrium flat plate boundary lyaer, the shear corre- 
lation coefficient is constant ( . 5 ) for a greater part of the boundary 
layer, but drops to zero rapidly near the edge of the boundary layer 
(measurements at Station l). For flow over convex surfaces, the shear 
correlation coefficient also remains constant (the constant varies from 
Station to Station) for a greater part of the "shear thickness", but 
drops to zero rapidly near the edge of the shear stress profile. The 
constant reached by the shear correlation coefficient is .38 at Station 9 
and .45 at Station 11. 

The turbulent energy equation can be obtained by adding equations 
(E-l4), (E-15) and (E-l6). The result is 
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From this equation and the three fluctuation components given by equations 
(E-l4), (E-15) and (E-l 6 ) or equations (E-l), (E-2) and (E-3), the following 


* 

For two-dimensional flow, because of symmetry about the center line 
plane, uw and vw are necessary zero. 



points can be noted. First, turbulent energy is not produced only in 
u 2 as in the case of a plane flow but also in v 2 . Second, the transfer 
of energy between u 2 and v 2 is due both to the pressure fluctuation 


and also to the term 


ku^v 

1+ky 


appearing in equations (E-l) and (E-2), 


Third, the total energy production on all fluctuation components is 
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oy 


A plot of the energy production is given in Figure 49- The smaller 
turbulent energy production explains the smaller turbulence intensities. 
Finally, the influence of the v 2 - production term, (+2k uv U) in 
equation (E— 2) or (E-15) on the flow should be pointed out. For flows 
over convex surfaces, k is positive, therefore 2k uv U is positive 
and this means a suppression of v 2 production. Hence, the radial 
movement of a fluid particle is suppressed. As a result, mixing 
activities between fluid layer are reduced. This means that radial dis- 
placement of a fluid particle in flows over convex surfaces is stable. 
Eskinazi and Yeh (1956) were the first to point out that this is in 
agreement with the stability criterion put forward by Rayleigh. 

The advection and the production by normal stresses are also cal- 
culated and shown in Figures 50 and 51 respectively. In flat plate boundary 
layer the advection is very much smaller than the production except in 
the outermost part of the layer. The present result shows that in the 
inner part of the layer (y /b < 0.4), the advection is one order of magni- 
tude smaller than the production (Figure 50). As to the outer part of 
the layer, curvature counteracts to reduce turbulence. The result is 
that the advection is nearly zero in this part of the layer. The production 
of normal stresses has been evaluated (Figure 51) 1° show that it is 
always small in comparison to the production. For the present experiment, 
the production of normal stresses is about two orders of magnitude smaller 
than the production. 
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3.2 Separating FIov Along Convex Wall 


The setting up of the adverse pressure gradient for this experiment 
has already been examined in detail in Section 2 . 3 . 2 , therefore, this 
will not be discussed again. However, it should be pointed out that a 
separating flow was obtained on the convex test wall. Although it was 
not possible to locate exactly the separation line, indications are that 
separation occurred around Station 11, at x = 66.0 inches. In this 
connection, it should be mentioned that later boundary layer calculations 
predict separation to occur at x = 67 inches (see Section 4.6.1). 

Because of the disturbing influence of the secondary flows near the end 
walls of the curved section, most likely separation would not occur 
along a line parallel to the z axis. Yaw measurements at Stations 8 
and 10 (Figure 35), incidate that there is significant cross-flow at 
x = 62.5 inches. This means that there is convergence in the flow at 
this plane, and as a result, the flow would most likely separate first 
near the end walls than along the tunnel center line. Velocity profile 
measurements at Stations 8 and 10 (Figure 52) lend evidence to this as 
do the measurements of wall static pressure at planes 14 inches above 
and below tunnel center line. Because of this, the flow in this region 
can no longer be said to be two-dimensional. This fact is borne out by 
the momentum balance calculation. 

In order to make use of equation ( 3 . 1.1) for the momentum balance 
calculation, the skin friction C„ and the boundary layer integral 
parameters & and 0 have to be known. Since the Law of the Wall is 
independent of the free stream conditions } and has previously been 
verified for flow over convex surfaces , it is assumed to be valid for 
the present experimental flow. This means that can again be deter- 

mined from the Clauser plot. The velocity profiles are numerically 
integrated to give 5* and 0 . Because of the errors involved in the 

velocity profile measurements near separation (i.e. at Stations 9 and 11), 
-x- 

the calculated 5 and Q at these two stations would be smaller than 



the true values. This should he borne in mind when examining the 
momentum balance calculation. 

The results of the momentum balance calculation using equation (3.1-1) 
is given in Figure 40. It can be seen that up to x = 59 inches, the flow 
is fairly two-dimensional, but the flow begins to deviate from two-dimen- 
sionality after that. This is consistent with all other measurements. 

The development of 6 * and 0 is shown in Figure 83 and the 

development of and H is given in Figure 84. It can be seen that 

* * 

6 and 0 increase steeply towards separation. However, 5 increases 

much faster than Q , resulting in an extremely large slope for H near 

separation. Again, this lends evidence to the fact that the flow 

separates around x = 66 inches. Because of the errors in the velocity 

profile measurements near separation, the determined from the 

Clauser plot would also be in error. Since near separation, the measured 

velocity is higher than the true value, the determined from Clauser 

plot would tend to be greater. This is borne out by the fact that the 

measured at Stations 9 and 11 are greater than the calculated 

(see Section 4.6.1 and Figure 84). 

The measured velocity profiles are plotted in Figure $2, and the 
results show that the measurements near the wall at Stations 9; U and 12 
are in error, because the measurements give a near constant velocity in 
this region. In actuality, this is not the case, especially at Station 12, 
where the flow is known to have separated. Also shown are the off center 
line velocity measurements at Stations 3 , 5 , 8 and 10. Again, the 
results substantiate the claim of two -dimensionality at Station 4, but 
not at Station 9- 

The semi-log plot of the measured velocity profiles at Stations 4, 

7, 9 and 11 is given in Figure 53* The friction velocity u is calculated 
from the measured Up W and . Again, the velocity profiles show the 

existence of a wall law region, and that all the data points in this 
region fall on the Law of the Wall, thus giving support to the claim that 
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the Lav of the Wall also applies to flow over convex surfaces in an 
arbitrary pressure gradient. This should come as no surprise, because 
for plane flow, the Law of the Wall has been shown to be independent of 
the free stream conditions. 

By comparison with turbulent boundary layers over plane surfaces, 
it can be expected that the Defect Law also holds for boundary layers 
over convex surfaces in arbitrary pressure distribution. A semi-log 
plot of (Up-Uj/oj. verses y/A is given in Figure 54. Again, the result 
shows that A is also a function of s/r . Therefore, the skin friction 
relation as given in equation (3-1-5) I s also applicable to convexly 
curved turbulent flows in adverse pressure gradient. 

In previous experiment, the straight probe with a straight wire 
P — p 

was used to measure u , uv and w . The slanted probe with a straight 
wire was used to obtain uv and v^ , while w was assumed zero. 

However, for the present experiment, due to the disturbing side influence 
of the secondary flows at x = 62.5 inches, it was felt that the assumption 
of vw = 0 may not be a good one. Therefore, a straight probe with a 
slanted wire was used to measure uv , v^ and vw (see Sections 2 . 2.5 
and 2.4.3); while the same technique was used for the measurements of 
u , uw and w . The results of the measurements at Stations 1, 4, 9> 

11 and 12 are given in Figures 55 to 62. Since the flow separates around 
Station 11 , not too much significance should be attached to the measurements 
near the wall at Station 12 . This is because there is probably a small 
reverse flow region at Station 12 . However, the measurements at the outer 
part of the boundary layer would provide some insight ijito the nature of 
separated flows, therefore the data at Station 12 should be examined in 
this light . 

The rreasurements at Station 1 serve to establish the reliability 
of the slanted wire technique to measure uv . At the same time , the 
measured vw would also indicate the validity of the assumption vw = 0 
made in the measurements of the previous experiment. The measured uw 
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and vw are shown in Figures 55 and 56 respectively. In Figure 57 
the distribution of u^ is shown. It can be seen that at Station 1, 
nowhere are uw and vw greater than ^ of u . Also, the measured 
turbulence intensities and shear stress at Station 1 correlate well with 
the measurements at Station 1 of the previous experiment and with the 
data of Klebanoff (1955)* This establishes the reliability of the 
slanted wire technique to measure uv . It also provides evidence to 
support the claim that the flow is truly two-dimensional in the straight 
section. 

As to the measured uw and "vw in the curved test section, the 
results show that toward separation, uw and vw are as high as 2(yf 0 of 
u^ . This is particularly true near the edge of the boundary layer 

— n 

where u is relatively small as a result of the effect of convex 
curvature. At Station 4, uw and vw are about 10$ of u for the 
greater part of the boundary layer. However, since the error of the 
rotating-wire method as determined in the fully developed pipe flow 
experiment (Appendix C) and the measurements at Station 1 is about 5$ , 
the flow at Station 4 can still be considered to be fairly two-dimensional. 
At Stations 9, 11 and 12, even discounting the error of the rotating-wire 
method, uw and vw still amount to about 1 5$ of u . To the best 
of the author's knowledge, no measurements of uw and vw are available 
for boundary layers in an adverse pressure gradient. Therefore, it is 
hard to estimate whether the uw and vw are the result of the effect 
of convex curvature or the three-dimensional nature of the flow. 

However, from the results of the measurements of the constant pressure 
flow experiment, it is very likely that uw and vw arise as a result 
of the convergence of the flow towards separation. 

From the flat plate experiment of Klebanoff (1955)> it i s known that 
the fluctuation velocities (u 2 , v 2 and w 2 ) reach their maximum very 
near to the wall. The separating flow experiment of Schubauer and 
Klebanoff (1951) shoved that under the influence of an adverse pressure 
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gradient , these maxima slowly move away from the wall. The same 
phenomenon is observed in the present experiment, but the effect of 
convex curvature tends to reduce the maximum reached by the fluctuation 
velocities (see Figures 57-59)- Again, this is due to the fact that convex 
curvature acts to reduce mixing, and hence suppresses turbulence activity. 
In spite of the large adverse pressure gradient (which promotes turbulent 
mixing) the turbulence intensities in the outer part of the boundary layer 
(y/& > 0.4) are still significantly smaller than the corresponding 

values at Station 1 (Figures 57, 58 and 59). Compared with the results 
of the constant pressure flow experiment, it can be said that the effect 


of convex curvature is very significant in the outer part of the boundary 
layer, and the data seem to indicate that the curvature effect is indepen- 
dent of the pressure gradient. This can be seen from a comparison of 

Figures 48 and 60. For the constant pressure flow, q^Aj^ ~ .002 at 

pw — 

y/s = .4 and this drops slowly to zero towards the edge of the "boundary 
layer. For the separating flow q/Upw ~ *002 at y/5 = .5 , and this too 
drops to zero in the same manner as in the case of the constant pressure 

O 

flow. The fact that q. is approximately the same for Stations 4, 9 and 
11 in the outer part of the boundary layer is an indication that the 


effect of curvature is far greater than the adverse pressure gradient 
effect . 


The shear stress distributions for Stations 1, 4, 9, 11 and 12 are 
plotted in Figure 6l. Under the action of an adverse pressure gradient, 
the point of maximum shear stress moves away from the wall. But, how far 
should the maximum point be away from the wall, and how great is the 
maximum in the absence of wall curvature cannot be determined, because 
of the lack of data for a corresponding plane flow. However, anticipating 
the results of the boundary layer calculations outlined in Chapter IV, a 
rough estimate can be obtained. The result of such a calculation is shown 
in Figure 86. It can be seen that convex curvature acts to reduce the 


maximum reached by the shear stress and prevent the point of maximum from 
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moving too far away from the wall. This will he discussed further in 
Section 4.6.1. In spite of the adverse pressure gradient, the shear 
stress again vanishes inside the boundary layer at y/& ~ *5 • This is 
another indication that convex curvature in the outer part of the boundary 
layer has a much greater effect on the flow than the pressure gradient. 

To provide a chock for the measured shear stress profiles, the 
boundary layer equations (4.2-15) to (4.2-17) are numerically integrated 
using the measured velocity profiles and the C^, determined from the 
Clauser plot. The resultant shear stress profiles at Stations 4, 9 and 
11 are then compared with tire measured profiles. It was found that the 
resultant profiles were quite similar to those calculated by the present 
prediction method, and for the sake of clarity they were not shown. 

From Figure 86, it can be seen that the measured shear stress is much 
lower than the shear stress calculated from the measured velocity profiles. 
The discrepancy could be explained by the fact that the flow is not quite 
two-dimensional, and by errors in measurements (it is pointed out by 
Bissonnette (1970) that in regions of very high velocity gradient, the 
rotating -wire method for the measurement of -uv becomes inaccurate). 

Much the same discussion about turbulent energy production in the 
previous section also applies to the results of this experiment, and 
hence will not be repeated. However, it should be pointed out that although 
the flow is under the influence of a severe adverse pressure gradient, 
there is very little turbulent energy production after y/s ~ • Even 

in the region near the wall, energy production is severely curtailed by 
the effect of convex curvature. Since the shear stress measurements are 
in error, therefore, the energy production calculations are only qualita- 
tively correct. 

3.3 Constant Pressure Flow along Concave Wall 

The z direction traverse at Station 1 serves two purposes. First, 
it serves to indicate the kind of disturbances created by the presence of 



the travel’ cing device in the tunnel, and second, it serves to provide 
a basis for comparison with subsequent traverses in the concave test 
section. Since the flow in the straight section is shown to be two- 
dimensional (see discussion in Section 2.3.1), the velocity inside the 
boundary layer and at a fixed distance away from the wall should be 
constant for all z planes, flier ef ore, if the traversing device 
(Figure 21) is set up properly inside the tunnel, the hot-wire would 
give a straight line on the X-Y plotter for the z direction traverse 
at Station 1. Indeed, this was the case, and the result of five different 
traverses at y = 0.2, 0.4, 0.6, 0.8 and 1.0 inches from the wall is 
shown in Figure 63. This substantiates the previous claim that the flow 
is two-dimensional in the straight section, and that the disturbances 
created by the presence of the traversing device do not disturbe the 
nature of the flow. However, it does increase the local flow velocity 
because of the blockage effect it has on the flow. To estimate the 
amount of increase, the traversing device was set at two different positions 
at Station 1 such that the hot-wire was 0.2 and 1 inches away from the 
wall. The hot-wire was then moved down to the tunnel center line plane 
and the mean velocities at these two y positions were measured. Compar- 
isons with the velocity profile measured in the same location show that 
the velocity increase varies from 1# to less than % of the free stream 
velocity, and that the increase is greater in the free stream than near 
the wall. The velocity increase is due to the relatively large cross- 
section of the strut (aerofoil shape with maximum thickness of l/2 inches 
and a chord of 2 l/2 inches) in a tunnel of 6 inches wide. Due to the 
peculiar design of this traversing device, it is not possible to make a 
quantitative study of the increase in velocity across the boundary layer. 

As a result, all traverses obtained in the z direction can only be 
interpreted qualitatively. 

The X-Y plotter trace of the z direction traverses at x = 70 
and 96 inches are shown in Figures 64 and 65 respectively. Again, traverses 
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verc taken at y = 0.2, 0.4, 0.6, 0.8 and 1.0 inches respectively. 

All five traces vere plotted on the same graph. Compared to the traces 
obtained at Station 1 (Figure 63), the results given in Figures 6 k and 
65 clearly indicate that there is significant velocity variations in 
the z direction. These velocity variations can be explained by 
assuming the existence of a system of longitudinal vortices similar to 
the Taylor -Gor tier type vortices inside the boundary layer. Then the 
positions of the high points (crests) la the trace could be taken to 
correspond to the positions between two vortices whose flow directions 
are directed towards the wall, and the positions of the low points (trough) 
could be taken to correspond to the positions between two vortices whose 
flow directions are directed away from the wall (see Figure 66). In the 
positions of the crests, faster moving fluid elements are being entrained 
into the boundary layer by the vortices, and as a result of the turbulent 
mixing process inside the boundary layer, the mean velocity becomes higher. 
On the other hand, in the positions of the troughs, the vortices sweep up 
slower moving fluid elements from the wall, and through mixing, the mean 
velocity becomes smaller. Thus, the "wave like" shape of the trace in the 

z direction. 

The data also show that there are more than one pair of vortices 
inside the boundary layer. For the 21 inches traversed in the central 
core of the tunnel, there are about 10 pairs of vortices at x = TO kiches 
and 7 pairs at x = 96 inches. This indicates that the vortices spread 
in the z direction and some vortex pairs end at the end walls of the 
tunnel as the flow moves downstream. The mean half wave length at x = 70 
inches is approximately 0.57 inches, and it is about 0.9 inches at x - 96 
inches. In addition, it can be inferred that the positions of the crest 


and trough are quite stationary, and that they remain at the same z 
position across the boundary layer. These results substantiate Tani’s (1962) 
findings that a system of Taylor -Gortler type vortices does occur an 
turbulent boundary layers over concave surfaces, and that the wave 
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amplitude does decrease as the radius of curvature increases. 

It should be pointed out that the above data only man age to show 
the presence of the vortices up to y = 1.0 inch. No traverses were mde 
above y = 1.0 inch at both x = 70 and 96 inches. This is most unfortunate 
because the maximum y position is limited by the width of the curved 
test section and the relatively large size of the traversing device. The 
traversing device is located outside the boundary layer for all the traverses 
at x = 70 inches, and partly inside the boundary layer at x = 96 inches. 
However, since the hot-wire is 2.75 inches away from the strut, the 
disturbances created by the strut will not be felt at the hot-wire, except 
that the flow between the strut and the wall will be slightly accelerated. 
This is certainly true at Station 1, and it is assumed to be true also 
at x = 70 and 96 inches respectively. 

Since the mean output of the hot-wire only gives tte normal com- 
ponent of the resultant velocity, and the flow direction inside the 
boundary layer is not known, it is not possible to construct a quantitative 
diagram of the vortex structure inside the boundary layer from the traces 
shown in Figure 64 and 65. As a result, only a qualitative diagram is 
given in Figure 66. The appearance of two layers of vortices at x = 96 
inches will be explained later. From this s ample-minded diagram, the z 
locations of the points of measurements can be determined. For Station 1, 
the measurements were taken at tunnel center line, i.e. z = 0 , while 
for x = 70 and 96 inches, measurements were taken at the position of 
the crest and trough of the wave system. The reason for this is obvious. 
Assuming the presence of the vortex system, there is a minimum of cross 
flow at the positions of the crest and trough. Those positions that are 
nearest to the tunnel center line and lie above it are chosen to be 
the points of measurements, because at these locations, the variation 
of the positions of the crest and trough is a minimum across the boundary 
layer. The locations of the trough and crest are labelled Stations 2 
and 3, 4 and 5 n't x = 70 and 96 inches respectively. The distance 


59 



between Stations 2 and 3 is O .672 inches and between Stations 4 and 5 
is 1.0g2 inches. The x and z coordinates of the measuring stations 
are given in Figure 6 . 

The velocity profiles measured at Stations 2, 3 , 4 and 5 are plotted 
in Figure 67 . Since the cross-flow is a minimum at these locations the 
measured velocities are not corrected for cross -flow. Instead V is 
assumed to be zero at these locations. The measured velocities are not 
corrected for V , which can be significantly greater than the corresponding 
flat plate flow. However, the results show that s/r varies between .07 
and .18, and since v/u can be at most of the order of ( 5 /u )(c)U/c3x) , 
therefore, the maximum resultant velocity Q ~ 1.0l8u . In this sense, 
the measured profiles can be considered as profiles of U . 

From the measured profiles, it can be seen that the boundary layer 
is thicker at the position of the trough than at the position of the 
crest, and the growth rate at these two positions is approximately the 
same. The thicker boundary layer at the trough can be explained as 
follows. Concave curvature enhances mixing, so when a fluid particle 
is displaced from a position nearer the wall to a position further away 
from the wall, it will tend to move even further out under the action 
of the centrifugal force. At the position of the trough, the particle 
is pushed even further out by the vortex motion, and as a result, the 
boundary layer becomes thicker. At the position of the crest, the 
vortex motion acts in the opposite direction to that of the centrifugal 
force, thus the boundary layer becomes thinner compared to that at the 
trough. The results also show that the boundary layer thickness at the 
position of the trough is about twice that at the position of the crest. 

To find out whether the Law of the Wall also holds for flow over 
concave surfaces, the velocity profiles are plotted in the Clauser form. 

The values of were obtained by extrapolating the measured profiles 

to y = 0 . This is because no wall static pressure measurements were 
made at the positions of the crest and trough. A sample plot of the 



velocity profiles at Stations 2 and 3 in the Clauser form is given in 
Figure 68. The plot shows the existence of a log region near the vail, 
but no Clauser lines match the data. The Clauser lines match the data 
points in the viscous region. However, this should not be taken as 
given the C^ at these locations. Therefore, the conclusion can be drawn 
that the Law of the Wall does not hold for flow along concave surfaces. 

As a result, the skin friction at the vail cannot be determined. 

Since convex curvature suppresses mixing, therefore, concave curva- 
ture should enhance mixing. The result of this would be a significant 
increase in the fluctuation velocities inside the boundary layer. 

Evidence to this fact is given by the turbulence measurements at Stations 
1 to 5* The results are shown in Figures 69 to 72. In analysing these 
results, one point should be noted, and this is the effect of the adverse 
pressure gradient on the flow at the entrance to the curved section. 
Bradshaw and Ferriss (1965) showed that in a relaxing boundary layer, 
the turbulence intensities decrease as the pressure gradient is removed. 
Also, in passing from a zero pressure gradient region into an adverse 
pressure gradient region, the turbulence intensities increase significantly. 
For the present case, if the concave curvature has no effect on the flow, 
then the turbulence intensities should increase from x = 24 inches to x = 

70 inches and decrease significantly from x = 70 inches to x = 96 

“0 / 2 

inches. The distribution of u /U^ does show a decrease in the region 
near the wall from x = 70 inches to 96 inches, but remains relatively 
constant in the outer part of the boundary layer (Figures 69). This is 
consistent with the findings of the convex wall, i.e. wall curvature has 
very little effect in the region near the wall. However, the same trend 
is not noted in the distribution of v^/u^ and w^/u^ (Figures 70 and 
71). There is no discernable decrease in v^ and w^ distribution 
across the boundary layer from x = 70 to 96 inches. 

According to the component of the turbulence energy equation 

(E-2), concave curvature promotes v^ production, and this is evident 
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from the distribution of \t/Up V (Figure 70). However, not all of 
the increase over that measured at Station 1 is due to the influence 
of concave curvature. Part of the increase is due to the adverse 
pressure gradient effect at the entrance to the curved section ; and 
part of the increase is due to the effect of the vortex system. At 
Stations 2 and the vortices have a resultant motion in the yz 
plane directed away from the wall. This contributes further to the 
outward radial movement of the fluid particles. At Stations 3 and 
the vortices have a resultant motion in the yz plane directed towards 
the wall, thus any radial movement of fluid particles outward would 
be hindered by the vortices. As a result, the v 2 at Stations 2 and 4 
are greater than at Stations 3 and 5 respectively. This is particularly 
true at x = 96 inches. The reason is that there may be more than one 
layer of vortices at x = 96 inches. This point will be discussed again 
later. The same is true of the v 2 /u 2 v distribution (Figure 71) < In 
other words, the system of longitudinal vortices contribute to the 
production of v 2 and w 2 fluctuations, but not u 2 . This is evident 
from the fact that near the wall, u 2 drops from x = 70 to 96 inches, 
but not v 2 and w 2 . As a result of this, the distribution of turbu- 
lence energy q. 2 /lJp V does not change significantly from x = 70 to 96 
inches (Figure 72). Note that the turbulence energy distribution is 
greater at Stations 2 and 4 than at Stations 3 and 5* This is especially 
true at x = 96 inches. 

The shear stress profiles at x = 70 inches are given in Figure 73; 
while those at x = 96 inches are shown in Figure 74. The shear stress 
profile at Station 1 is also shown for purpose of comparison. The fact 
that the shear stress profiles at x = 70 inches show a maximum is due 
to the adverse pressure gradient effect at the entrance to the curved 
section. The maximum remains even at x = 96 inches. This indicates 
the influence concave curvature has on the mixing piocess. However, 
the shear stress profiles at x - 96 inches show two points of maximum. 
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This could be due to the fact that there is more than one layer of 
longitudinal vortices at x = 96 inches. Turbulence intensity measure- 
ments also seem to lend evidence to such a postulate. However, the 
strongest evidence comes from the profiles of uw and vw at x = TO 
and 96 inches (Figures 75 and 76). Measurements of uw and vw at 
Stations 2 and 3 show that there is only one maximum (Figure 76), while 
the measurements at Stations 4 and 5 distinctly show two maximum points 
for both uw and vw . It can be postulated that at x = 70 inches, 
there is only one layer of longitudinal vortices, while at x = 96 inches, 
there are two layers. This is depicted in Figure 66. Also shown are 
the velocity traverse at a constant distance away from the wall and 
the boundary layer thickness distribution in the z direction. Of 
course, within each bigger vortex, there are smaller vortices, but these 
are not shown. 

From previous discussion, the directions of rotation of the vortices 
at x = 70 inches can be postulated to be as shown in Figure 66. Based 
on the same evidence, the directions of rotation of the first layer of 
vortices (the layer nearer to the wall) at x = 96 inches are assumed to 
be the same as that at x = 70 inches. The directions of rotation of 
the second layer of vortices can be postulated after an examination of 
Figures 75 and 76. At Station 2, the distribution of uw is positive 
across the boundary layer, while that of vw is negative (Figure 75). 

At Station 4, the distribution of uw is also positive across the 
whole boundary layer, while that of vw is also negative (Figure 76). 

This seems to indicate that the second layer of longitudinal vortices 
would also have the same direction of rotation as the first layer. 

At Station 3 , the uw distribution is negative , and the vw 
distribution is positive. The same is true at Station 5 . This 
again lends evidence to support the postulate that the directions of 
rotation of the second layer of vortices at x = 96 inches are the same 
as the first layer. The fact that there are two layers of vortices at 
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x - 90 ir-chec in very clear from Figure 76. However, the vortex structure 
A ~ inches in not stable. Therefore, further downstream, it may 
revert back to that of a single layer structure. 

The distributions of uw and vw at Stations 4 and 5 all show 
the existence of two waves with the waves going to zero at the wall, at 
y/b ~ .6 and at the edge of the boundary layer. Based on all these 
evidences, the resultant vortex structure as shown in Figure 66 is 
postulated for the boundary layer flow at x = 70 and 96 inches 
respectively 1 '. 

Finally’ the shear correlation coefficient and the turbulent energy 
production are calculated and plotted in Figures 77 and 78. As expected, 
the shear correlation coefficient does not remain constant across the 
boundary’ layer. This is due to the effect of the longitudinal concave 
curvature which gives rise to the Taylor -Go rtler type vortices. As to 
the peculiar behaviour of the turbulent energy production at x = 96 
inches, this is probably due to the peculiar velocity profiles at Stations 
4 and 5- 



Chapter TV 

THEORETICAL PREDICTIONS OF 
CURVED TURBULENT BOUNDARY LAYERS 


4.1 Introduction 

The ultimate goal in boundary layer research is the understanding 
of the physics of such flow and at the same time to be able to develope 
a procedure for pi'edicting, both quantitatively and easily, the piactical 
features (e.g. H, C f , and the shear stress) of boundary layers. In 
the past, investigators have been quite successful in the development 
of prediction methods for two-dimensional boundary layers. In these 
calculation methods, be they based on integral or differential equations, 
it has been assumed that the static pressure variation across the layer, 
due to either rapid boundary layer growth, flow injection at the surface, 
or surface curvature (when j o /r j « l) has very little effect on the 
flow and hence can be neglected. Such an assumption neglects the effects 
of the curvature of the mean flow streamlines. These effects are shown 
to have a considerable influence on the mixing process in the flow. 

Convex curvature in the mean flow streamlines reduces the mixing between 
the fluid layers and this leads to an early cut off of the shear stress 
(see Fig. 59). In contrast, concave curvature enhances the mixing pro- 
cess, and the shear stress profile remains quite full until the edge of 
the boundary layer where it drops to zero steeply (see Figs. 73 and 7^-). 

Thompson (1963) was the first author to point out the importance 
of incorporating a curvature parameter in calculation methods. In a 
critical review of existing two-dimensional calculation methods, he found 
that the small curvature that existed on most of the test surface where 
measurements were nade , was not large enough to cause substantial 
difference between the measured velocity profiles and that calculated 
assuming the flow to be two-dimensional plane flow. However, there was 

between the measured and calculated shape factor 


a consistent difference 



case vas the 


development in the streainwisc direction. A well-known 
measurements on an aerofoil by Schubauer and Klebanoff (1951)- In 
their measurements , the curvature parameter o/r , as defined by 
Thompson, vas ~ 1/150 downstream of tic pressure minimum. Most pre- 
dictions methods gave fairly good correlation for the velocity pro! lies 
and the momentum thickness, but all predicted the chape factor too low. 
Thompson argued that the effect of curvature vas on the entrainment 
process. He modified Head's entrainment function by the inclusion of 
an empirical factor vhich vas assumed to be a simple function of 
the curvature parameter &/r . This improved the agreement betveen the 
experimentally measured H and that predicted by his calculation 
method. Although the method used by Thompson was not very satisfactory, 
his results brought to light the very important conclusion, namely that 
the primary influence of the streamline curvature vas on the mechanism 
of the turbulent motion and the entrainment of free stream fluid into 
the layer. Therefore, ar.y attempt to extend existing two-dimensional 
calculation methods to include curvature effects should be directed to 
the modification of some parameter, e.g. Head's entrainment function 
or eddy viscosity, that will take the physics of the flow into account. 

More recently, Bradshaw ( 1968 ) has drawn the analogy between the 
Richardson Number, vhich is a meteorological parameter, and a curved 
flow parameter "L" vhich describes the effect of streamline curvature on 
turbulent flow. Using this analogy to apply meteorological data to curved 
turbulent flows, Bradshaw showed that the apparent mixing length was 
affected appreciably even though 5/fc ~ l/300 . Incorporating this 
modified mixing length into the Bradshaw, Ferris and Atwell calculation 
method, Bradshaw managed to obtain better agreement between H , 5 and 
C as measured by Schubauer and Klebanoff (1951) and that predicted by 
his calculations. Again, this points to the need to modify the eddy 
viscosity or mixing length function, such that the resultant form will 
be general enough to predict both curve and plane turbulent flows 
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accurately. 

Hence , the object of the present chapter is to seek a set of self 
consistent turbulent boundary layer equations for flows along curved 
surfaces, and the corresponding momentum integral. This is done in 
Section 2 and Appendix D. 

In Section 3 & modified eddy viscosity function which consists of 
the product of the flat plate eddy viscosity and a factor that is a 
simple function of the curvature parameter is derived. This modified 
eddy viscosity approaches the eddy viscosity function put forward by 
Mellor and Gibson (1966) as the surface curvature becomes zero. A 
technique similar to that proposed by Mellor and Herring (1970) is 
adopted for the derivation of the modified eddy viscosity function. 

Assumption s a re made to simplify the turbulence energy equations of 

2 2 2 * 

u , v , w and uv for curved turbulent flows (for derivation of 

these equations see Appendix E) so that all terms in these equations, 

including the pres sure -velocity correlations and the triple velocity 

correlations terms, can be expressed in terms of the double velocity 

correlations. These equations are further simplified by assuming the 

advection and diffusion terms to be small compared to the dissipation 

and production terms, and hence can be neglected. This is equivalent 

to assuming that energy production balances dissipation and that a state 

of equilibrium is reached as far as the energy distribution is concerned. 

This is a plausible assumption because existing flat plate (Klebanoff 1955) 

and pipe flow (Laufer 195^) data do indicate such a state of equilibrium 

the flow near the wall. The resultant equations are algebraic and 

can easily be solved for the shear stress. 

In relating the various terms in the turbulence energy equations 
to u^u f length scales are introduced as proportionality constants. 

Therefore, in order that the shear stress be specified completely by 
mean flow quantities only, empirical statements for these length scales 
should be derived. Mellor and Herring (1970) have shorn that the four 



length scales are not truly independent and from experimental correlation 
two relatione can be found for the four length scales. This leaves 
two length scales to be specified. If the eddy viscosity function for 
curved turbulent flow as derived in Section 3 is truly a general one, 
then the length scales appearing in the expression should be independent 
of curvature and pressure gradient. Hence, it is sufficient to obtain 
empirical statements for the length scales by only considering the plane 
flow turbulence energy equation and the Lav/ of the Wall. This is done 
in Section 4. 

The boundary layer equations together with the eddy viscosity 
function are then solved using a finite difference method adopted by 
Herring and Mellor (1970). Reduction of the set of equations to an 
ordinary differential equation is given in Section 5* The calculations 
are then compared with the measurements on both the convex and concave 
surfaces. Finally, a brief discussion of the calculation method and 
the comparisons are given in Section 6. 

4 . 2 Boundary Layer Equations for Curved Flow and the Momentum Integra l 

The Navier-Stokes equations for a constant density incompre ssible 
flow can be written as : 

= 0 (4.2-1) 

+ q*Vq = - ~ VP + VV q. (4.2-2) 

Consider a flow over a two dimensional curved surface. Use general 
orthogonal coordinates with x measured along the surface, y normal 
to the surface and z at right angles to the x-y plane which is the 
plane of the motion. If k(x) = -Vr(x) i s "the curvature of the surface, 
(k is taken to be positive for convex curvature, and negative for con- 
cave curvature) the elements of length along the parallel curves and along 
the normal are h^ = 1 + ky and h^ = 1 . The element of length along 
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the z direction is = 1 . Taking u , v , and v to be the 
velocity components along the x , y and z directions respectively, 
and working out the components of the vector equation of motion (see 
Owczarck 1964), the following equations for u , v and w are obtained: 
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arc the shearing streamer.. It should be pointed out that the pressure 
p and the stresses T. . have the dimension of velocity square in these 
equations . The continuity equation can ho vrit-ben 
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Consider a turbulent flow that is two dimensional in the mean. 
Therefore, the mean fj.ow in the z direction is zero, and all ^ of 
mean flow quantities vanish. The field equations of motion for such a 
flow when decomposed into ensemble mean velocities U , V plus fluctuatxng 
velocities u , v , w and their concomitant pressure P and p are: 
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dw 1 d , ^ 
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^ij arC fluctuating viscous stresses given by the same 

expressions as T except that all velocities in T are now replaced 

ij * 

by the fluctuating velocities. 

The only component of vorticity that is non-zero is that normal to 
the x-y plane and is given by: 

r _ _1 dV _ dU _ k 

^ 1+ky dx dy 1+ky U (4.2-14) 
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With the field equations of motions given by equations (4.2-7) to 
(4.2 -9) and the vorticity defined by equation (4.2-14), the next step 
is to simplify these equations to obtain a set of self-consistent boundary 
layer equations for a turbulent flow that is two -dimensional in the mean 
along a curved surface. The assumptions to be made are that the boundary 
layer thickness 5 is very small compared to the x dimension and the 
radius of curvature R of the surface is of the same order as the x 
dimension. Therefore, &/r is very much smaller than 1 . The geometry 
of the surface is specified by R(x). The function R(x) has to be 
smooth but otherwise it can be a general one. 

Recently, Mcllor (1970) has demonstrated that the method of matched 
asymptotic expansion can also be applied to turbulent boundary layers. 

Two small parameters, e - u. /[J and e e = v/u £ , (where u is any 
characteristic turbulent velocity, U Q is any characteristic free stream 
velocity and £ is a characteristic length) appear as a result of 
making the equations non-dimensional with respect to u^_ , and £ . 

However, Mellor (1970) has shown that expansion in one parameter, e , is 
sufficient and that (e/c 11 )^ 0 as £ 0 for arbitrary n . Unlike 

the case of laminar boundary layers, three layers exist, (an outer 
layer, a middle layer and an inner layer) due to the presence of the 
two small parameters. Also, the law of the Wall is shown to be a con- 
sequence of the matching between the inner and the middle layers. In 
view of this, the set of boundary layer equations for a curved flow is 
obtained through the method of matched asymptotic expansion instead of 
through dimensional argument. This is carried out in Appendix D for 
equations (4.2-7) to (4.2-9) and equation (4.2-14). The resultant 
set of equations in dimensional form Is given by equations (d- 22) to 
(D- 25 ) and equation (D- 26 ), which is 
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(4.2-16) 
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(4.2-17) 

(4.2-18) 

(4.2-19) 


where is the potential velocity which is, as yet, undefined. 

The inner boundary condition is given by the no slip condition. 

The outer boundary condition can be defined so that the velocity 
approaches the free stream velocity as y -* 5 - The free stream velocity 
is obtained from the Bernoulli equation and the condition of zero 
vorticity. This immediately gives the free stream velocity distribution 
as : 


u (x,y) = u (x)e 

p' pw 


-ky 


(4.2-20) 


where U^(x,y) is the potential velocity and U^^x) is the potential 
velocity at the surface. 

It should be noted that, inserting (4.2-20) into (4.2-19) shows 
that x/p does not vanish as y 

3 + 2k i 


if v is maintained constant, 
e 


On the other hand, the terms 


0 as y 


This 


6y v P' " P 

circumstance is identical to that obtained in laminar flow. This is 

pointed out because, later on, an eddy viscosity model where does 

stay constant for large y is assumed. This is, o'f course, not 

particularly realistic but it is known that predicted velocity profiles 

are insensitive to the detailed behavior of at large y . 

Equation (4.2-17) must be used together with the other approximations 

of the boundary layer equations (see Appendix D). If — — - = k is 

P oy 

used instead, the pressure term will not balance the advective terras 
in equation (4.2-16) when it is integrated over the whole layer. Thus, 
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integrating equation (4 *2 -IT) snd then differentiate the resultant 
expression with respect to x it can be shown that 


1^P_ i a fl 2 -2ky. 
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into equation (4.2-16) and then integrate from 0 to 
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*The :Lnt egrat ion is somewhat complicated. It is helpful to first rewrite 
the shear stress terms in equation (4.2-16) as p _1 e" 2 ^ S(e v x)/dy and 
then set U = c" hy u(x,y) where k = k(x) and u ~ U pw as y - » . 
Solving the continuity equation gives 


V = 
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-ky 


r r 6u _ ^ dk r 

" ^o Sx ^ dx o 


y u dy] 


Substituting the above U, V terms and equation (4.2-21) into equation 
( 4.2-16) the integration way be performed and the result cleared of 
exponential functions. Finally, noting that u = U pw (u/u p ), the 
result can be written in the form of equation (4.2-22). 


7 ? 


Note that the two integrals in equation (4.2-22) are finite since u/u -* 1 

1 ^ P ^ 

at the edge of the bounaary Jayer. If - — = kU^ had. been uoed instead 
of equation (4.2-17), a term vould be introduced which would blow up 
ao the limit of integration approached infinity. 

The set of equations (4.2-15), (4.2-16), (4.2-17) and (4.2-18) 
together with the eddy viscosity to be derived in Section 3 form 

a closed set, and can be solved when the appropriate initial and boundary 
conditions are specified. The boundary conditions are given by: 


U(x,0) = v(x,0) = 0 


(4.2-26) 


at the wall, and 

U(x,y) = U p (x,y) = U pv (x)e' ky (4.2-27) 

at the edge of the boundary layer. To solve this set of equations , the 
method of Herring and Mellor (1970) is used. The reduction of these 
equations to a single ordinary differential equation is given in Section 


^ -3 The Eddy Viscosity Hypothesis 

As pointed out by Mellor and Herring, an expression for the shear 
stress ( -uv) can be obtained from the turbulence energy equations when 
further assumptions arc* made to simplify the equations to a set of algebraic 
equations. Invoking the assumptions that in the region near a wall, the 
advection and diffusion of turbulence energy are very much smaller than 
the production and dissapation of turbulence energy, the advection and 
diffusion terms in equations (E-l4) to (E-17), as derived in Appendix E, 
can be neglected. The result- is 
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where A and £ are length scale as defined in equation (E- 6 ) and (E-7) 
respectively. 

The equations are correct to order 8 because of the inclusion of 
the curvature terms. Since the objective of the present Section is to 
find an eddy viscosity function that includes curvature, this purpose 
will not be served if the curvature terms in the turbulence energy 
equations are excluded. 

Omitting all the algebra, the result of solving equations (4.3-1) 
to (4.3-4) for ( -uv) is 
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It should be noted that in deriving this expression for (-uv), isotropy 

is not assumed. If isotropy is assumed, the same expression will still 

3 /2 

be obtained, except that the factor (l - 6 ^^/a) becomes 1 . In 

this case, when curvature is assumed zero, the expression. 
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is obtained. Immediately, it can be recognised that the eddy viscosity 

. 1/3, 3/2 Su . . 

can be identified with the expression (iUA ) r— • Similarly, one 

2/0 3/2 0/2 xy 

can identify (i^A ' ) (l - 6z/h) ' as the eddy viscosity in a 
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plane flow when isotropy 
by v ti and after some 
as 


is not assumed. Denoting the latter expression 
rearrangement , equation ( 4 . 3 - 9 ) can be written 
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The shear stress is given by equation (4.2-19) 
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where v is the eddy viscosity. Comparing equation (4.3-6) with 
e 

equation ( 4 . 2 - 19 ) can be written 
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This then is the desired eddy viscosity function. Note that for 
convex curvature k is positive; and the modifying factor within the 
curly brackets is always smaller than 1 , except at the wall. Hence; 
the physics of convexly curved flow, which in simple terms can be 
characterized by the decreased mixing; is embedded in the eddy viscosity 
function. For flow over concave surfaces, the present investigation has 
demonstrated that a system of longitudinal vortices exists inside the 
boundary layer. As a result; the flow is no longer two-dimensional, 
and the present approach to characterize the flow by a scalar eddy 
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viscosity function is in doubt. However, the eddy viscosity function 
as given in equation ( U . 3 -T ) affords the practical engineer a simple 
approach to estimate, though approximately, the various features of 
the turbulent boundary layer on a concave surface. This would involve 
some kind of assumptions about the flow or the prediction method (This 
will be discussed in detail in Section 4.2-6). For the present, it is 
enough just to point out the limitation of the eddy viscosity function. 
Now, it remains to define p which involves the ratio {zjk), and this 
is given in the next section. 


4.4 Determination of (f^/A) 

The eddy viscosity for a turbulent flow along any smooth curved 
surface is completely defined by equation (4.3-5)- Theoretically, v g 
can be computed once the variation of Z ± and A across the boundary 
layer is known. This requires emperical statements for Z ± and A very 
near the wall, which at present cannot be obtained with certainty due to 
the lack of reliable turbulence data in this region. Since the present 
objective is to find a modifying function that includes curvature for 
the flat plate eddy viscosity, it suffices to identify 
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with the eddy viscosity hypothesis put forward by Mellor and Gibson (1966). 
Therefore, only the ratio (^/a) retrains to be determined in the expression 

for v 

From the emperical statements of Z ± and A given by Mellor and 
Herring (1970), it can be seen that {z ± / A) is practically constant in the 
overlap and defect regions of the boundaiy layer. Due to the fact that 
the flow near the wall is predominantly viscous in nature, the curvature, 
according to Van Dyke (1962), is a second order effect in this region. 
Hence, it is sufficient only to determine the ratio Uj/a) in the outer 


77 



part of the boundary layer, and this can be obtained, most conveniently, 
by considering the Law of the Wall region only. 

Near the wall, Mellor and Herring (1970) assume that 

& ~ ( y > 1 > v ) 

or 

‘-f - % to ' X - “ 

S imilarly 


Aa 

v 


- V x) 


As yet, the functions ^(x) and ^(x) are undefined. For the Law 
of the Wall region, it is further assumed that these functions can be 
a PProximated by a power law namely 



A £ 
v 


m 

B (32 

V 


(4.4-2) 


where A and B are constants and n and m are integers to be 
determined. 

The eddy viscosity function should be applicable to both curved 
and plane flows with arbitrary pressure distribution. Therefore, it is 
necessary only to determine z ± and jv for the flow along a flat plate. 
The turbulence energy equations for ? , ^ ^ and «y , neglecting 

the diffusion and advection terms are; 


i a 

3 i. 




12 2 q 3 

3 0 + 3f 


+ 2 


uvf- = 0 
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(4.4-3) 
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In this region, the Lav of the Wall also hold and 


Su _ 

Sy Ky 




2 T w 


vhere ^ - — is the friction velocity and K is the von Karman 
constant. Making the velocities non-dimensional with respect to u 
and with the substitution of equations (4.4-1), (4.4-2) and (4.4-7) 
equations (4.4-3) to (4.4-6) can be written as 


u -n +2 

JL /£E\ q , +2 1+2 


— m 

A ' v ' 
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where 
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It is further assumed that the shear stress is constant, i.e. ( -uv ) - 1, 
in the Law of the Wall region, and the function 4*^ and <t>^ are universal, 
hence } v<*>^ and V< ^A are independent of viscosity. This necessary 

implies that n = m = 1 • With these simplifications, equations (4.4-8) 
to (4.4-11) reduce to: 


Bq_ + (u +2 - ~ q +2 ) + 2Aq +3 - 6 ^ 


3 


K 


t 1 +2 \ o A +2 p 

B(v - ~ q ) + 2Aq. = 0 


B(w +2 -iq +2 )+ 2Aq +2 = 0 


+ , _+s ^ 3Av 
q (-uv ) + — — 


+2 


= 0 


Solving these equations give: 


+2 +2 
v + 
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B 



4 

K 


= a 3 b (i - 6 ^) 3 


(4.4-12) 

(4.4-13) 

(4.4-14) 

(4.4-15) 

(4.4-16) 

(4.4-17) 

(4.4-18) 


A value of iJk which gives the best agreement in C f between the 
calculated values and the measured values for the case of constant pressure 
flow along convex wall is chosen and used for calculations of the other 
two cases. However, whatever value one chooses for l^jk , it should 
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satisfy equation (4.4-lB), so that the von Kamian constant can always 
be recovered. 




• j 


Reduction Of the Set 
Different ia.i Eg- j.a t ion 

The set of equations 


o f Boundary Layer Equations to an Ordinary 
(02-15), (4.2-16), (4.2-17) and (4.2-18) 


together with boundary conditions (4.2-26) and (4.2-27) can be reduced 
to an ordinary differential equation through the following substitutions. 
Following Herring and Mellor (1970), assume 


6f 


f’(x,n) 



(4.5-1) 


where 


JL 


Differentiate equation (4.5*1) with respect to x and y and integrate 
the continuity equation (4.2-15) to obtain V y it can be shown that 


u = e" 2aT1 [u_ul (1-f’) 2 - d-f’) If 
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(4.5-3) 


+ <»'" - v f s + ^ i 


where a = k5* and G = fdtj . The primes in a , 5 and refer 
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to ordinary differentiation with respect to x . 

The pressure term can be obtained by differentiating the integral 
of equation (4.2-17). The result is; 


1 cfP 
p £)X 


-2an 
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[U u' 

pw pw 


U pw , * ' 
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' *S 1 

- a 5 )t)] 


(4.5-4) 


Bearing in mind that 1 + ky is taken to be 1 only whenever 

it appears as a coefficient and noting that v is a function of n 

e 1 ’ 

the shear stress terms can be written as 




dy'P 


+ e &Tl a(<ff ) - e dT1 a(2<t>' + a<|>)(l - f * ) ) 




(4.5-5) 


where 
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U 6 
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Substitute equations (4-5-2) to (4.5-5) into equation (4.2-16) 
and after some rearrangement, the following equation for f is obtained, 
namely: 


bV . 
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(4.5-6) 


+ a(4>f ) - a(2<!> *+a<)> )(l-f ' ) 


This equation (4.5-6) can be further reduced to an ordinary differential 
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equation through the use of the identities 
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for the differential with respect to x . 

The resultant equation as function at x. can be written as 
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The overbar denotes average value and the CVs are all calculated from 
the previous iteration. Rearranging the equation further, the final 
equation to be solved can be written as: 


< V 4 = \ ♦ Vi + Vi + Vi 


(It. 5-9) 


and 
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(It. 5-10) 


b 2 C 2 “ C 3 


b l = - C 4 


The boundary conditions at the wall are: 


f (x,0) = 1 
f(x,0) = 0 


(4.5-11) 


and the free stream boundary condition is 


Lim '/ s 

f (x,t)) = 0 


(4.5-12) 


The unknown function b r is given by equation (4.3-6). 

5 

it can be vritten as 
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(4.5-13) 
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is the flat plate eddy viscosity function. 


where $ 

F 


U. 


pw 


F_ 

6 * 


The equations (4. 5 “9) snd (4.5-13) together with boundary conditions 
(4.5-11) und (4.5-12) are solved using as initial conditions the measured 
data at Station 1 of the three experimental cases , namely: (l) constant 

pressure flow along convex wall, (2) separating flow along convex wall 
and (3) constant pressure flow along concave wall. The results are then 
compared with the present data. A discussion of this comparison is 
given in the next section. 


4.6 Comparison with Present Data 

The boundary layer program of Herring and Mellor (1968) is used to 
numerically integrate the ordinary differential equation (4.5-9) with 
the set of boundary conditions (4.5-11) and (4.5-12). The program uses 
a fourth order Runge-Kutta technique for the numerical integration. 

Although such a technique works well for all types of boundary layer 
development on a flat surface, it fails to give a solution that converges 
to the required accuracy on both f and f* when the surface is curved. 
The difficulties occur at the point where curvature begins, and in part 
is due to the fact that the shear stress vanishes at about half the 
boundary layer thickness at this point. In order to overcome this short- 
coming, another integration technique is used. The method is discussed 
in Richmeyer and Morton (1964). Essentially, it reduces the ordinary 
differential equation to a set of algebraic equations and these are 
then solved for f simultaneously with the boundary conditions. 

Instead of shooting out from the wall, the method proceeds inward, thus 
eliminating completely the initial guess on f"(x,0) which is required 
in the case of the Runge-Kutta technique. The shear stress is then 
obtained by differentiating f * numerically. The whole calculation is 
repeated until f and f f converge to the required accuracy. 

Initially, fqA = *0136 (obtained by setting q + = 3) is used. 
However, this gives a that is too high compared to the data of constant 
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pressure flow along convex wall. Various values of (^/a) are tried until 
one ( i -^/ a ) is found such that the calculated C f and the measured C f 
for the case of constant pressure flow along convex wall correlate with 
each other. Such a value of (i / A) is found to he .OUlT and this is 
used for the boundary layer predictions of the other two cases. 

4.6.1 Turbulent Boundary Layers along Convex Surfaces 

The results for the case of constant pressure flow are given in 

Figures 79 to 8l. The agreement among the various integral parameters 

'ft 

0,5 and H are very good (Figure 79)- In order to show more 

explicitily the curvature effect, two more calculations are made. One 

is without curvature , i.e. R -* oo , while the curvature for the second 

calculation is halved, i.e. R =20 inches. For both calculations, 

o 

the same initial condition and the same distribution of potential 

velocity at the wall are used. It can be seen that the result of convex 

* 

curvature is to reduce Q and increase 6 ) hence H is very much 

different from the corresponding flow along a flat plate. The calculation 
with no curvature shows that H = 1-35 when the flow is in equilibrium 
again after passing through the favorable pressure gradient. However, 
with the designed curvature (R q = 10 inches) H still keeps on 
increasing and this is an indication that the flow is not in equilibrium 
even at x = 75 inches. On the other hand, with the curvature halved, H 
approaches a constant at x = 75 inches. The prediction of C f is 
of course excellent, since i-J , A is chosen so that the calculated and 
measured matched each other at x = 71 inches, i.e. Station 11. 

For the corresponding flat plate flow , is about 13% higher than the 

measured C „ . This strongly indicates that turbulent flows along 
convex surface cannot support as high an adverse pressure gradient as 
the same flow over a flat surface. This fact will again be borne out 
in the separating flow case. Calculated velocity profiles at Stations 
2, 7, 9 and 11 are shown in Figure 80 together with the measured profiles. 
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In general, agreement is good, and the present method predicts the 
velocity pi of lie cpiite well even at Station 2 which locates downstream 
of the strong favourable pressure gradient. For purpose of clarity, 
the velocity profiles from the other two calculations are not shown. 

The present method also predicts the shear stress profile very 
well (Figure 8l), especially the point where the shear stress vanishes. 
Near the wall, the agreement is off, but in this region, the measurements 
arc in error because the rotating -wire method is not accurate in region 
of high shear as pointed out in Appendix D. Aside from this region 
which is about one displacement thickness, the agreement is very good. 

At Station 7j the shear stress profiles of the zero wall curvature 
(H 0 "*■ “)> ana half wa H curvature (R = 20 in.) are also shown. This 
shows clearly the effect of convex curvature in "cutting off" turbulence . 
With no curvature , the shear stress vanishes at the edge of the boundary 
layer (t) ~ 11), with = 20 inches, the shear stress vanishes at 
1 I 5 j while with R^ = 10 Inches (the designed curvature), the shear 
stress vanishes at T) ^ 4 . Hence it can he seen that even if the 
surface has a very small curvature, the point of zero shear will be 
somewhere Inside the boundary layer (the boundary layer is assumed to 
have been established previously). 

In the calculation of the separating flow case, instead of using 

the measured distribution of potential velocity at the wall all the way 

up to x = 75 inches, the extrapolated velocity distribution from x = 62.5 

inches to x = 75 inches is used (Figure 82). This gives a linear 

decelerating potential velocity distribution at the wall, and this is 

the desired velocity distribution to he set on the convex surface (see 

Section 2.3.2). The results are shown in Figures 83 to 86. As in the 

case of constant pressure flow, two more calculations (one with R -► 00 

o 

and the other with R^ =20 in.) are made for the sake of parametric 
study of the curvature effect. The calculated and measured 0 and 6* 
are given in Figure 83. Under the influence of adverse pressure gradient, 
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the differences in the development of 0 is unnoticeable . The 0 

development for the case of R q = 20 inches is practically the same 

as that of R -+ m . hence not shown on the graph. However, the 
o 

differences in &' are discemable, especially near separation. 

Although the location of the point of separation is not known, 
it can he said with confidence that the flow separates somewhere 
between x = 66 inches and x = 69.5 inches. If separation is defined 
as the point where C f -* 0 , then the separation point as calculated by 
the present method falls right within the range of x = 66 inches to x = 
69.5 inches, and is at x ~ 67 inches (Figure 84). With R q = 20 inches, 
separation occurs at x ^ 73 inches, and for zero wall curvature, the 
flow does not separate until x ~ 80 inches. This supports the previous 
conclusion that under the same distribution of potential velocity at 
the wall, the flow separates earlier when the surface has a convex 
curvature . 

The calculated C agrees well with the measured C f up to 
x = 59 inches, and after this the calculated value is generally lowered. 

The reason is that in the actual flow, the invisc id -viscous interaction 
causes the velocity distribution to level off, and therefore has a 
delaying effect on separation. This is borne out by the fact that if 
the measured velocity distribution is used instead of the linear 
decelerating velocity distribution, the present method gives good agreement 
with C„ up to x = 66 inches. However, separation is predicted to be 
at x = 78 inches, and no separation is predicted for the other two 
cases. This points out one of the difficulties in trying to predict 
near separation flow. Unless a way can be found to account for the 
invisc id -viscous interaction near separation, the difficulty remains. 

The agreement between calculated and measured velocity profiles at 
Stations 4, 7, 9 and 11 is excellent (Figure 85). However, the agreement 
is poor between calculated and measured shear stress profiles at 
Stations 4 and 7 (Figure 86). In spite of the poor agreement, the point 
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where x -* 0 is accurately predicted. Even though the flow is under 
the influence of a strong adverse pressure gradient, x still vanishes 
at about half the boundary layer thickness, and that by decreasing the 
curvature to half, it pushes the point where x *♦ 0 by about the 

sane amount as in tlje case of the constant pressure flow. This seems 
to indicate tliat curvature operates independently of the pressure gradient 
and this may also be the reason why the value of (£^/ A) > obtained by 
considering flat plate data, works so well in this case. The designed 
convex curvature not only causes to decrease by about 25$ (Figures 

79 and 84), it also causes the maximum shear stress reached inside 
the boundary layer to decrease by the same amount (Figure 86). 

In conclusion, the value of (fp/h) so obtained is also good for 
pressure gradient flow and that the present method can be considered 
rather successful in the prediction of boundary layer developments 
along convex walls with arbitrary pressure gradient. 

4.6.2 Turbulent Boundary Layers along Concave Surfaces 

The present experimental investigation has shown that the turbulent 
boundary layer on a concave surface is different from that on a convex 
surface because of the presence of a system of longitudinal vortices. 
However, it is most desirable, from the practical engineer's point of 
view, that turbulent boundary layers over concave surfaces can be predicted 
by simple method like that discussed above . This would be possible and 
therefore some approximations of the boundary layer growth could be 
obtained if the following assumptions were made about the flow. Firstly, 
assume the vortex pairs that consitute the vortex system are similar. 

(This can be evident from Figures 64 and 65). Secondly, assume the flow 
to be completely characterized by parameters which are obtained by 
averaging over anyone vortex pair. Finally, assume the eddy viscosity 
function given by equation (4.3-7) is capable of describing the resultant 
"average" flow. That this assumption is feasible can be seen from 


equation ( 4 . 3 - 7 ). For concave curvature, k is negative, so the factor 
within the curly brackets in equation ( 4 . 3 - 7 ) is always positive, in 
simple terms, this can be taken to characterize the increased mixing 
of flows over concave surfaces. 

It is only in this sense that the present calculation: arc compared 
with the measurements on the concave surface. Since only two measure- 
ments at two different z planes in a vortex pair are made, the averages 
of the various integral parameters can at best represent approximately 
the "average" flow. 

Having established the premise for the comparison between the 
present calculations and the measured data, the results are presented 
in Figures 87 to 90 - As a result of the vortex system, the Law of the 
Wall does not hold. Since direct measurements on the skin friction have 
not been made, the wall shear stresses at the four measuring stations 
are not known. Because of this, no comparisons between C f are made. 

The results for the velocity profiles are given in Figures 87 and 
89. Presented in Figure 87 is the results at the x = 70 inches position. 
The average values of 6* , 0 and H are included and these compared 
quite favourable with the calculated values . The shape of the profile 
is also quite similar. However, the calculated boundary layer thickness 
is not the same as the average boundary layer thickness. At the x = 96 
inches position, the calculated values of 6* and 0 are approximately 
greater than the mean values between Stations 4 and 5, but the shape 
factor H agrees to witttfn IDfa (Figure 89). This apparent difference in 

,y 

5 and 0 is probably due to the presence of more than one system of 
vortices at this position. 

The calculated shear stress profiles at x = 70 inches and 96 inches 
are given in Figures 88 and 90 respectively. The calculated shear stress 
at the wall at both these locations seem to agree well with the measured 
profiles, and at x = 70 inches the calculated profile and the average 
of the measured profiles correlated very well near the wall (Figure 88). 
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Even the maximum reached and the position where this mxiinum occurs are 


predicted quite correctly by the calculated profile. However, far away 

from the wall, substantial difference between the calculated and the 

average of the measured profiles begin to show. The difference increases 

towards the edge of the boundary layer, and become constant in the free 

stream. At x = 96 inches, the average of the measured profiles show 

a rather constant shear stress across the boundary layer, and towards 

the edge, the shear stress decreases steeply. However, the calculated 

profile does not display such a shape at all (Figure 90)* It has the 

same shape as that at x = TO inches with the maximum shear stress 

located away from the wall. (The appearance of the shear maximum away 

from the wall under a zero pressure gradient is due to the memory of 

the adverse pressure gradient which the flow has gone through at the 

entrance to the concave wall). Although the calculated value compares 

favourably with the average of the measured values near the wall, the 

calculated value is generally much lower than the average of the 

measured values far away from the wall. However, this difference 

decreases towards the edge of the boundary layer. Even then the 

calculated profile does not vanish but approaches a constant in the 

free stream. The fact that the calculated shear stress does not vanish 

can be accounted for by the shear stress relation (4.2-19). As explained 

_2ky 

in Section 4.2, t/p approaches the value of v( -2k U e ) in the 

free stream when v is taken to be v outside of the layer. This 

e 

is a small viscous stress, and in the actual flow, the shear stress would 

probably approach this value. On the other hand, in the present method, 

no provision is made to allow the eddy viscosity to slowly decrease to 

v , the molecular viscosity, as the edge of the boundary layer is reached. 

An a result, the value of approaches 2 v e ^ (see eq_. (4.5-13)). For 

flows along flat plates, the condition of zero vorticity in the free 

stream implies zero shear stress too. Therefore, it is immaterial what 

value v« takes. For flow along convex surfaces, the shear stress goes 
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negative somewhere inside the boundary layer, and although vanishing 

vorticity does not imply vanishing shear in the free stream, the 

eddy viscosity at the point where the shear stress goes negative is 

taken to be v so that the shear stress will again approach 
—2k y 

v( -2k e ) in the free stream. Hence, in this case adjusts 

itself. This is evident from equation (4.5-10), since the quantity 

inside the curly bracket decreases as r\ increases. However, the same 

quantity increases with r) if the curvature is concave, and this is 

the reason why the calculated shear stress remains large towards the 

edge of the boundary layer. This points to the need of modifying (i^A) 

in such a way that it will effect a rapid decrease in v towards the 

e 

edge of the boundary layer. 

It can be argued that the condition of small viscous shear can be 
satisfied by writing: 


T — ,au , „ x 

P e 1 


in accordance with Prandtl's (1929) mixing length argument. With this 
expression for (-uv), the eddy viscosity hypothesis becomes: 


T kU(c)U /dy + ku) 

i -/3 


(dU/fey - kU) 


3/2 


(dU/&y - ku Y 


SU/feyOu/fey + ku) 


Such an expression for v is undesirable because of the singular 

e 

behavior of v as 5b /dy 
e 


0 inside the boundary layer and 
(5U/fey + WJ) -+ 0 towards the free stream. Therefore, in spite of the 
fact that the condition of small viscous shear cannot be satisfied for 
the flow over concave walls, the present approach is adopted. 

In closing, it can be said that the present method can be quite 
useful in predicting the mean characteristics of flows over concave 
walls if (i^A) i s modified so that the condition of small viscous 
shear in the free stream is satisfied. 
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Chapter V 
CONCLUDING REMARKS 


5 .1 Conclusions 

As a result of this investigation, the following conclusions can 
be drawn. For the sake of clarity, attempts are made to divide the 
conclusions into the following three categories. The first two 
categories are concerned with the experimental investigation and 
the first one is limited to turbulent boundary layers along convex 
surfaces, while the second one dwells only on constant pressure turbu- 
lent boundary layers along concave surfaces. The third and final 
category includes those conclusions that are drawn as a result of the 
theoretical investigation. 


5.1.1 Turbul ent Boundary Layers along Convex Surfaces 

(i) In spite of the secondary flow which arises as a result of 
the longitudinal curvature of the test wall, a nearly two dimensional 
flow is established along the central plane of the convex test section. 
However, due to the rapid growth of the boundary layer thickness near 
separation in the case of separating flow, the secondary flow becomes 
quite significant. 

(ii) The Law of the Wall holds for turbulent flows along convex 
surfaces. It's validity can also be demonstrated through the method of 


matched asymptotic expansion applied to the curved turbulent boundary 
layer equations. The skin friction obtained from Clauser's plot of the 
velocity profiles correlates well with the momentum integral. 

(iii) A skin friction relation in which A is both a function of 
and 5/ R is obtained. The relation reduces back to the flat 

v t dx' 

plate skin friction relation as &/r -* 0 . It is applicable to turbulent 
boundary layers over plane or convex surfaces with arbitrary pressure 
gradient so long as A (&/R ? 



(iv) Measurements in the case of constant pressure flow show that 
the flow towards the end of the convex test section is not quite in 
equilibrium although separately the velocity profile and the shear 
stress profiles at the last two stations are quite similar. Also , the 
measurements indicate that even if equilibrium is reached, the defect 
part of the profile will be quite different from that of the corresponding 
flat plate profile . 

(v) The initial decrease in the intensities of turbulence is due 
partly to the favourable pressure gradient and partly to curvature. The 
inability of the intensities to increase further downstream is a clear 
indication of the effectiveness of convex curvature in reducing mixing 
between fluid layers. 

(vi) As a result of the reduced mixing, the boundary layer 
growth is retarded on convex surfaces. 

(vii) Also the ability of the flow to support adverse pressure 
gradient is reduced, hence , under the same wall static pressure distri- 
bution, the flow would separate earlier than the corresponding plane 
flow. 

(viii) Turbulence energy production is drastically reduced and is 
quite small at about half the boundary layer thickness. 

(lx) For the case with zero pressure gradient, the shear stress 
decreases steeply outside the viscous region and approaches zero at 
about half the boundary layer thickness. However, under the influence 
of strong adverse pressure gradient, the maximum shear stress occurs 
at some distance away from the wall, but it still goes to zero at 
about half the boundary layer thickness. 

(x) For the case with no pressure gradient, the shear correla- 
tion coefficient again remains constant for the greater part of the 
shear thickness, but drops to zero steeply towards the edge of the 
shear stress profile. 
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5.1.2 Turbulent Boundary Layers along Concave Surfaces 

(i) A system of longitudinal vortices similar to the Taylor- 
Gortler type vortices exists. For a given geometry of the concave surface 
and a given Reynolds number per unit length of the flow, the system of 
longitudinal vortices is stationary and does not depend on the upstream 
conditions . 

(ii) The limited data lends evidence to the fact that there are 
more than one system of vortices. At least two such systems are detected 
at stations 4 and 5; x = 9 6.0 inches. 

(iii) The boundary layer at the position of the trough of the wave 
is approximately twice as thick as the boundary layer at the position of 
the crest of the wave. Initial evidence also indicates that the growth 
rate of the boundary layer is the same in both these positions. 

(iv) Instability of the fluid particles, as a result of the 
concaye curvature, promotes mixing between fluid layers, hence a sub- 
stantial increase in the turbulence energy all across the boundary layer. 

(v) The shear stress profiles at stations 4 and 5, x = $ 6.0 inches 
show two peaks, one near the wall and one far away from the wall. The 
location of these peaks for the measurements at the position of the 

crest are different from that at the position of the trough. The 
existence of two peaks indicates the existence of two systems of 
vortices . 

(vi) The Reynolds stress uw is positive when measured at the 
position of the trough and negative when taken at the position of the 
crest. The opposite is true for vw . However, at x = 9^*0 inches, 
the distributions of uw and vw are definitely wave like and show 
two peaks. Again this indicates that there are at least two systems 
of longitudinal vortices at this stream position. 

(vii) Unlike the case of convex curvature, there is no region 
inside the boundary layer where the shear correlation coefficient is 
constant. This is another indication of the three dimensional nature 
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of the flow . 

(viii) For a concavely curved turbulent flow, turbulence energy 
production is not confined to the region very close to the wall, rather 
it extends nearly to the edge of the boundary layer. 

5.1.3 The optical Predictions 

(i) The set of equations (4.2-15) to (4.2-18) is a self consis- 
tent set of curved turbulent boundary layer equations. 

(ii) The eddy viscosity hypothesis as proposed by Mellor and 
Gibson is generalized to include the effect of curvature . 

(iii) The value .0417 for the ratio (^/a) is found to give 
excellent correlation for all three cases including the concave curvature 
case . 

( iv) The results of the present calculations supports the 
assumption that the value of (i^/A) obtained by just considering the Law 
of the Wall region can be used for the whole layer. 

(v) For convex surfaces, the present method predicts the 
boundary layer growth and the point of zero shear stress accurately. 

The present method also predicts the separation point fairly well. 

(vi) The present method can be used to predict the mean boundary 
layer development on a concave surface. 

(vii) For concave surfaces, because of the existence of the 
longitudinal vortices, the constant B in the Law of the Wall is no 
longer a constant. Although a certain log region still exists, the 
skin friction deduced differs considerably from that predicted by the 
present method. 

(viii) With convex curvature in the mean flow stream lines, the 
maximum reached by the shear stress is about 75 % of the corresponding 
flow over a flat surface. 
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^ . 2 Roconxxridatlonr; Tor Future Work 

Tar- following extensions of the present investigation arc 
rocorruncnclr.-d for future vork. 

(i) In view of the similarity of the equations that govern curved 
turbulent boundary layers and those tliat describe stratified flows, a 
frequency corresponding to the Brunt -Vaisala frequency may exist inside 
the boundary layer. Both Eskinazi and Yeh's findings and a preliminary 
measurements by the author do not seem to indicate the existence of such 
a frequency. However, due to the limited time available, this v/as not 
pursued in depth, but the author believes that this aspect of the flow 
should be investigated further. 

(ii) Efforts should be devoted to establish a formal analogy between 
the centrifugal force effect and the buoyancy effects in stratified flow. 

(iii) Spectral measurements should be made in order to better 
understand the detail structure of curved turbulent boundary layers, 
and the results compared with the data of Eskinazi and Yeh. 

(iv) The system of longitudinal vortices should be studied in 
more detail., so that its effect on the turbulence structure can be 
better understood. 

(v) Measurements of shear stress at the vail axe necessary for the 
clarification of the question of the validity of the Law of the Wall 

in concavely curved turbulent boundary layers . 

(vi) More measurements should be made in the transverse plane to 
see if similarity of the flow exists. Also, more measurements are 
required between the crest and the trough of the wave to establish the 
fact that in the mean, such a flow can still be treated assuming two- 
dimens ional ity . No doubt, this will be of value to the practical engineer. 

(vii) Parametric study of the curvature affect using the present 
calculation method will help improve design of such things as aerofoils, 
turbine blades, etc. 
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APPENDIX A 

CONSTRUCTION DETAILS OF THE CURVED WALL TUNNEL 

The curved wall tunnel used for the present investigations is an 
open -return, suction type wind tunnel. It is powered by a two speed, 

10 II. P. fan. The lower speed is capable of giving a flow with a Reynolds 
number per unit length of 4.37 x 10 ^ while the higher speed gives a 
corresponding Reynolds number of 7 *5^ x 10 . Since all the runs on both 
the convex and concave test walls are conducted with the lower Reynolds 
number, the higher speed capacity of the fan was never used. The tunnel 
consists of six different units, and these are; (l) the entrance section 
( 2 ) the contraction section, ( 3 ) the straight section, ( 4 ) the curved 
test section, (5) the exit section, and ( 6 ) the composite unit of the 
9CP bend, the diffuser, the fan and its housing. The tunnel is assembled 
by bolting these six units together as shown in Figure 1. With this 
arrangement, the curved test section can be changed with relative ease. 

The entire top wall and floor of the tunnel, from the inlet screen to 
the diffuser, is flat, thus giving the tunnel a rectangular cross-section. 
The depth of the tunnel is 4 feet and this gives a nearly two dimensional 
flow on the side walls of the tunnel. The quality of this two dimensional 
flow has already been commented on in both Sections 2.3 and 3.1 

The entrance section of the tunnel consists of five individual 
compartments 3 feet wide by 4 feet deep. Each has one layer of 2 l/2 
inches thick Hexcel aluminum honey-comb of l /4 inches cell size, and 
each compartment is separated by a single layer of l 8 -mesh screen. These 
compartments are bolted together with two more layers of screen, one in 
front of the first honey-comb and one downstream of the last honey-comb. 
This entire unit is then bolted on to the contraction section. 

The contraction section is designed to give a contraction ratio 
of 6:1. The side walls of this section are made of l /8 inches thick 
masonite, bent to the required curvature and kept in place by five ribs 
of 1 inch thick fir equally spaced on each wall. Tlie top and bottom 
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vails of the tunnel are made of l/2 inches thick paper faced fir plywood, 
and are “bolted on to the ribs on the side vails. Both side walls are 
lined with l/l6 inches thick vhitc panelite to give the required smooth 
surface . 

The straight section is 4 feet long and has a cross-section of 6 
inches by 48 inches. This cross-sectional geometry leads to nearly two 
dimensional flov everywhere on the two side vails. The side vails 
are made of l/k inches thick paper faced fir plyvood lined with l/l6 
inches thick white j>anelite. In order that the joint between the con- 
traction section and the straight section be smooth, the white panelitc 
lining is of one piece . This can be accomplished by cementing the lining 
on to the side walls after the straight section has been bolted on to the 
contraction section, thus making these two sections a composite unit. 

One side wall of the straight section is hinged to the contraction section, 
thus rendering it adjustable. The adjustment is l/2 inches at the end 
of the straight section. The original intent was to achieve added control 
on the free stream velocity distribution especially at the entrance to 
the convex test section. Later, it was found that the free stream 
velocity remained rather constant up to about 8 inches to the exit end 
of the straight section, and a flexible flap installed in the last 8 
inches will provide the necessary control. A diagram showing the flap 
and its control mechanism is given in Figure 31- With the installation 
of the flap, there is no need for the wall to be adjustable any more. 
Therefore, it is positioned to give a uniform width of 6 inches with the 
fixed wall and then securely clamped to the top and bottom wall of the 
straight section. Except for a section of 1 feet vide, l/2 inches thick 
paper faced fir plywoods are installed as the top and bottom walls of 
this section. The side walls are strengthened with five ribs of 1 inch 
thick fir equally spaced between the top and bottom to prevent them from 
wobbling. A l/2 inches thick by 1 feet wide lucitc is placed mid-way 
between the entrance and exit end of the straight section, so that flow 



observation can be made. 

Together two curved test sections are made, one has a convex test 
wall while the other has a concave test wall. The geometry of these 
curved walls have already been discussed in Section 2.1.2 and their speci- 
fications given in Figures 3 and 4. Hence, only the details of their 
construction will be described. 

The convex test wall is made of cold rolled steel 1/32 inches thick, 
with surface finish specified as 63 micro inches. The wall is rolled to 
the specified curvature and its shape is maintained by four ribs of 1 
inch thick fir also cut to the same curvature (see Figure 2). It is 
then bolted to the straight section and to the fixed wall of the exit 
section. The wall opposite the convex test wall is also made of the 
same material. Since it is necessary to be able to control the shape of 
the outer wall so that the desired pressure distribution on the convex 
wall can be obtained, three rolls of struts, with six in each roll, are 
installed. The linear adjustment is large enough to give both a constant 
velocity distribution and a linear decelerating velocity distribution on 
the convex wall. Details of the control mechanism is shown in Figure 91* 

As explained in Section 2.3.3, end wall jets are required for control 
of the secondary flow. Four tangential jets of l/8 inches I.D. are 
installed at the entrance to the curved test section on both the top and 
bottom walls of the tunnel. Two pressure regulators graduated to 1 psig 
are used to regulate the air supply which is from the shop’s main com- 
pressed air supply o The large fluctuation of the main compressed air 
supply made the installation of the pressure regulators necessary. The 
details of the end wall jets together with the side wall jet and the 
flap are shown in Figure 31* The side wall jet has already been dis- 
cussed in Section 2.2.3* Due to the presence of the side wall jet, the 
entrance end of the adjustable wall is secured to the straight section 
frame. The exit end is then screwed onto the adjustable wall of the exit 
section. 
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For no other reason "but the availability of material in the shop, 
the concave wall is rolled out of l /32 inches thick stainless steel with 
a surface finish specified to 32 micro inches. The shape of the wall is 
maintained by two ribs of l/2 inches thick aluminum alloy spaced 28 
inches apart. The opposite wall is made adjustable through two rolls 
of struts, with nine struts in each roll. It is bolted to the straight 
section and the fixed wall of the exit section. The details of the 
adjustable inner wall is shown in Figure 92. The same number of tangential 
wall jets are installed at the entrance to the curved section. However, 
no side wall jet is required because the inner wall is adjusted to mini- 
mize flow deceleration at the entrance to the concave test section. 

Because of this, the entrance end of the adjustable inner wall can be 
bolted to the straight section. Again, the exit end is screwed onto 
the adjustable wall of the exit section. 

The bottom wall of the curved test section is the base which 
also serves as the support for the tunnel, while the top wall is made 
of l/2 inches thick lucite. The lucite top is not clamped or bolted 
onto the side walls. The pressure difference between the ambient air 
and the flowing stream will press the top onto the side walls. Leakage 
can be prevented when rubber seals are cemented onto the side walls as 
shown in Figures 91 and 92. 

The exit section is so constructed that one adjustable wall is 
all that is required to accommodate both the convex and concave test 
sections . Then come the 90° bend, exit diffuser and the fan and housing 
unit which are all bolted together to form one composite unit. The 
support for the whole tunnel can be adjusted to give a level bottom wall. 
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APPENDIX B 


POSITIONING AND ACCURACY OF TIE VARIOUS PROBES 
1 ) Conrad Probe 

The Conrad probe is used as a nulled direction probe. The null 
ancle can be found by aliening the probe in the direction of flow in the 
pipe flow system (Figure 93 )* The Conrad probe is carried by a rotary 
probe carrier (Figure 9 ) which in turn is fixed to the probe carrier 
(Figure 10). The whole unit is then mounted on a horizontal platform, 
and the Conrad probe introduced into the pipe flow horizontally. Every - 
tjjne the Conrad probe is introduced into the rotary probe carrier, it 
is set into the same position by a set screw with a conical tip that 
fits right into the conical dent in the probe . This arrangement allows 
a reference for the null angle measurement to be established, and does 
not depend on the relative location of the whole unit as long as the 
probe carrier is always mounted horizontally. Once the null angle is 
found, the whole unit is then transferred to the curved wall tunnel. 

The Conrad probe is introduced into the test section from the wall 
opposite the test wall, and the probe carrier is mounted as shown in 
Figure 10. In order that the axis of the probe be normal to the wall 
where measurements are to be taken, radial lines are inscribed on the 
bottom wall and the lucite top of the tunnel. The probe axis is then 
aligned with these radial lines by sighting with a telescope placed on 
the lucite top. The Conrad probe is shaped like a goose-neck (Figure 8) 
so that the edge of the probe is aligned with the axis of the probe. 

Hence, once the axis of the probe is aligned with the radial lines, the 
edge of the probe will traverse along the normal to the wall. Since 
both the curved wall and the Conrad probe are made of metal, an electrical 
circuit can be so arranged that when the probe is touching the curved 
wall, a light goes on. When the probe is brought back to such a position 
that the light just goes off, the leading on the vernier mounted on the 
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probe carrier is noted. This gives the zero reading for y . This way 
of positioning the probe has an accuracy in y of the order of + 0.001 
inches. Since the traverse is towards the free stream, and because of 
the small uncertainty in the measurement of y , the measurements near 
the vail will be much less accurate than those far avay from the wall. 

The rotation of the Conrad probe is controlled by a set of gears so 
arranged that the smallest angle which can be measured accurately on the 
dial is one tenth of a degree. 

The reason for choosing this particular probe (see Figure 7 for 
probe geometry) for the yaw measurements is because of its high sensiti- 
vity and its zero scale effect (see Bryer, Walshe and Garner 1958) . 

The pressure difference between the two openings is measured by a Pace 
Model CPflD- .1PSID pressure transducer and the :read out is on a DI8A 
Digital Voltmeter. Because of the pressure fluctuations, the Conrad 
probe can be rotated to read to +.01 volts only. The sensitivity of 
the probe decreases with the free stream velocity (see Figure 7), 
therefore the accuracy of the yaw measurement decreases as the wall is 
approached. To determine the effect of the free stream velocity, the 
null angle of the probe is measured for two different velocities in 
the pipe flow. For a velocity of 72.6 ft/sec., it was found possible to 
determine the null angle to + .2° , and when the velocity is reduced to 
28.9 ft/sec., it was only possible to determine the null angle to + 0.5° • 
Even though the turbulent fluctuating velocities in the curved tunnel 
are not the same as that in the pipe flow, much the same accuracy would 
be expected when the Conrad probe is used to measure the secondary flow. 

2) Total Head-Probe 

The total-head probe is also shaped like a goose-neck with the 
edge of the probe aligned with the axis of the probe (Figure 7). With 
this configuration, the total-head probe can be positioned in the tunnel 
in the same manner as the Conrad probe, therefore, the same accuracy 
in the measurement of y . 



The major source of errors of a total-head probe used for measure- 
ment on a low velocity stream are the effect of turbulence , the effect 
of yawing and the effect of a wall. The errors duo to the effect of 
turbulence on total -pres sure readings are not well understood. However, 
assuming the frontal part of the probe to be a stagnation point, 
Goldstein (1936) found theoretically that the total -head probe measured 
the total head of the total velocity vector, i.e., 

Pt ± = P s + i p{(U + uf + v 2 + w 2 ) 


This expression is correct only if the frontal area can be considered a 
true point. However, due to the finite size of the total-head-tube 
hole, deviations from the above expression may be expected. Aside from 
the lateral velocity gradient effect, which is neglected in the above 
expression, there is also the lateral velocity fluctuations effect. 

This may not produce an impact pressure as given by Goldsteins expression, 
but rather appreciably smaller. Hinze and Van der Hegge Zijnen (1949) 
neglected the effect of the lateral turbulence velocities and arrived at 

p t . = p s + i pu 2 + i P u 2 


if 


is the true total pressure, then 


p *i 


1 c 

& pu 



Hence, for a turbulence level less than 10$, the error in total -pres sure 
measurement will be less than 1$ of the dynamic pressure. For the 
present investigations, the turbulence level in the flow along convex 
surfaces is always less than 10$, therefore the readings of the total- 
head probe are not corrected for turbulence effect. As for the velocity 
profile measurements on concave surfaces, hot-wires are used, thus 
eliminating the need to correct for turbulence effect, since in this 
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c?if'.e, the turbulence .level is expected to be higher than the corresponding 
flow along convex surface. 

Gracey, et.al. ( 1951 ) found that moat simple total-head probes arc 
rather insensitive to yaw. They presented their data in terms of the 
"critical angle" of a given probe. The "critical angle" is defined as 
that angle at which the error in reading the total pressure amounts to 
0 f fhe indicated dynamic head. They found that for all simple probes, 
ranging from cylindrical heads to ogive heads, the critical angles arc 
of the order of + 15 ° in subsonic flow. Rogers and Berry ( 1950 ) also 
found that the critical angles for a round nose probe and a flat nose 
probe are of the order of + 15 ° • Therefore, the flat nose total-head 
probe used in the present investigation is quite insensitive to yaw. 
Although there is a certain amount of secondary flow in the test section 
due to the longitudinal curvature of the test wall, the flow in the 
plane of symmetry is not expected to deviate too much away from the 
stream direction. Yaw measurements show that the angle of deviation is 
less than a degree. Again, there is no need to correct the total pressure 
measurements to account for probable yaw of the flow. 

Very little work has been done to investigate the probe-wall effect. 
By observing the response of a one-sided rectangular tube with the wall 
as the bottom surface, in a known laminar profile, Stanton, et.al. (1920) 
were able to plot "effective position of tube" against "opening of tube" 
and used this plot as a calibration curve in interpreting measurements 
near the wall in turbulent profiles. The correction for wall effect 
was found to be in the form of an "effective displacement" of the probe 
centre. Corrections to the measurements near the mil are made using the 
curve of Stanton, et.al., but the corrections are found to be negligible. 
In using the curve of Stanton, et.al., the assumption was made that the 
effective displacement was equal in the calibrated laminar boundary 
layer and in the measured turbulent layer. However, this assumption 
does not hold even in the viscous sub-layer of the turbulent boundary 
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layer. As a result, the wall effect is not fully accounted for by 
the curve of Stanton, et.al. on top of thin, the transverse velocity 
gradient effect id also important. These two effects together would 
account for the observed shifts of the measurements near the wall. 

3 ) Hot -W ire Probes 

The hot-wire probe is introduced into the tunnel much the same way 
as the Conrad probe and the total -head probe. However, in this case, 
the distance between the hot-wire and the wall cannot be measured by 
the same method. This is due to the fact that the hot-wire probe, 
unlike the Conrad probe or the total -head probe, is very fragile, hence 
the hot-wire will break once it touches the wall. To remedy this, a less 
accurate method for determining y is used. 

A broken hot-wire probe is used as a dummy probe. The length 
(from the tip of the prongs to the edge of the adaptor) of this dummy 
probe is measured accurately (to one tenth of one thousand of an inch) 
using a telescope. The length (from the hot-wire to the edge of the 
adaptor) of the hot-wire probe is also measured to the same accuracy. 
First the dummy probe is used and the probe is traversed inwards towards 
the test wall until the tip of the prongs touches the wall. This closes 
the electi ic circuit and the light goes on. The probe is then brought 
back until the light just goes off. The reading on the vernier in the 
probe carrier is noted. The whole probe is then retracted until it is 
outside of the tunnel. The dummy probe is removed from the adaptor 
and the hot-wire probe is put in its place. Knowing the backlash of 
the probe carrier, the length of the dummy probe, the length of the 
hot-wire probe and the vernier reading at the point where the dummy 
probe just touches the wall, the hot-wire probe can be set at any distance 
away from the wall by moving the probe towards the test wall, v/ith this 
arrangement, the boundary layer is traversed inward from the edge and 
not outward from the wall. 

The miniature TGI hot-wire probe has a mark on it so that every 
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time it in snapped into the adaptor , good contact is ensured only when 
the mark on the probe is flushed with the edge of the adaptor. This 
also ensures that the length between the hot-wire and the edge of the 
adaptor is always the same. This fact is borne out by repeated measure- 
ment of the length when the miniature hot-wire probe is snapped in and 
out of the adaptor. The naxinuun variation between five different 
measurements is less than 0.001" . The backlash of the probe carrier 
is measured by a dial gauge to 0.001". Repeated measurements show that 
the backlash varies between 0.006" and 0.010". The backlash was also 
checked from time to time and in no case was it found to exceed 0.010". 

A mean value of 0.008" is used for all hot-wire measurements. 

With this arrangement , it is possible to determine y to an accuracy 
of + 0.008". Since the accuracy of the hot-wire measurement near the wall 
depends on the size of the wire other than the velocity gradient , heat 
transfer to the vail, etc., it is not very meaningful to try to measure 
the turbulence velocities any closer to the wall than the length of the 
hot-wire. The miniature TSI hot-wire has a sensitive length of 0.050", 
therefore an error in y of + 0.008" will not contribute very much to 
the overall error of the turbulence velocities measurements. Besides, 
the measurements at Station 1 serve as a check not only in the two 
dimensionality of the flow, but also in the accuracy of the above 
arrangrnennt for the determination of the position of the probe away from 
the wall. 

The various effects associated with the rotating-vire method that 
affect the accuracy of the hot-wire measurements have been discussed by 
Bissonnette (1970), therefore they will not be repeated here. As for 
the reliability and the overall accuracy of the rotating -wire method, 
Appendix C discusses this in some detail. 
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APPENDIX C 


FULLY DEVELOPED TURBULENT PIPE FLOW EXPERII-TENT 


In order to test the accuracy of the rotating -wire method, a 
fully developed pipe flow experiment is set up. The pipe assembly lias 
a convergent nozzle at the entrance. This is then followed by a diffusor 
and two layers of 2 l/2 inches thick Hexcel aluminum honey-comb flow 
straighteners of l/2 inches cell size. The flow straighteners are pro- 
ceeded and followed by two l8-mesh screens. This is followed by two 
18 feet long aluminum pipes of 6.005 inches I.D. The working section 
is downstream of the aluminum pipes and is made up of one 3 feet long 
lucite tubing which is machine bored to 6.005 inches I.D. Downstream 
of the lucite tubing is another aluminum pipe of 3 feet long. In order 
to minimize vibration in the flow in the working section, the exit diffuser 
fan housing unit is connected to the aluminum pipe by a bellow. The pipe 
flow tunnel is of the open-return, suction type and the air is drawn in 
by an axial flow fan driven by a variable speed d.c. motor. A picture 
of the pipe flow tunnel together with the hot-wire equipment is presented 
in Figure 93 • 

Fully developed turbulent flow is obtained at about 60 diameters 
downstream of the entrance. The axial pressure gradient for the next l8 
diameters is plotted in Figure 9 U. The friction velocity can be calcu- 
lated from the measured pressure gradient and is given by: 


u 


r . 

1 


P 


dp 

dx 


(C-l) 


dp 


where is the inner radius of the pipe and ^ is the axial pressure 

gradient. If the flow is fully developed turbulent flow, the Reynolds 
stress uv is related to the friction velocity u T by the expression: 


uv 


2 , r , 

u t (— ) 


(C-2) 


However, very close to the wall, this does not hold. 


108 



The measured uv is compelled- with "the iheoreiical prediction as 
given by equations (c-l) and (C-2). The comparison is shown in Figure 95 
and the general agreement is good. The measured normal stresses u , v<- 
and w 2 are shown in Figures 9 6 , 97 and. 98* Together with the present 
data is also plotted the data of Laufer (195*0, Sandborn (1955), and 
Patel (1963 and 1968) for comparison purposes. It is seen that the 
present data falls well within the evident scatter between the various 
sets of data* 

If the flow were truly axisymmetric and fully developed , then the 
off-diagonal stress components uv and vw are zero. The values uv 
and vw as measured by the rotating -wire method are less than % of 
the local u 2 • Also, these values are randomly distributed across the 

pipe . 

The above comparison shows that the rotating -wire method is rather 
reliable. It also indicates that, in similar flow conditions, the 
rotating-vire method will permit the determination of each component of 
the Reynolds stress tensor to about % of the locally measured turbulence 
energy per unit mass. This estimate will no longer be true in regions 
of very high shear such as boundary layer measurements very close to 

the wall. 

The above work on the verification of the rotating-wire method was 
performed in cooperation with Mr. Luc Bissonnette, who is also a graduate 
student in the Department of Aerospace and Mechanical Sciences. 



APPENDIX D 

CURVED BOUNDARY LAYER EQUATIONS 
THROUGH METHOD OF MATCHED ASYMPTOTIC EXPANSIONS 

Equations (4.2-7) to (4.2-9) are made non-dimensional by dividing 

all velocities by U q , all coordinates and the radius of curvature R 

by £ and the pressure by pU . The Reynolds stresses are made non- 
2 ° 

dimensional by u t where u t can be any velocity characteristic of 

the turbulence field. After some algebra the equations can be written 
as : 
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(D-M 


t = 1 Sv _ cRJ K 

^ 1+KY oX ’ dY 1+KY 

vhere T. . are the non-dimensional components of the Reynolds stress 
J 

tensor and 



Mellor (1970) has shown that — 0 as e -» 0 for arbitrary n . 

e 11 

Therefore, the technique of matched asymptotic expansion in terms of e 
is applied to equations (D-l) to (d- 4) for the sole purpose of seeking 
a set of self-consistent second order boundary layer equations. 

Due to the presence of the two small parameters e and £ , three 
regions exist in the flow field and these are: (l) the Outer or Inviscid 
Region where the length scale is i , (2) the Middle or Defect Region 

where the characteristic length is a = ei , and (3) the Inner or 

T> 

Viscous Region vhere the length scale is £& = £e£ . Since the radius 

of curvature is assumed to he of the same order as i , K(x) = 0(l) . 

Mellor 1 s (1970) procedure is followed closely in the subsequent derivation 
of the boundary layer equations. 

Outer or Inviscid Layer 

In this region the free stream turbulence and the free stream 

vorticity are assumed to be zero, hence T. = 0 . t = 0 and the following 

ij b 

expansions are assumed 

U = U^CXjY) + eU 2 (X,Y) + e^X^Y) + 

V = V-jUjY) + eV 2 (X,Y) + e 2 V 3 (X,Y) + (D-5) 

P = P-l(X,Y) + eP 2 (X,Y) + e 2 P 3 (X,Y) + 

Substitute (E~5) into equations (D-l) to (D-3) and collecting terms, 
the following is obtained, 
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to first order 
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1 SP 1 


1+KY BX 


(d-6) 


1 Bv x Bv-l 
1+KY U 1 Bx” + V 1 bY~ 


-JL-u 2 

1+KY 1 


Bp 3 

By 


to second order 


Bu 2 * 

ST* Si (1+ct)v 2 ' 0 


i Bu 2 Bu x Bu 2 Bui v 

1+ky < u i b 3T + u 2 b~ 3 + ( v i bT” + v 2 bF 3 + ife (W¥i) 

i Sp 2 


1+KY BX 


(U-, 


Bvc 


bV, 


BVi 


Bv- 


1+KY IU 1 Bx + u 2 Bx 3 + fv 2 By + v i By 3 


k 


1+KY a i U 2 


BP2 

By" 


(D-7) 


to third order 


BUo * 

Sx + 37 ((ltIOf)V 3> - 0 


2 BUo Bu 2 Bu x BUo 

1+ky * u i Bx” + u 2 Bx - + U 3 bF 3 + ( v i By" 


Bu 2 Bu-, K 
v 2 3— + V 3 3— 5 + I7ky '“1-3 


(u n v. + U 2 V 2 + U^J 


1+KY BX 


Bv 3 Bv 2 Bv x Bv 3 BVo Bv-, 

1+KY IU 1 Bx" + u 2 Bx” + u 3 bF 3 + {V 1 By” + v 2 By - + v 3 5F 3 


(u-. 


K 2 ^ p 3 

lfet U 2 + 2U X u 3 }= - w - 


(D-8) 
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Middle or Defect Layer 


In this region, the length scale is A t = e£ , therefore Y = ey 
where y is the normal coordinate in this region. In addition, the 
following expansions are assumed 

2 

U = u (X,y) + eu 2 (X,y) + e u^(x,y) + 

2 

V = e{v (X,y) + ev 2 (X,y) + e v^(X,y) + ) 

2 (D " 9) 

P = P x (X,y) + eP 2 (X,y) + e P^(X,y) + 

T. . = t .(X,y) + et .(X,y) + e 2 t (X,y) + 

10 lij 


Hence 

1 + KY = 1 + eKy 

and 

(1 + eKy) -1 = 1 - eKy + £ 2 (Ky) 2 - e^(Ky)^ + e (Ky) 

Rewrite equations (D-l) to (D-3) i n terms of X and y . Substitute 
(D”9) into the resultant equations and collecting terms, (since e - 
0(e n ), all terms that have e as coefficient can he neglected) the 
following is obtained; 

to first order 

Suq dv-j^ 

ST + By” 

Su l Bui _ Spi ( D _io) 

u i SF V 1 ay " " sx 

bvi 
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to second order 


du 2 Bv 2 SCKyv-^) 
bX by + by 


u 


3^2 0U2 u U 1 ou l 




Su-i 


1 bX 


+ v. 


+ u 


1 by 2 bX 1 by 


dp 2 St-, 

KV 1 (U 1 + * sr> - - ST * 17 2 


bv. 


Ku n 


1 by 

to third order 

d(Kyv 2 ) 

bX by + by 


0 


du- 




3u-, 


5u n 


c)u^ 


3 ^“3 UUL 1 UU 1 ou 2 

A 1 + V 1 ar + u 3 aF + v 3 sT + u 2 


du 2 

+ V 2 ST + K U 2 T 1 + K “l v 2 + 


Ky (v - 1 + v, 


ldy tv 25 


du, dt-, 

1 lx^ 


Sy 


} = - 


h l2 

bX 


^ t lxx ^ t 2 xy 
+ 37— + 3—^ + 2 Kt 
oX ay lay 

SV 1 SV 1 dp_ dt_ 

U 1 ST 4 T 1 sr - K u i< 2 V ■ - ST + sr* 


Inner or Viscous Layer 
In this 
Y = eey , and 


In this region, the length scale is eA } hence y = ey 

t 


— = 1 - ee (Ky) + (el') 2 (Ky) 2 


1+KY 


(P- 11 ) 


(P-12) 
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The following expansion is assumed; 


U = 

V = 

P = 

T. . = 
3-J 


U-p^y) + eu 2 (x,y) + e 2 u (x,y) + 

ee{v i (x,y) + ev 2 (x,y) + v (x,y) + ) 

P x (Xy) + g P 2 ( X , y ) + e 2 p (X,y) + 

W x,x) + e W x,x) + e2 w x ^ } + •••• 


(D-13) 


Rewrite the equations (D-l) to (D-3) in inner variables. Substitute (D-13) 
into the resultant equations and collect terms, the following is obtained: 

to first order 


bX + by 

0 


0 


0 

*2- 
o Uj_ 



Spi 


to second order 


du£ 

SX + by 


0 


0 


0 


A * 

ay hxy + ay •* 
dp 2 
By 


to third order 


Sun bvj 

ET + sF 


0 


(D-14) 


(D-15) 
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0 



8G3 

8y 1 2xy 

+ 8y 


Bt iyy 

, 

By 

8y 


The aim here is to obtain a set of first order curved boundary 
layer equations, therefore the solutions of these various equations will 
not be discussed. Instead, the reader is referred to the paper by 
Mellor (1970). The boundary conditions for the three layers are very 
well discussed in Mellor's (1970) article and these will not be repeated 
here. However, it should be pointed out that the conditions of zero 
vorticity in the free stream is satisfied by matching the vorticity in 
the middle layer as y -» 00 to that of the outer layer as Y -* 0 . 

This gives: 


to first order 


Lim 

y -*■ » 


SUf 

8y 


0 


(D-17) 


to second order 


<$U ? / , o \ 

Lim {- r — - K u } = 0 (D-l8) 

dy 1 

y -* » 

The set of second order boundary layer equations is given by the sets 
(D-ll) and (D-15). The free stream velocity is given by (D-17) and (D-l8). 
Since the first order equations in the Middle and Inner layers are the 
same as that for a flat plate (see Mellor 1970 ), the solutions to u^ 
and v are : 


u = U x (X,0) 


v 


1 


= - y u^x^) 
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With these substitutions, the boundary layer equations become 


Bu^ c)v 2 

bx~ + by 



5 u 2 

u i ST * (u 2 


5u 2 dUj 

y sT ) 5T 



dU x 

1 dF” 


} y + 


— t 1 ) 

Sy V 


(D-19) 


(D-20) 


p dl>2 

KUi(X,0) = 


(D-21) 


where 


T_ 

P 


Su 2 

= - UV + V -7 

py 


and 


y -> oo 



+ KU 1 ) 


0 


The rather surprising result is that equations (D-19), (D-20) 
and (D-21) are linear. It is believed that they are the self consistent 
turbulent boundary layer equations (also for the plane case obtained by 
simply setting K = 0). However, turbulent boundary layer researchers 
have long been accustomed to making some of the approximations involved 
in these equations (and justified consistently in a laminar like boundary 
layer approximation), but not others. Notably the non-linear terms in 
the x-momentum equation are retained. For the present working equations, 
the same practice is adopted here. This is to conform to convention and 
hopefully to improve accuracy although in the light of the present 
asymptotic analysis this practice means that some higher order terms are 
being reinstated while other terms of the same higher order are being 
left out. Thus, reinstating some terms found in the original equations. 
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the set of boundary layer equations become: 

(D-22) 
(D-23) 


k U 2 = 

P P dy 

(D-24) 

Lim + kU) = 0 

(D-25) 


y - oo 

Here the full non-linear advective terms in the x-momentum equation 

are reinstated. Also, in order to satisfy equation (D-25), i.e. 

-ky 

U ~ Up = e as y oo , it is necessary to include both viscous 
terms in equation (D-23) where 

i = v e ( |^ " kU) ( D - 26 ) 

In this way the combined viscous term will limit to zero and the outer 

boundary condition can be satisfied even if approaches constant as 

y 00 • The constant may be the actual molecular viscosity or it can 

be the eddy viscosity artificially maintained constant for large y . 

In the case of plane flow it is well established that whether v main- 

e 

tains a constant value or decreases for large y makes little difference 
in the results. 

p 

As explained in the text, it is necessary to keep the term kU 
on the left of equation (D-2^) rather than replace it with kU 2 , for 
example. Otherwise, the pressure term will not quite balance out the 
advective terms and the integral of equation (D-23) will not exist. This 
is not only uncomfortable with regard to the von Karman integral, but is 
equally uncomfortable in the process of numerical integration of the 
full equations* 


5u 5v 


+ c— + k V = 0 


x by 




- - jr - I- t 1 ) ♦ 1 

p d* dy v p y p 
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APPENDIX E 


DERIVATION OF TIE TURBULENCE ENERGY EQUATIONS FOR 
CURVED BOUNDARY LAYERS 

2 “o 

The equations for u , v and v can be obtained by nultiplying 
equation (4.2-11) by u , equation (4.2-12) by v and equation (4.2-13) 
by v . Similarly the equation for uv is obtained by adding v times 
equation (4.2-11) to u times equation (4.2-12). Making use of both 
the continuity equations for the mean and fluctuating quantities, and 
taking the time average of the equations, the results are: 


5u 2 i a_ 

bt 1+ky bx 



+ u^ + 2u(p~x )] 
v xx 


+ (u 2 V + 
by 


2 

u v 


- 2u 


T ) 

xy 


+ gl [u " 


2 u T 


xz 


} 


+ — {2uvU + 3 u 2 V + 3 U 2 v} = 
1+ky 


2 , 2u 2 5U 

1+ky dx 


— 5U 2 Su 
2 uv < — + rr ~, — P 


2 Su 

T 


dy 1+ky bx 1+ky xx <3x 


au 4 k — au 

2 x - — + — ux - 2 x c — 

xy by 1+ky xy xz ^z 


(E-l) 



1+ky bx 


{■v^U 


+ u 


v 2 


2 VT } 

xy 


+ {v 2 V + + 2v( p 

oy 


T 

yy 



fv^v - 2vx ) 

yz 


k 

1+ky 


{4 uv U 


2 

v 


V + 2 


2 

u v 



2uv bV 
1+ky bx 


(E-2) 


0 -2 bV 
- 2 v — + 

Sy 


2p 


dv 

'by 


2t 

xy <3v 

1+ky bx 


- 2 


bv 
r yy by 


- 2 


bv 
yz b Z 


2k J 
Vt T 

1+ky xx 


- T ) 

yy 
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dw 1 d , 2 . 2 


bt 1+ky bx 


{w U + uv - 2 v t ] + 
xz 


^ {V2V + ^ + ^ ^ + ^P-T^J} 


ZZ 


2r 


+ r~r — (v 2 V + w 2 } = 2p 3w XZ 5w 


1+ky 


dz 1+ky dx 


(E-3) 


o dw 2k „ dw 

- 2 x -r — + - — - — w T - 2 x r — 
yz dy 1+ky yz zz dz 


duv 1 d , — . 


d 2 _ 


St ' 1+ky dx (uvU + U v ' u Txy + v(p- Txx )} + — (uw + uv - vx xy 


+ u(p-T )) + r- {uvw - UT - VT } 

yy oz yz xz J 

k. 2 2 2 o — 2 

+ 1+kv -2u ) u + u v - u + uv + 2uvV) = 


u(p- 

T 

yy 

k 

1+ky 

(0 

UV 

dU 

1+ky 

bX 

T xy 

bu. 

1+ky 

bx 


uv dU 2 dU du 1 dv 
V dy + ^ + 1+ky 


2 . 
u dV 


1+ky dX 


dV 

UV dT 

(E-4) 


du du T xx dv 

T w -n T \ 


bv 


1+ky dX yy dy T ^dz 1+ky dx ^dy 

k 


dv 2k — 

xz dz 1+ky VT xy 1+ky l xx _l yy 


u ( t -t ) 


In accordance with the mean turbulence field, closure scheme of 

Mellor and Herring ( 1S*T0 ) _? assumptions are made to relate all terms in 

equations (E-l) to (E-k) in terms of u . d . . For the sake of completeness 

J 

these assumptions are: 


U.T, . + U T . = v { . 
i kj j ki dx 


du u du u du .u 
J k A , i j 


dx. 


(E-5) 
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\j Sx + T ki ox- k 


Su J . *a\. 


3 A 


(E-6) 


du i 1 q 


— * {1 j-Tr - -4 4) 

dx. 6 in 1 J 3 
J 


pu. 


3 Sx. 


5u^u k 5^1% ^ u i u j 


u.u.u, = - q£_ (~ rr ^ + -x + J 


*i”j _ k ^ '3 1 ax. ax. 


ax,. 


(E-T) 

(E-8) 

(E-9) 


With these 
reduce to 


s implications and after much algebra, equations (E-l) to (E-4) 


DU 2 la _1 2 ja : \ 

= - 3 a ' 3 A 1 +ky ax 
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3 


-2^4- 44 + u2 v} 


by 1+ky 

i a 
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a 


2 2 
x i au 2 j_ aq, 
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(e-ii) 


+ - 3 ^L )} 

by By 
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(E-13) 


Invoking the two dimensional boundary layer approximations, and 
in addition, the boundary layer thickness is assumed to be very much 
smaller than the radius of curvature of the surface, such that 6/k(x) « 
then the coefficient 1 + ky can be approximated by 1 . with these 


12 ? 



iraplif icationo, equations (E-10) to (E-13) reduce to 
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figure 1. Curved wall tunnel. 
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FIGURE 3 ££M£1RV OF COAfl/fX JEST SECJIOhL (scale 



Geometry of concave test section. (Gcale 1" = 4 ") 
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Figure 4. 




STATIONS 3 AND 10 ARE 14 (IN) ABOVE PLANE OF SYMMETRY 
STATIONS 5 AND 12 ARE 14 (IN) BELOW PLANE OF SYMMETRY 

FIG. 5a LOCATION OF MEASURING STATIONS FOR CONSTANT PRESSURE 
FLOW ALONG CONVEX WALL 



STATIONS 3 AND 8 ARE 14 (IN) ABOVE PLANE OF SYMMETRY 
STATIONS 5 AND 10 ARE 14 (IN) BELOW PLANE OF SYMMETRY 

FIG. 5b. LOCATION OF MEASURING STATIONS FOR SEPARATING FLOW 
ALONG CONVEX WALL 
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FIG. 6. LOCATION OF MEASURING STATIONS FOR CONSTANT PRESSURE 
FLOW ALONG CONCAVE WALL 










Total -head probe and Conrad probe 
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Figure 9 * Rotary probe carrier 




} 


7 f; 


Photograph of hot-wire equipment. 
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ROOM TEMPERATURE RE 
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FIG. 14 

SAMPLE SET OF CALIBRATION CURVES FOR CONSTANT TEMPERAT 
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C 13 r> -1 *~~ c ■■ 

FIGURE lb 

DlAG&AttAllC LAYOUT OF HOJ-WlUE 

£Q U/PMEUJ 
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A SAMPLE CALIBRATION CURVE FOR LINEARIZER OUTPUT 

FIG. 17 


1M 


STRAIGHT - WIRE 


u 


SLANTED - WIRE 





(p) PROJECTION OF SENSOR IN PLANE XOZ 
(m) NORMAL TO LINE (p) IN PLANE XOZ 

FIG 18 . COORDINATE SYSTEM FOR THE ROTATING STRAIGHT -WIRE AND 
SLANTED -WIRE PROBE 
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PROBE STRUT CROSS-SECTION' 


SCALE : 2 x FULL SIZE 


PROBE CARRIER 

FIGURE 22 CROSS ~ SECT! 0 UAL DRAIV/UG OF STRUT 
AUD PROBE CARRIER 
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DEFECT PLOT OF VELOCITY PROFILES AT STATION 
NO. I, X (IN) = 24. 
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FLAT PLATE DATA, KLEBANOFF (1955) 
CONSTANT PRESSURE FLOW ALONG 
CONVEX WALL 
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FIG. 26 





FIG. 27 DISTRIBUTION OF u'/U DW AT STATION NO. I, X(IN)=24. 
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Figure 31. Blowing and flap control mechanisms. 
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WALL STATIC PRESSURE DISTRIBUTION FOR CONSTANT PRESSURE FLOW ALONG 
CONCAVE WALL 
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& °= tan'(^f-) 


FIGURE 34- YAW DISTRIBUTION FOR CONSUUT 

PRESSURE FLOW ALOUG CONl/FX 
WALL 
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(N,)A Aft 5 ’ 2 ON NOI1V1S (NI)X 



FIG. 35. YAW DISTRIBUTION FOR SEPARATING FLOW ALONG CONVEX WALL 
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FIG. 39 VELOCITY PROFILES FOR CONSTANT PRESSURE FLOW ALONG CONVEX WALL 
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PIG 41 SEMI -LOG PLOT OF VELOCITY PROFILES TO SHOW THE LAW OF THE 
UUAI I rfoion FOR CONSTANT PRESSURE FLOW ALONG CONVEX WALL. 








FIG. 43 DEFECT PLOT OF VELOCITY PROFILES FOR CONSTANT 
PRESSURE FLOW ALONG CONVEX WALL 
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FIG. 45 DISTRIBUTION OF v 2 /U| w AND w 2 /U?L FOR CONSTANT PRESSURE FLOW 
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ALONG CONVEX WALL 
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X (IN) 52.5 545 56.5 

STATION NO. 2 3 4 5 6 
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FIG. 52. VELOCITY PROFILES FOR SEPARATING FLOW ALONG CONVEX WALL 
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FIG. 53. SEMI -LOG PLOT OF VELOCITY PROFILES TO SHOW THE LAW OF THE WALL 
REGION FOR SEPARATING FLOW ALONG CONVEX WALL 






FIG. 55. DISTRIBUTION OF uw/U?, 

P 


FOR SEPARATING FLOW ALONG CONVEX WALL. 





FIG. 57 DISTRIBUTION OF u z /U| w FOR SEPARATING FLOW ALONG CONVEX WALL 
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FIG. 59. DISTRIBUTION OF wVU£ w FOR SEPARATING FLOW ALONG CONVEX WALL 
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Reference voltage for all traces 
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Z( in. 





VELOCITY 
TRACE AT 
CONSTANT y 




VELOCITY 
TRACE AT 
CONSTANT y 



Z 


X = 96.0 IN. 


FIG. 66. ASSUMED VORTEX STRUCTURE INSIDE BOUNDARY LAYER 
FOR CONSTANT PRESSURE FLOW ALONG CONCAVE WALL 
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FIG. 69. DISTRIBUTION OF u 2 /Up vv FOR CONSTANT PRESSURE FLOW ALONG CONCAVE 
WALL. MEASUREMENTS WERE TAKEN AT BOTH THE CREST AND TROUGH 
OF THE WAVE. 
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FIG. 71. DISTRIBUTION OF w 2 /U| w FOR CONSTANT PRESSURE FLOW ALONG CONCAVE WALL. 
MEASUREMENTS WERE MADE AT BOTH THE CREST AND TROUGH OF THE WAVE. 
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STATION X (IN) 
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FIG. 73 

SHEAR STRESS DISTRIBUTION FOR CONSTANT PRESSURE FLOW ALONG CONCAVE WALL. 
MEASUREMENTS WERE MADE AT BOTH THE CREST AND TROUGH OF THE WAVE. 
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FIG. 75. DISTRIBUTION OF uw/U AND vw/U^ w FOR CONSTANT PRESSURE FLOW 
ALONG CONCAVE WALL. MEASUREMENTS WERE MADE AT BOTH 
THE CREST AND TROUGH OF THE WAVE 
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DISTRIBUTION OF uw/Up W AND vw/Up W FOR CONSTANT PRESSURE FLOW ALONG CONCAVE 
WALL. MEASUREMENTS WERE MADE AT BOTH THE CREST AND TROUGH OF THE WAVE. FIG. 76 
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WALL. MEASUREMENTS WERE MADE AT BOTH THE CREST AND TROUGH OF THE WAVE. 
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PRESENT DATA 
PRESENT CALCULATIONS 

CALCULATIONS WITH HALF WALL CURVATURE 
CALCULATIONS WITH ZERO WALL CURVATURE 




CONVEX CURVATURE 
BEGINS AT X (IN) = 48 


) 8 

FIG. 79 

COMPARISON OF SKIN FRICTION, SHAPE FACTOR, DISPLACEMENT 
AND MOMENTUM THICKNESS CALCULATIONS AND PRESENT DATA 
FOR CONSTANT PRESSURE FLOW ALONG CONVEX WALL. 




COMPARISON OF VELOCITY PROFILE CALCULATIONS AND 
PRESENT DATA FOR CONSTANT PRESSURE FLOW AJ ONG 
CONVEX WALL 
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FIG 81. COMPARISION OF SHEAR STRESS CALCULATIONS AND PRESENT DATA FOR 
CONSTANT PRESSURE FLOW ALONG CONVEX WALL. 





DISTRIBUTION OF POTENTIAL VELOCITY AT THE WALL 
FOR SEPARATING FLOW ALONG CONVEX WALL. 

FIG. 82 





COMPARISON OF DISPLACEMENT AND MOMENTUM THICKNESS 
CALCULATIONS AND PRESENT DATA FOR SEPARATING 
FLOW ALONG CONVEX WALL. 
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COMPARISON OF SKIN FRICTION AND SHAPE FACTOR 
CALCULATIONS WITH PRESENT DATA FOR SEPARATING 
FLOW ALONG CONVEX WALL. 
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STATION 4 7 9 II 

X (I N )= 54.5 59 62.5 66 



COMPARISON OF VELOCITY PROFILE CALCULATIONS AND 
PRESENT DATA FOR SEPARATING FLOW ALONG CONVEX 
WALL. 
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SEPARATING FLOW ALONG CONVEX WALL. 




AND PRESENT DATA FOR CONSTANT PRESSURE FLOW 
ALONG CONCAVE WALL AT X ON) = 70 
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COMPARISON OF VELOCITY PROFILE CALCULATIONS 
AND PRESENT DATA FOR CONSTANT PRESSURE FLOW 
ALONG CONCAVE WALL AT X (IN) =96 
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FIG. 90 '8 

COMPARISON OF SHEAR STRESS CALCULATIONS WITH PRESENT DATA FOR CONSTANT 
PRESSURE FLOW ALONG CONCAVE WALL AT X(IN)=96. 



Figure 91* Adjustable outer wall test section. 
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r Photograph of the Equipment 

rtGUAS 9 3 for Pipe Flow Experiment 



DISTRIBUTION OF MEAN PRESSURE ALONG THE AXIS OF THE PIPE 

Fie use 94- 
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Figure 95 
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Figure 97 
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STATIC) S 
DELS (IN) 


Y\Ul V 

CONSTANT PkKS5UI<2 PLOW 
X (IN) - 2 4.00 

0 . J 8 3 THETA (IN) = 0 

y (r:i) 

0.00 4 
0 . 0 9 6 
0. OPH 
0,010 
0.012 
0.014 
0.010 
0.024 
0.029 
O.C 04 
0.039 
0.044 
0.054 
0.004 
0.074 
0.034 
0.094 
0. 104 
0.124 
0.144 
0.164 
0.184 
0. 204 
0.224 
0. 244 
0.264 
0. 284 
0.304 
0. 324 
0.344 
0. 364 
0.384 
0 , 404 
0.454 
0.504 
0.564 
0.604 
0.0 54 
0.704 
0.7 54 
0.804 
0. 854 
0.90 4 

0. 954 
1.00 4 
1.15 4 
1. 204 
1.304 

1. 404 
1.504 
1.60 4 
1.704 


ALONG CONVEX ALL 

!■’ (IN) - 0.0 

05 9 UPV ( FT/S EC) = 7v) . 2 7 


II (!* T/S EC) 

2 3, 

. 8 1 

26 , 

.29 

30 , 

. 1 3 

33, 

. 6 1 

35, 

■ 1 9 

36, 

> 9 3 

39, 

. 3 1 

40, 

, 9 2 

42 , 

, 14 

43, 

, 1 5 

43. 

, 9 2 

44 , 

,72 

45. 

,81 

47. 

, 1 6 

48. 

12 

4 9 , 

24 

50. 

1 3 

50. 

69 

52. 

, 4 1 

54. 

15 

55. 

96 

57. 

0 6 

58. 

10 

59. 

82 

61 . 

1 1 

62. 

10 

63. 

20 

64 , 

20 
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0 9 
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90 
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7 6 

67. 

9 8 

67. 

61 

68. 

77 

69. 

5 0 

69. 

97 
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12 
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27 

70. 

27 
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70. 
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70. 
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TAT i O M 2 
2LS(IN) - 


TABU' 

CONSTANT BkKSS'Jl' E TLOW 
/ (IN) - 52.5 3 

0.112 THETA (IN) = 0 


Y (IN) 

0 

. 0 0 4 

0, 

. 00 6 

0 

.008 

0, 

.010 

0, 

.0 12 

0. 

.014 

0 

.019 

0, 

. 024 

0, 

. 079 

0, 

.034 

0, 

.0 39 

0, 

.044 

0, 

.0 54 

0, 

. 064 

0, 

.074 

0 . 

. 084 

0, 

. 094 

0. 

, 104 

0. 

. 124 

0. 

, 144 

0. 

, 164 

0. 
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0. 

204 

0. 

224 

0. 
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0. 
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0. 

284 

0. 

304 

0. 
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0. 

344 

0. 

364 

0. 
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0. 

404 

0. 
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0. 

504 

0. 

554 

0. 

604 

0 . 

654 

0. 

7 04 

0. 

754 

0 . 

8 04 

0 . 

854 

0. 

904 

0. 

9 54 

1 . 

00 4 

1. 

10 4 

1 . 

204 

1. 

30 4 

1 . 

4 0 4 

1 . 

50 4 

1 . 

h » 4 

1. 

70 4 


ALONG C’UNVKY WALL 

I (IN) = 1 

085 UP W ( FT AS SC) = 

u (»• r/s !, c) 

26. 6 8 
2 9 . h 1 
32. 9 6 
35.87 

38.52 
40.3a 
43. 38 
4 5.39 
46.77 
48.19 
4 9.24 

50.26 
52. 1 7 
53.8 5 

55.22 

56.49 
57.7 3 
58.6 3 

60.13 
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62.23 
62.85 
6 3.45 
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64. «0 

65.26 
65.74 
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66.4 7 
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67.17 
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68. 3 4 
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70.46 
71.15 
71.70 

72.05 
72. 1 3 
72. 1 3 

7 1.93 
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71.40 
7 0.76 

70.14 
6 9.55 
68.95 
6 8.11 
6 7.6 7 
67.0 3 


; . 5 9 5 

7 9.7 4 
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■; i Luo A! 0:0 co»\ v T * k ali. 
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60.6 9 
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0.20 4 

62.47 
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6 3.03 
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6 3 . 6 8 
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64 . 2 9 
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0. 324 
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6 9 . 3 6 
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70.24 

0.604 
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0.654 

7 1.55 

0.704 
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0.754 

7 2.41 
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72.59 

0.854 
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0.90 4 
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0.954 
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1.00 4 
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STATION 6 
Dr:LS(n:j = 


’AiiLt. 


CONSTANT i’i KSSlil . fi FLO* 

ALONG 

CON vr X A ALL 


X ( I N) = 'j b . L> J 


:• ( i N ) - 

125 

TlilTA(TN) - 0 

. 09 1 

UP W (F':’/SEC) : 

Y (IN) 



u (h ;/si c*) 

0 . 0 <; 4 



22.98 

U . 006 



75. 46 

0.008 



28.6 0 

0.010 



31.94 

0.012 



34.64 

0.0 14 



7 6.89 

0. 019 



39.4 1 

0.024 



4 1 . 53 

0.020 



42.97 

0.044 



44.0 5 

0. 0 19 



4 5.19 

0.044 



4 5.92 

0.054 



47 .56 

0.064 



48.92 

0.074 



50.22 

0.0R4 



51 . 56 

0.094 



52.74 

0. 104 



53.81 

0. 124 



55.9 0 

0.144 



57.77 

0.164 



59.02 

0.184 



60. 45 

0.204 



61.41 

0. 224 



62.30 

0.244 



63. 1 1 

0.264 



63.74 

0.284 



64. 39 

0. 304 



6 4.84 

0. 324 



65 . 3 9 

0. 344 



65.92 

0. 364 



66.3 3 

0. 384 



66 . 86 

0.404 



67.21 

0.454 



68.19 

0.504 



69.11 

0. 554 



6 9.80 

0. 604 



70.47 

0. 65 4 



71.37 

0.704 



71.91 

0.754 



72.24 

0.804 



72.46 

0. 8 54 



72.55 

0.904 



72.44 

0.9 54 



72. 19 

1.004 



71.93 

1 . 104 



71 . 34 

1.204 



70.74 

1. .30 4 



70.13 

1.404 



69.50 

1.50 4 



68.99 

1.604 



68.27 

1.704 



67.67 


11 ..: C ■' 

7 5.7° 
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STATION 
DULS (IN) 


1 ABLE 


CONSTANT FhESSOF 

3 FLOrf 

ALONO CONVEX SALJ, 

x ( r n> = 

59.00 

K(IN) = 

. 1 30 TfiLTA (IN) 

= 0. 

09 3 1)8 W (FI /SEC) 

Y (IN) 


11 (FT/S PC) 

0.00 4 


21 . 84 

0.0Of> 


> 23.05 

0.008 


27.93 

0.01 0 


32.00 

0.012 


35. 1 5 

0.014 


36.8 3 

0.019 


39. UO 

0.0 24 


40. 3 9 

0.029 


4 1.82 

0.034 


42.74 

0.039 


43.31 

0.044 


44.21 

0.054 


45.68 

0.064 


47.27 

0.074 


48. 43 

0.004 


49.65 

0.094 


50.84 

0. 104 


51 .77 

0. 124 


53.99 

0. 144 


55.62 

0. 164 


57.61 

0.184 


58.70 

0.204 


60. 16 

0.224 


61.45 

0. 244 


62.21 

0.264 


62.97 

0. 284 


63.74 

0. 304 


64. 22 

0. 324 


64.83 

0.344 


65. 31 

0. 364 


65.88 

0. 384 


66 .25 

0.404 


66.77 

0.454 


67.67 

0.504 


68.63 

0.554 


69.48 

0.604 


70.28 

0.654 


71 . 00 

0.704 


71.63 

0.754 


72. 10 

0.804 


72. 3 7 

0.854 


72.41 

0.904 


72.23 

0.954 


71.99 

1.004 


71 . 72 

1 .104 


7 1.15 

1. 204 


70.5 3 

1.304 


69.93 

1.404 


69.35 

1. 504 


b 8 . 76 

1.604 


68.1b 

1.704 


67.57 


1 1 . r i.>0 

78.34 


?30 



STATION « 
DELS (IN) = 


TAHI.K 6 


CONSTANT I'RESSnrS P LOW 

A LG NT, 

CONVEX WALL 


X (IN) = fi-J.C-J 


b: (IN) = 

133 

THFTA(TN) - 0. 

09 4 

OP W (FT/SEC) = 

Y (IN) 



U (FT/SEC) 

0.004 



18 . n 

0.006 



20.0 5 

0.008 



23.50 

0.010 



27.51 

0.012 



30.4 8 

0.014 



32.82 

0.019 



36 . 27 

0.024 



38.28 

0. 029 



39.82 

0.034 



40.91 

0.039 



42.29 

0.044 



42.93 

0.054 



44. 47 

0.06 4 



45.66 

0.074 



46.85 

0.084 



48.02 

0.094 



49.32 

0. 104 



50. 16 

0. 124 



51.86 

0. 144 



54.33 

0. 164 



56.13 

0. 104 



57. 64 

0.204 



59.15 

0.224 



60.40 

0. 244 



62.17 

0. 264 



63.09 

0.284 



64.23 

0. 304 



65. 03 

0. 324 



65. 53 

0. 344 



66.09 

0.364 



66.57 

0. 384 



67.2 0 

0.404 



67 .64 

0.454 



68.50 

0.504 



69.60 

0.554 



70.19 

0.604 



71.08 

0.654 



71 .83 

0.704 



72.41 

0.754 



72.77 

0.804 



72.86 

0.8 54 



72.78 

0.904 



72.61 

0.954 



72.3 6 

1.004 



72.08 

1. 104 



71.52 

1.204 



70.97 

1. 304 



70.41 

1.404 



69.8 3 

1. 504 



69.28 

1.604 



68.69 

1.704 



6 8.11 
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STATTON 9 
DELS (IN) = 


TABLE 7 


CONSTANT TRESS (THE FLOW 

ALONG 

CONVEX WALL 


X (IN) - 67.00 


iMIN) = 

. 146 

THETA (TH) = 0. 

100 

UP w ( FT/S EC) 

y (in) 



U ( KT/S EC) 

0.00 4 



25. 35 

0.006 



25. 80 

0.008 



27.80 

0.0 10 



30.17 

0.0 12 



3 2 . 1 3 

0.014 



33.54 

0.019 



35.86 

0.024 



37.39 

0.029 



38.79 

0.0.34 



39.71 

0. 039 



40.46 

0.044 



41.49 

0.054 



42. 53 

0.064 



43.80 

0.074 



4 4.53 

0.084 



45.58 

0.094 



46.66 

0. 104 



47.80 

0. 124 



49.44 

0. 144 



51.27 

0. 164 



52.97 

0. 184 



54.53 

0.204 



56.42 

0.224 



57.86 

0. 244 



59.36 

0.264 



60.73 

0.284 



61 . 93 

0. 304 



63.23 

0.324 



64.12 

0. 344 



65.02 

0.364 



66.00 

0.384 



66.45 

0 . 404 



67. 10 

0.454 



68.47 

0.504 



69.63 

0.554 



70.48 

0.604 



71.29 

0.654 



71.85 

0.704 



72.36 

0.754 



72.8 3 

0.004 



73.03 

0.854 



73. 10 

0. 904 



72.98 

0.954 



72.77 

1.004 



72.52 

1 . 104 



7 2 . 0 2 

1.204 



71. 53 

1.304 



71.04 

1.404 



70.55 

1.504 



70. 0 5 

1.604 



69.53 

1.704 



68.9 9 


12. R2 0 
7n . 7 9 


232 



'C3 M 


r A I 1 1 O N 11 
KLS (IN) - 


CONSTANT I'N Ei>S(Jft IS 
X (IN) = 

0 . 155 TSIKTA (IN) 
Y(TN) 

0.004 
0 . 00 6 
0.008 
0.010 
0.012 
0.014 
0.019 
0.024 
0.029 
0.034 
0.039 
0.044 
0.054 
0.064 
0.074 
0.08 4 
0.C94 
0. 104 
0.124 
0.144 
0. 164 
0.104 
0.204 
0.224 
0. 244 
0.264 
0.284 
0. 30 4 
0.324 
0. 344 
0. 364 
0. 384 
0.404 
0.454 
0.504 
0.554 
0.604 
0.654 
0.704 
0.754 
0.804 
0.854 

0. 904 
0 . Q 5 4 

1.004 

1. 104 
1.204 
1.304 
1.404 
1 . 504 
1.604 
1. 704 


8 

ALONG CONVEX WALL 

H (IN) = 
0.105 II P W (FT/SEC) = 
l! ( FT/S KC) 
22.8 1 
2 4.40 
2b. 50 
30. 32 

32.23 
33. 30 

35.3 3 

37.01 
38.07 

39.23 

40 . 04 
40. oO 

42.04 
42.98 
43.9 3 
44.91 

45.84 
46 .89 
48.56 
50.13 

51.84 

53.01 

54.73 
56. 17 

57.4 5 
59. 18 
60.11 
61.39 
6 2.66 
63.61 
64.66 

65.6 0 
66 . 33 
68.00 
69.31 
7 0.23 

71.02 
71.86 
72.47 
72.82 
73. 13 
7 3.20 
73. 12 
72. 95 

72.74 

72.22 
71.72 

71.22 

70.7 3 
70.21 
6 9.69 
69.15 


T A 8 1. K 
FLOW 
7 1 . 00 


12.82 0 
7 8 . Q 8 
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TABLE 


9 


STATION 1 
Y (IN) 


CONSTANT PRESSURE FLON ALONG CONVEX WALI. 

X(IN) = 24.00 R (IN) 

0**2 V**2 H**2 UV UW 


0.70 0.176 
0.60 0.256 
0.50 0.893 
0.40 4.978 
0.35 8.278 
0.30 11.396 
0.25 14.743 
0.20 17.117 
0.15 19.664 
0.13 20.520 
0.11 21.443 
0.09 21.917 
0.07 23.005 
0.06 23.466 
0.06 23.46b 
0.05 23.926 
0.04 24.293 
0.03 25.527 
0.02 27.932 


0.053 0.026 

0.145 0.080 

0.913 0.526 

2.388 2.755 

3.113 5.065 

4.586 6.813 

6.243 0.568 

8.408 10.023 

9.890 11.007 

9.613 11.447 

10.156 11.979 

9.651 11.620 

8.512 12.649 

10.602 12.320 

12.320 
12.488 
12.752 
13.105 
13.151 


-0.001 0.004 

-0.024 -0.006 

-0.379 -0.005 

-1.902 0. 108 

-3.026 0.155 

-4.183 0. 247 

-5.492 0.328 

-6.626 0.107 

-8.035 0. 235 

-7.923 0.195 

-7.842 0.447 

-8.057 0.452 

-7.781 0.218 

-8.210 0.384 

0. 384 
0. 305 
0. 51 1 
0.023 
0.051 


0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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TABLE 10 


STATION 7 

Y (IN) 

1.50 

1.30 

1.10 

0.90 
0.70 
0.60 
0 . 50 
0 . 40 
0 . 36 
0.30 
0.25 
0.20 
0.15 
0 . 13 
0.11 
0.09 
0 . 07 
0.06 
0. 06 
0.05 
0 . 04 
0. 03 
0.02 


CONSTANT I'hESSUKE FLOW ALONG CONVEX WALL 


X(TN) = 

U *<:2 V **2 

0.138 0.018 

0.151 0.020 

0 . 1 H 4 0.127 

0.326 0.665 

1.867 2.079 

3.575 2.281 

4.309 1.912 

4.988 1.913 

5.495 2.015 

6.051 2.486 

7.309 4.916 

9.795 4.293 

11.967 7.664 

1 2 . 967 7 . 435 

14.001 8.662 

14.532 9.469 

14.977 9.015 

15.200 8.120 

15.200 

15.513 

15.479 

16.089 

16.533 


59.0 0 

W**2 UV 

0.065 - 0.002 

0.063 - 0.004 

0.098 - 0 . O 06 

0.294 - 0.012 

1.213 - 0.210 

2.175 - 0.395 

3.211 - 0.490 

3.855 - 0.538 

4.106 - 0.721 

4.625 - 0.935 

4.998 - 1.444 

6.244 - 2.458 

7.389 - 3.781 

8.250 - 4.166 

8.336 - 4.519 

8.534 - 4.884 

8.759 - 4.999 

8.784 - 5.053 

8.784 
8.642 
8.614 
9.021 
8.653 


K ( TN) = 11.520 


uw 


0.003 

0.0 

- 0.007 

0.0 

- 0.01 3 

o.c 

- 0.014 

0.0 

- 0.005 

0.0 

0.072 

0.0 

0 . 189 

0.0 

0.224 

0.0 

C . 315 

0.0 

0. 329 

0.0 

0.470 

0.0 

0.670 

0.0 

0 . 593 

0.0 

0.673 

0.0 

0.437 

0.0 

0 . 520 

0.0 

0.322 

0.0 

0 . 34 5 

0.0 

0 . 345 


0. 295 


O . 371 


0 . 325 


0 . 252 
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TABLE 11 


CONSTANT PRESSURE FLOW ALONG CONVEX WALL 


STATION 9 


Y (IN) 

U**2 

1 . 50 

0.428 

1 . 30 

0.3 80 

1.10 

0.400 

0.90 

0. 536 

0. 70 

1 . 994 

0. 60 

2.973 

0.50 

3.512 

0. 40 

4.013 

0.3 5 

4-. 737 

0. 30 

6.081 

0.25 

8.646 

0.20 

10.827 

0. 15 

13.563 

0. 13 

14.514 

0. 1 1 

15.564 

0.09 

16.109 

0.07 

16.506 

0.06 

16.520 

0.06 

16.520 

0.05 

16.868 

0.04 

17.063 

0.03 

1 7 .696 

0. 02 

18.699 


X (IN) 

= 67.00 

V**2 

W**2 

0.016 

0.364 

0.043 

0.307 

0.200 

0.3 33 

0.888 

0.415 

1. 708 

1 .254 

1.816 

1 . GOB 

1 . 952 

2. 124 

3.196 

2.631 

5. 047 

3.153 

6.516 

4.037 

9.081 

6.024 

8.723 

7.505 

10.292 

9.388 

9.408 

9.955 

10.424 

10.240 

9 .404 

10.203 

8.441 

10.179 

7.357 

10. 442 


10.442 


10.417 


10.445 


10.379 


10. 167 



V (IN) 

UV 

uw 

0.059 

-0. 136 

0.06 4 

0.001 

0.05 1 

-0.016 

-0.002 

0.027 

0.027 

0 .095 

0.200 

0.209 

0. 340 

0. 24 5 

-0.230 

0. 205 

-0.867 

0. 274 

-1.046 

0.131 

-3.242 

-0. 122 

-4. 197 

-0. 125 

-4.526 

-0. 244 

-4.641 

-0. 180 

-4.857 

-0.212 

-4 . 697 

-0. 122 

-4.950 

0. 105 

-5.843 

0- 094 


0.094 


0.007 


0.009 


0.011 


0.063 


12 


0 . 
0 . 
0. 
0. 
0 . 
0* 
0.0 
0.0 
0.0 
0.0 
0.0 
0. o 
0.0 
0. 0 
0.0 
0.0 
0.0 
0.0 


820 

vw 
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TABLE 12 


CONSTANT PRESSURE FLOW ALONG CONVFX WALL 


STATION 11 


X<1N) = 

71 .00 


i'(IN) - 

12. 

Y(IN) 

0**2 

V**2 

W**2 

UV 

UW 


1 . 50 

C • ft 76 

1.911 

1.121 

0.426 

-0 .015 

0.0 

1 .30 

0.733 

1.620 

0.857 

0.404 

0.046 

0.0 

1.10 

0.710 

1 . 361 

0.681 

0.401 

0.072 

0.0 

0.90 

0.820 

1.117 

0.761 

0. 182 

-0.023 

0.0 

0. 70 

2.C98 

1.539 

1.424 

0.115 

0.052 

0.0 

0.60 

2.632 

1.775 

1.746 

0.242 

0. 133 

0.0 

0. 50 

2.900 

1.920 

2.100 

0.209 

0. 134 

0. 0 

0.40 

3.207 

2.604 

2.363 

-0.313 

0. 055 

0.0 

0.35 

3.676 

2.757 

2.758 

-0.854 

-0. 102 

0.0 

0.30 

4.681 

5. 028 

3.097 

-1.941 

-0. 323 

0.0 

0.25 

6.497 

6. 950 

4.157 

-3.097 

-0.451 

0.0 

0.20 

9.927 

0.916 

5.606 

-4.328 

-0.446 

0.0 

0.15 

12.951 

9.478 

8.078 

-5.383 

-0.488 

0.0 

0.13 

13.903 

10.100 

8.505 

-5.478 

-0.699 

0.0 

0.11 

15.193 

9.762 

9.329 

-5.870 

-0. 382 

0.0 

0. 09 

15.515 

9.956 

9.565 

-5.817 

-0.500 

0.0 

0.07 

1 5. 915 

10.581 

9.098 

-5.741 

-0.298 

0.0 

0.06 

16.086 

8.444 

9.433 

-5.518 

-0.332 

0.0 

0.06 

16.086 


9.433 


-0. 332 


0.05 

16.174 


9.607 


-0. 200 


0,04 

16.337 


9.696 


-0.006 


0.03 

16.743 


9.263 


-0.215 


0.02 

17.635 


9.603 


0. 169 
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STATION 1 
DELS (IN) = 


TABLE 13 


SEPARATING FLOW ALONG CONVEX WALL 


0.095 
Y(IN) 
0.004 
0.006 
0.008 
0.010 
0.012 
0.0 14 
0.019 
0.024 
0.029 
0.0 34 
0.039 
0.044 
0.054 
0.064 
0.074 
0.084 
0.094 
0. 104 
0.124 
0. 144 
0. 164 
0. 184 
0.204 
0. 224 
0.244 
0. 264 
0.284 
0. 304 
0.324 
0. 344 
0.364 
0. 334 
0.404 
0. 454 
0.504 
0. 554 
0.604 
0.654 
0.704 
0.754 
0.804 
0.854 
0.904 

0. 954 

1.004 

1. 104 
1.204 
1. 304 
1.404 
1. 50 4 
1.604 
1.704 
1.804 
1.904 

2.004 


X (IN) = 24.00 

THETA (IN) = 0.068 


R(IN) = 
UPW (FT/SEC) = 
U (ET/SEC) 
21.91 
25. 64 
29.57 
32.47 
34.71 
36. 13 
38. 99 
40.39 

41.82 
42.93 
43.75 
44.73 

45.82 
47.12 
48. 15 

48.96 
49.69 
50.50 

51.82 
53.28 
54.43 

56. 11 

57. 31 

58. 96 
59.79 

60.95 
62. 18 
63.26 
64. 31 
64.98 

6 5. 68 
66. 36 

66.97 
68.14 
69.01 
69.63 
70.0 5 
70.4 9 
70. 76 
70.88 

70.95 

7 0.95 
70.95 
70.95 
70. 95 
70.95 
70.95 
70.95 
70.95 
70.95 
70.95 
70.95 
70.95 
70. 95 
70.9 5 


0 . 9 

70.9 5 
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TAULE 14 

SEPARATING PI.C-W ALONG CONVEX WALL 
STATION 2 X (IN) - 52. SO K (IN) = 

r > F. L S (IN) = 0.1 1« THETA (IN) = 0.na7 UPW(FT/SEC) = 


v (IN) 

U (FT/S FC) 

0.004 

25. 17 

0.006 

27. 52 

0.008 

30. 71 

0.010 

33. 30 

0.012 

35.3 3 

0.014 

36.91 

0.019 

40.72 

0.024 

42.21 

0.029 

42.92 

0.034 

44.42 

0.039 

45.13 

0.044 

45.86 

0.054 

4 7.33 

0.064 

49.24 

0.074 

50.42 

0.034 

51.62 

0.094 

53. 26 

0. 104 

5 4.31 

0.124 

56.3 0 

0. 1 44 

57.06 

0. 164 

59.07 

0.184 

60. 11 

0. 204 

60. 94 

0.224 

61 . 72 

0.244 

62.18 

0. 264 

62.79 

0.284 

63.25 

0. 304 

63.81 

0. 324 

64.22 

0. 344 

64.74 

0. 364 

65.25 

0. 384 

65.7 1 

0.404 

66.18 

0.4 54 

67.20 

0. 504 

68.00 

0.554 

68.99 

0.604 

69.73 

0.654 

70. 32 

0.704 

70.76 

0.754 

70.99 

0.804 

71.14 

0.9 54 

71.0 7 

0.904 

70.91 

0.9 54 

70.63 

1.004 

70.3 4 

1.104 

6 9.68 

1. 204 

6 9.02 

1. 304 

68. 37 

1. 404 

67.73 

1.50 4 

67.0 9 

1.604 

66.46 

1. 704 

65.8 4 

1.804 

65.22 

1.904 

64 . h 1 

2.004 

64.00 


10. c .°0 

77. J '1 
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T AbLK 15 

SEPARATING FLOW ALONG CONVEX WALL 
STATION 4 X (IN) = 54*50 K (IN) = 10.890 


DSLS (IN) = 0 . 153 TfILTA (IN) 

Y(IN) 

0.094 
0.005 
0.008 
0.010 
0.012 
0.014 
0.019 
0.024 
0.029 
0.034 
0.039 
0.044 
0.054 
0.064 
0.074 
0.084 
0.094 
0. 104 
0. 124 
0. 144 
0. 164 
0. 184 
0. 204 
0. 224 
0.244 
0.264 
0,284 
0. 304 
0. 324 
0.344 
0.364 
0. 384 
0.404 
0.454 
0.504 
0. 554 
0.604 
0.654 
0.704 
0.754 
0.804 
0.854 
0.904 

0. 954 

1.004 

1. 104 
1.204 
1. 30 4 
1 . 404 
1. 504 
1.604 
1. 704 
1.804 
1.904 

2.004 


0.105 UP W ( FT/S EC) - 74.61 

U (FT/S EC) 

20. 95 
22.33 
23. 37 
24 . 77 
26.98 

28.55 
31.82 

3 3.49 
34.93 
35 . 9 8 

37.09 

37.9 8 

39.13 

40.7 4 
41.91 

4 3.18 
44.29 

45.50 
47.57 

49.80 
51.44 

53. 07 

54.59 

55.73 
56 . 68 

57.48 
58. 17 
58.85 

59.49 
60.01 

60.51 

61.06 

61.55 

62.74 

63.75 

64.70 

65.59 
66.46 
67.12 

67.71 

68. 14 

68.3 3 
68. 39 
68.27 

68.04 
67.42 

66.80 
66.19 
65.53 
6 4.9 9 
64.39 
63.00 
63. 22 
62.64 

62.07 
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STATION r> 
DELS (i N) = 


TABLE 1<S 

SEPARATING FLO* ALONG CONVEX WALL 

X (IN) = r )6. r )0 r { I N ) - 11. TOO 

0.197 THETA (IN) = 0.127 UiH ( FT/S EC) = 72.32 


Y (III) 

U { FT/S OC) 

0.004 

18.60 

0.0 06 

20. 85 

0.000 

2 1 . 0 0 

o.oio 

2 3.-16 

0.0 12 

24.26 

0.014 

25.3 b 

0.014 

26.7 3 

0.0 24 

28. 38 

0.029 

29.36 

0.0 34 

29.69 

0.039 

30 .47 

0.044 

3 1.07 

0.054 

32.2 4 

0.064 

33.41 

0.074 

34. 82 

0.084 

35.6 6 

0.094 

36.59 

0. 104 

37.62 

0.124 

3 9.94 

0.144 

4 1.43 

0. 164 

43. 1 9 

0.184 

45.09 

0. 204 

4 6. 68 

0.224 

48.53 

0. 244 

49.69 

0.264 

51.07 

0. 284 

52.09 

0.304 

53.08 

0. 324 

54.14 

0. 344 

64.7 3 

0. 364 

55.57 

0.3 84 

66. 09 

0. 404 

56.7 0 

0.454 

58. 08 

0. 504 

59.26 

0. 554 

60. 34 

0.604 

61.28 

0.654 

62. 28 

0. 704 

63. 17 

0.754 

63.98 

0.804 

64.7 3 

0.854 

65. 36 

0.904 

65.67 

0.954 

65.85 

1.004 

65. 91 

1.104 

65.53 

1. 204 

64.94 

1.304 

64. 37 

1 . 404 

63.79 

1.504 

63.23 

1. 604 

62.67 

1 . 704 

62. 1 1 

1.804 

61.56 

1 . 904 

61.01 

2.004 

60.47 


?41 



TABLE 17 

SEPARATING FLOW ALONG CONVFX WALL 
STATION 7 X (I N) = 59.00 P(IN) = 11.520 

DELS (IN) = 0.270 THETA (IN) = 0.162 UP W ( FT/S EC) = 68.76 


Y(IN) 

o (t r/sr:c) 

0.004 

1 6 . 0 1 

0.006 

17. 10 

0.008 

17.54 

0.010 

18.78 

0.012 

19.51 

0.014 

20.4 0 

0.019 

21.13 

0.024 

21.90 

0.029 

22. 42 

0.084 

22.9 3 

0.039 

23.69 

0.044 

23. 91 

0.054 

24.86 

0.064 

25.5 3 

0.074 

26.18 

0.084 

26.96 

0.094 

27.13 

0. 104 

27.97 

0.124 

29.62 

0.144 

31 .06 

0.164 

32.78 

0. 104 

33. 92 

0. 204 

36.05 

0. 224 

37.15 

0.244 

38.99 

0.264 

ao. 22 

0.284 

41.60 

0. 304 

43.02 

0.324 

a4. 14 

0. 344 

45.40 

0.364 

46.39 

0. 384 

47.48 

0. 404 

48.05 

0.454 

50. 1 3 

0.504 

51.47 

0.554 

53.09 

0.604 

54.31 

0.654 

55.28 

0.704 

56.49 

0.754 

57.35 

0.804 

58.29 

0.854 

59. 33 

0.904 

60.14 

0.954 

60.87 

1.004 

61 . 52 

1. 104 

61.89 

1.204 

61 . 90 

1. 304 

61.40 

1.404 

60.87 

1. 504 

6 0.34 

1.604 

59.82 

1.704 

59.3 0 

1.804 

58.79 

1.904 

58.28 

2.004 

57.78 
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T\HL 1-; 1 ‘i 

SEPAhATJNG FLOW ALONG rONVFX WALT. 

STATION 0 X (IK) = 6P.50 T; ( IN ) = 11. 

DLLS (IN) = O.J°1 THf'TA (IN) = 0 . 2 1 2 UPHfFT/SLC) = f>5 


Y (IN) 

11 ( FT/S EC) 

0.004 

15.33 

0.006 

15.28 

0.008 

15.71 

0.010 

16.13 

0.012 

16.61 

0.014 

16. 66 

0.019 

17. 96 

0.C24 

18.81 

0.029 

18.97 

0.034 

19.46 

0.039 

19.36 

0.044 

19.71 

0. 054 

20. 1 5 

0.064 

20.28 

0.074 

20.82 

0.084 

20.95 

0.094 

21.42 

0. 104 

21.71 

0. 124 

22.56 

0. 144 

23. 38 

0. 164 

24. 32 

0. 194 

25.53 

0.204 

26.16 

0.224 

27.51 

0. 244 

28.54 

0.264 

29.24 

0. 284 

30.69 

0. 304 

32.01 

0. 324 

33.24 

0.344 

3 4.31 

0. 364 

35.60 

0. 384 

36.74 

0.404 

37.52 

0.454 

40.2 3 

0.504 

42.72 

0.554 

44. 52 

0.604 

46. 16 

0.654 

47. 30 

0.704 

48.54 

0.754 

49.58 

0.804 

50.53 

0.854 

51.31 

0.904 

52.20 

0.9 54 

52.99 

1.004 

53.65 

1. 104 

55.5 ) 

1. 204 

57.07 

1.304 

58.00 

1. 404 

57.99 

1 . 504 

57.76 

1.604 

57.27 

1.704 

56.79 

1.804 

56.32 

1.904 

55.84 

2. 004 

55. 37 


<iao 

.so 
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TAt'.LE 19 

SEPARATIN',; FLOW ALONfl CONVEX WALL 
STATION 11 X (IN) = b 6 . 0 0 I- (IN) = 12. 

DELS (IN) = 0.485 TiiETA (IN) = 0.23? UPW (FT/SEC) = t>2 


Y(IN) 

U (FT/SKC) 

0.004 

11.39 

0 . 0 0 6 

11.34 

0.008 

11.28 

0.0 10 

11.23 

0.012 

11.17 

0.014 

11.12 

0.019 

10.98 

0.024 

11.93 

0.029 

12.01 

0.034 

12.59 

0.039 

13. 34 

0.044 

13. 86 

0.054 

14.00 

0.064 

14.65 

0.074 

14.44 

0.084 

14.66 

0.094 

15.04 

0. 104 

15.25 

0. 124 

15.85 

0, 144 

16.37 

0. 164 

17.48 

0. 184 

18.40 

0.204 

19.52 

0.224 

20.^2 

0.^44 

21.20 

0. 264 

22.30 

0.284 

23.24 

0. 304 

24.35 

0.324 

25.26 

0. 344 

26. 38 

0.364 

27 . 78 

0. 384 

28. 97 

0.404 

29.83 

0.454 

32.44 

0. 504 

35. 39 

0.554 

37.82 

0.604 

39.89 

0.654 

41 .89 

0.704 

43. 52 

0.754 

45.00 

0.804 

45,88 

0.854 

47.03 

0.904 

47.93 

0.954 

48. 88 

1.004 

49.50 

1. 104 

51.1 1 

1. 204 

52.53 

1. 304 

54.01 

1.404 

55. 08 

1.504 

55.6 9 

1.604 

55. 54 

1.704 

55.31 

1.804 

54. b8 

1.904 

54.26 

2.004 

53.84 
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STATION 12 
DSLS (IN) = 


TAULR 20 

SEPARATING FLOW ALONG CONV 
X ( I N) = b 9 . 5 0 

0.608 THETA (IN) = 0.270 

Y (IN) 

0 . 004 
0.006 
0.008 
0.010 
0.012 
0.014 
0. 019 
0.024 
0.029 
0.034 
0. 0 39 
0.044 
0.054 
0.064 
0.074 
0.004 
0-094 
0. 104 
0. 124 
0. 144 
0.164 
0. 184 
0. 204 
0.224 
0. 244 
0.264 
0. 284 
0.304 
0.324 
0.344 
0. 364 
0.384 
0. 404 
0.454 
0.504 
0.554 
0. 604 
0.654 
0.704 
0.754 
0.804 
0.854 
0.904 

0. 954 
1.004 
1 . 104 

1. 204 
1.304 
1. 404 
1.504 
1.604 
1.704 
1.804 
1 . 904 
2.00 4 


X WALT. 

K ( I N) = 
UP W (FT/SEC) = 
U (FT/S EC) 

13.04 
1 2 . 9 9 

12.95 

12.90 
12. 83 

12.91 

12.70 

12.58 
12 . 66 
13.03 
13.10 

13.18 

13.05 

1 3 . 9 3 
13.72 

13.51 
13.75 

14 . 07 
14.53 
14.48 

14.51 

14.96 
15. 31 

15.58 
16.24 
16.86 

17. 19 

17.71 

18.63 

19.05 

19.91 

20.85 
21. 34 
23.38 

25.64 
27.95 
30. 32 
33. 17 

35.58 
37.69 
39.40 

41.07 

42.59 
43.90 
4 5.00 

46.7 7 

4 8.25 
49.88 
50 . >3 5 

51.9 3 

52.85 

5 3. 29 
53.22 

52.92 

52.51 


12. 8 ? 0 
61.39 
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TABLE 21 


STATION 1 

V (IN) 

3.00 

1.00 

0 . 00 
0.80 
0.7 0 
0. 60 
0. 55 
0. 50 
0.45 
0.40 
0.3 5 
0.30 
0.25 
0 . 20 
0. 15 
0 . 10 
0.05 
0. 05 
0.04 
0.0 3 
0.02 


S h'P AP AT I NCi FLOW ALONG CONVEX NALI. 


X (IN) 

U**2 V**2 

0.146 0.057 

0.158 0.057 

0.173 0.080 

0.214 0.097 

0.408 0.286 

1.591 1.094 

3.294 1.643 

5.559 2.883 

7.743 3.361 

9.800 4.411 

13.524 5.400 

15.546 6.608 

17.27o 7.102 

18.470 8.608 

21.140 9.102 

23.306 9.192 

24.760 8.902 

24.760 
25.414 
28.137 
20.788 


24 . 00 


W**2 

UV 

0.000 

-0.010 

0.006 

-0.010 

0.020 

-0.031 

0.060 

-0.051 

0. 189 

-0. 102 

1 . 067 

-0.536 

1.820 

- 1 . 236 

2.696 

-2.052 

4.619 

-2.878 

6.4 56 

-3.750 

7.712 

-4.649 

9.443 

-5.559 

10. 091 

-6.516 

1 0.667 

-7.223 

11.097 

-8.116 

12.176 

-9.245 

12.206 

-6.135 

12.20b 


12.859 


13.503 


14.528 



(IN) 

it 

o 

o 

uw 

vw 

-0.001 

0.019 

-0.001 

0.019 

0. 0 

0.023 

0.003 

0.024 

0.006 

-0.0Q4 

0.020 

-0.116 

0. 034 

-0.277 

0.113 

-0. 190 

0.104 

-0.264 

0.045 

-0.222 

0.085 

-0.371 

-0. 063 

-0.409 

-0.091 

-0.476 

-0. 070 

-0.496 

-0. 046 

-0.326 

-0. 074 

-0.236 

-0. 139 

-0.191 

-0.139 


0. 029 


-0.030 


-0.081 
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TABLE 22 



s 

EPAK ATING 

FLOW ALONG 

CONVEX 

WALL 


STATION 4 


X (IN) 

= 54.50 


P(IN) = 

1 0 . 89 0 

Y (IN) 

U**2 

V**2 

W**2 

ov 

UW 

vw 

1 . 40 

0. 159 

-0 . 384 

0 . 0 b J 

0.249 

0.007 

0.019 

1.20 

0.250 

-0. 212 

0.193 

0.193 

0.029 

0. 045 

1 . 0 0 

0 . 9l 1 

0.248 

0. 701 

'.) . . 9)5 

0 . 16 9 

0.44 2 

0.00 

2.026 

1.159 

1 .609 

0.318 

0. 671 

0.391 

0.00 

3.011 

1 . 574 

1 .859 

0.413 

0.505 

1.194 

0.70 

4.019 

2. 180 

2.656 

0.355 

0.772 

1.000 

0.60 

4 . 954 

2.793 

3.254 

0.161 

0.58 3 

1.393 

0.50 

5.460 

3.295 

3.651 

-0.234 

0.532 

1.132 

0.4 0 

5.899 

2.74 1 

3.943 

-0.035 

0.48 0 

1.082 

0.30 

6.132 

2.841 

4.067 

-0.107 

0. 446 

1.035 

0.20 

6.945 

2.558 

4.754 

-0.544 

C.239 

0.907 

0. 15 

8.450 

3.232 

5.750 

-1.307 

-0. 127 

1 .420 

0.10 

12.718 

2. 622 

7.542 

-2.113 

-0. 335 

0.948 

0.05 

14.740 

5. 1 36 

9.608 

-3.933 

-0.778 

1.205 

0.05 

14.740 


9.608 


-0.778 


0.02 

12.602 


5.734 


-0. 265 
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TABLE 23 


STATION 9 

y (in) 

2.00 
1.80 
1.60 
1. 40 
1 . 20 
1.00 
0. 90 
0.80 
0.70 
0.60 
0.50 
0.40 
0. 30 
0.20 
0 . 10 
0.05 
0.05 
0.02 


SEPARATING FLOW ALONG 



X(IM) - 

62.50 

U**2 

V**2 

W**2 

0. 350 

- 0 . 347 

0. 188 

0.358 

-0. 172 

0. 243 

0.361 

-0. 283 

0.156 

0.532 

0.065 

0.250 

1.22 1 

0. 394 

0.813 

2.402 

1 . 101 

1.643 

2.630 

0.797 

1 . 955 

3.250 

0.713 

2.152 

3.653 

0.992 

2.372 

4. 351 

1. 213 

2.863 

6.795 

1. 262 

4.220 

11.144 

2. 660 

7.169 

14.963 

3.885 

9.820 

15.452 

4.560 

9.944 

12.169 

4. 325 

7.226 

9.967 

3.011 

6.137 

9.967 


6.137 

9. 166 


5.105 


CONVEX WALL 



R(IN) - 

1 1 . 84 

UV 

u w 

vw 

0.414 

0.00 6 

C. 1 88 

0.427 

0 . 007 

0.172 

0.352 

0.009 

0,163 

0.401 

C. 019 

0.283 

0.146 

-0. 139 

0.549 

-0.086 

-0.275 

0.456 

-0.082 

-0 • 4 4 4 

0.290 

-0.073 

-0. 548 

0.196 

-0.273 

-0. 608 

0.0 86 

-0.442 

-0.616 

0.0 05 

-0.709 

-0. t>2 5 

0.1 89 

-2. 142 

-0. 483 

0.7 06 

-3.366 

-0. 408 

1.100 

-4.672 

-0.257 

1.781 

-3. 187 

-0. 186 

0.549 

-2.295 

-0.060 
-0.060 
-0. 102 

0.452 
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TABLE 24 


SEPARATING FLOW ALONG CONVEX WALL 


STATION 11 


X (IN) - 

66 . 00 


MIN) = 

12. 820 

y (in) 

U**2 

V**2 

W* + 2 

UV 

UK 

VW 

2.00 

1 . 207 

-0. 302 

0.565 

0.625 

-0.025 

0.501 

1.80 

1.137 

-0.272 

0.454 

0.87 8 

-0.029 

0.3 83 

1.60 

1 .044 

-0.541 

0.331 

1.012 

0.003 

0.227 

1. 40 

1 .463 

-0.557 

0.634 

1.399 

0.00 2 

0.270 

1.20 

2.569 

0.909 

1 .388 

0.651 

-0. 286 

0.4 82 

1.00 

3.317 

1 . 097 

1 .934 

0.191 

-0.618 

0. 341 

0.90 

3.526 

1.203 

2.165 

0.045 

-0.708 

0.471 

0. 80 

3.942 

0. 979 

2.570 

-0.079 

-0.745 

0.417 

0.70 

5.196 

1 . 787 

3.480 

-0.648 

-0.623 

0.631 

0. 60 

7.905 

1 . 828 

5.018 

-0. 902 

-0. 345 

0.9 67 

0.50 

1 1 .354 

2.793 

7.103 

-1.874 

-0.063 

1.196 

0.40 

14.703 

4.580 

9.557 

-3. 325 

0.022 

1.702 

0.30 

15.509 

5.090 

10.366 

-3.874 

-0.057 

1.664 

0. 20 

14.560 

4.263 

8.410 

-3. 356 

0.026 

1 .040 

0. 1 0 

10.682 

3. 460 

6.967 

-2.503 

0. 150 

0 .999 

0.05 

8.857 

2.561 

5.543 

-1.727 

0. 187 

0.262 

0.05 

8.857 


5.543 


0. 187 


0.02 

7.778 


4.280 


0. 165 
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TABLE 25 


SEPARATING PLOW ALONG CONVEX WALL 


STATION 12 


X (IN) 

= 69.50 


P(TN) = 

12. 82 

Y (IN) 

IJ**2 

V* *2 

W* *2 

UV 

UK 

VW 

2. 00 

4 . 1 1 j 

1. 2 3t> 

1 .892 

2.32 5 

-0. 147 

1 .839 

1.80 

3.60 3 

1 . 489 

1 .244 

1.993 

0.023 

1.428 

1.6 0 

3.10 3 

1 .450 

1.169 

1.530 

-0.082 

0.980 

1.40 

3 . 4 99 

0 . 664 

1.522 

1 . 88 1 

-0. 345 

0.774 

1.20 

3.907 

1 . 925 

2. 132 

0.557 

-C. 69 1 

0.7 32 

1 . 00 

4.495 

2.405 

2.991 

-0.299 

-0.685 

0.613 

0.90 

5.858 

2. 389 

4.045 

-0.859 

-0. 584 

0.74 4 

0.80 

8.502 

2. 153 

5.826 

-1 . 176 

-0.411 

0.560 

0.70 

1 1.985 

3. 142 

7. 322 

-2.208 

-r. 153 

0.885 

0. 60 

15.402 

2.710 

9.762 

-2.816 

-0.493 

2.177 

0. 50 

17.019 

5.309 

1 1. 581 

-4.279 

-C. 507 

2.0 8C 

0.40 

17.899 

3.601 

1 1.068 

-3. 190 

-0.921 

2.218 

0. 30 

16.634 

2.981 

11 .354 

-3.215 

-1. 254 

2.786 

0. 20 

1 4.534 

3.286 

8.707 

-2.726 

-0. 925 

2.603 

0.10 

10.651 

4.211 

6. 083 

-2.973 

-0 .623 

1.311 

0. 10 

10.651 


6.083 


-0. 623 


0.05 

9.295 


5.005 


-0.514 


0.02 

5.511 


2.411 


-0.287 
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STATION 1 
DELS (IN) = 


TADLE 2‘) 


CONSTANT PPLSSEJK 

rl FLOW ALONC 

COLCAV5 ‘’ALL 


x (IN) = 

2 4.0 :) 

P (IN) = 

087 

TUFT A (IN) 

= 0.062 

npw (FT/SEC) 

Y (IN) 



U (FT/SKC) 

0.006 



25.50 

0.0 08 



29.0 9 

0.0 10 



32.77 

0.015 



36.9 6 

0.020 



39. 15 

0.025 



40. 64 

0.0 30 



41 .84 

0.015 



42.3 3 

0.040 



4 3.93 

0.050 



45. 33 

0.060 



4 6.52 

0.070 



47.82 

0.080 



48.91 

0.090 



50.01 

0. 100 



50.80 

0. 120 



52.20 

0. 140 



53.89 

0. 160 



55.39 

0. 180 



57.18 

0. 200 



58 . 10 

0. 220 



59.77 

0. 240 



61 . 06 

0. 260 



62.0b 

0. 280 



63. 46 

0. 300 



64.05 

0. 320 



65. C5 

0. 340 



65.75 

0. 360 



6b. 64 

0. 380 



67.14 

0.400 



67. 64 

0.450 



68.54 

0.50 0 



69 . 23 

0. 550 



6 9.73 

0. 600 



70.13 

0. 650 



70.2 3 

0.7 00 



70. 3 3 

0.750 



70.53 

0.800 



70.6 3 

0.850 



70.63 


0. n 

70.*) 
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STATION 
DSLS (IN) 


TALUK 2/ 

CONSTANT n.FSSHKr’ FLOW ALOtM CONCAVt * ALL 

X (IN) = 70.00 L (IN) = -17. 0 30 

0.209 TH F TA ( IN ) = 0.2.3‘i Ui'W (FT/SEC) = 0 1.20 


Y (IN) 

U (J-T/SPC) 

0.0 08 

3 3.19 

0.010 

34 . 53 

0.015 

3 6.68 

0. 020 

38.11 

0.025 

39.0 3 

0.030 

39.85 

0.035 

40.47 

0.040 

4 0.98 

0.060 

4 2.41 

C.080 

43. 20 

0. 100 

4 3.9 5 

0. 130 

44.68 

0. 160 

44.40 

0. 190 

45.49 

0. 220 

46. 00 

0. 250 

46. 51 

0. 280 

46.82 

0. 310 

47.23 

0. 340 

47.54 

0, 370 

4 7.9 5 

0. 400 

48. 36 

0.450 

48.77 

0. 500 

49.48 

0. 550 

50.10 

0. 600 

50. 40 

C. 650 

5 1 . 0 2 

0.700 

51.43 

0.750 

52.0 4 

0. 000 

52.45 

0.050 

52.86 

0.000 

53.40 

J. 950 

54 . 0 9 

1.000 

54. 5 0 

1.100 

55.94 

1. 200 

57. 37 

1 . 300 

'.1 8.91 

1.400 

6 0.3 4 

1 . 500 

6 1.8 3 

1.600 

6 3.31 

1.700 

64. 54 

1 . 800 

65.* 7 

1.930 

6 6.69 

2.000 

67 . 

2.500 

70. 09 

3.000 

72. 33 
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STATION 
DHLS (IN) 


TABLE 26 

CONSTANT PH ESS OR E FLOW ALONG CONCAVE WALL 

X ( IN) = 70.00 t'(IN) = -17.830 

0.200 THETA (IN) = 0.153 UPW ( FT/S EC) = 61.20 


Y (IN) 

U (Y T/S EC) 

o.ooo 

31.18 

0.C10 

32.70 

0.015 

35. 52 

0 . 0 20 

37.44 

0.025 

* 8 . 4 5 

0.030 

39.2 6 

0.035 

39.76 

0.040 

40. 37 

0.060 

41.70 

C.OOO 

42.89 

0.100 

43.60 

0. 130 

44.41 

0. 160 

44.91 

0. 100 

4 5.52 

0. 220 

46.02 

0.250 

4 6.42 

0.280 

47.03 

0.310 

4 7.43 

0.340 

47.84 

0. 370 

48. 24 

0. 400 

48.85 

0.450 

49.55 

0. 500 

50.66 

0.550 

51.64 

0. 600 

52.68 

0.650 

53.99 

0.700 

55.10 

0.750 

56.-21 

0.800 

57.43 

0.850 

58.43 

0. 900 

59.24 

0.9 50 

60.45 

1.000 

61.46 

1.100 

63.28 

1. 200 

6 4.59 

1.300 

65.50 

1. 400 

66.00 

1.500 

66.31 

1 . 600 

66.71 

1.700 

67.01 

1.800 

67 . 32 

1 .900 

67.7 2 

2. 000 

68. 32 

2.500 

70.5 5 

3.000 

72. 97 
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STATION <4 
DELS (IN) = 


TABLE 2<> 


CONSTANT I'KL’SSBKE FLOW 

ALONG 

CONCAVE WALL 


X ( I N ) = 96.00 


8 (IN) =• 

.250 

THETA (IN) = 0. 

227 

UPW (FT/SEC) 

Y (IN) 



0 (KT/SFC) 

o.oo r> 



37. 69 

o.roo 



39.4 0 

0.010 



40.74 

0.015 



43.17 

0.020 



4 4.43 

0.025 



45 . 30 

0. 0 30 



4 6.22 

0.035 



47.0 6 

0.040 



47 . 5 9 

0.060 



49.69 

0.080 



51.17 

0. 100 



52. 22 

0. 150 



54.33 

0.200 



55.27 

0. 250 



55. 40 

0. 300 



56. 54 

0. 350 



56. P 5 

0.400 



57.38 

0.500 



57. 91 

0. 600 



58.42 

0.700 



58.75 

0.0 00 



59.17 

0.900 



59. PO 

1.000 



60. 1 1 

1.200 



61 . 17 

1.400 



62.22 

1.600 



63.17 

1 . 800 



6 3.80 

2.000 



64.33 

2.200 



64.85 

2. 400 



64. 85 

2.600 



65. 17 

2. 800 



65. 5 9 

3.000 



6 6.43 

3 . 20 0 



67.17 

3.400 



6 0 . n 4 

J.600 



70.12 

3.800 



71.70 

4.00 0 



73.3 8 

4.200 



74.0 5 

4.400 



75.70 

4.6 00 



76 . 75 

4.800 



77. 30 

5.00 0 



70. 9 1 


tt 7 () 
. h 0 
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STATION 5 
DELS (IN) = 


TAULE 

CONSTANT PL ESS OR E FLOW 
X (IN) = 96.00 

0.124 TiiETA ( I M) - 0 

Y(IN) 

0.006 
0.008 
0.010 
0.015 
0.0 20 
0.025 
0.0.10 
0.0 35 
0.040 
0.060 
0. 080 
0. 100 
0.150 
0 . 200 
0. 250 
0. 300 
0. 350 
0.400 
0. 500 
0. 600 
0.700 
0.800 
0.900 
1. 000 
1.100 
1. 200 
1. 300 
1. 400 
1.500 
1 . 600 
1 . 700 
1 . 800 
1 . 9o0 
2.000 
2.200 
2. 400 
2. 600 
2. 800 
3.000 


ALONG CONCAVE WALL 

K (IN) = -21.870 
113 UPW (FT/SEC) = 61.60 


U ( FT/S 

SC) 

31 

. 26 

35 

. 17 

37 

. 94 

41 

. 8 5 

43 

. 91 

4 5 

.24 

46 

. 39 

47 

.09 

47 

.81 

50 

.28 

51 

.93 

52 

.96 

54 

. 7 0 

55. 

.39 

56. 

. 35 

56. 

. 86 

57. 

. 17 

57. 

, 4 8 

57, 

. 99 

58. 

,40 

58, 

, 82 

59, 

,29 

59. 

. 74 

60, 

, 96 

60. 

, 9 8 

61 . 

8 0 

6 2. 

62 

6 3 . 

34 

64 ♦ 

16 

64. 

78 

65. 

50 

6 6. 

2 2 

66 . 

7 3 

67. 

2 5 

68. 

28 

69. 

1 0 

69. 

92 

70. 

5 4 

71. 

16 
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TA8LE 31 


CONSTANT PRESSURE FLOV ALONG CONCAVE WALL 


STATION 1 
Y (IN) 


. V 'J u . 

1.00 0 . 

0.90 0. 

0. 80 0. 

0.70 0. 

0.60 1. 

0.55 3. 

0.50 6. 

0.45 8. 

0.40 11. 

0.35 14. 

0.30 16. 

0.25 17. 

0.20 19. 

0.15 22. 

0.10 24. 

0.05 25. 

0.05 25. 

0.04 26. 

0.03 29. 

0.02 30. 


X (IN) = 

0**2 V**2 

129 0.G57 

134 0.057 

143 0.063 

160 0.119 

286 0.455 

840 1.188 

406 2.445 

139 3.348 

670 4.462 

278 5.318 

142 6.259 

071 7.110 

841 8.274 

697 8.919 

263 9.773 

178 9.132 

615 9.339 

615 

440 

048 

763 


24.00 



nv 

0.007 

-0.008 

0*007 

-0.008 

0*008 

-0.024 

0.034 

-0.039 

0.306 

-0.093 

1.274 

-0.716 

2.058 

- 1.46 1 

3.281 

-2.343 

5.340 

-3.306 

7. 102 

-4.159 

8.898 

-4.929 

9.920 

-5.743 

10.290 

-6.694 

11.368 

-7.410 

11.623 

-8. 073 

12.920 

-8.246 

12.900 

-8.341 

12.900 


13.135 


13.8 92 


14.816 



P (IN) = 0.0 

uw vw 

0.300 -0.001 

0.001 - 0.001 
0.001 - 0 . ('0 2 
0.002 -0.004 

0.014 -0.022 

0.024 -0.040 

0.049 -0.00) 

0.105 -0.175 

0.123 -0.205 

0.128 -0.383 

0.116 -0.432 

0.021 -0.395 

0.166 -0.276 

0.170 -0.469 

0.101 -0.432 

0.068 -0.502 

0.049 -0.582 

0.049 
0.103 
0. 115 
0.003 
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TABLE 32 


CONSTANT PRESSURE FLOW ALONG CONCAVE WALL 


STATION 2 


X (IN) 

= 70.03 


fl (IN) 

= -17. R 3 0 

7 (IN) 

U**2 

V**2 

W**2 

II V 

uw 

VW 

2. 5l) 

0.265 

0 . 27 4 

0.171 

-0.007 

o . o 1 2 

-0. 1 1 '3 

2.3 0 

0.428 

0.633 

0.428 

-0.123 

0.064 

-0.133 

2. 10 

1 . C 2 3 

1 . 256 

1.189 

-0.237 

C. 27 9 

-0.152 

1.90 

2.729 

2. 855 

3.109 

-0.979 

0.381 

-0.523 

1. BO 

4.5 07 

4.117 

4.300 

-1.818 

1. 607 

- 1 . 3UO 

1.70 

10.465 

7. 393 

9.803 

-4.241 

3. 80 5 

-2.063 

1 . 60 

12. 165 

4. 9 66 

7.945 

-3.117 

4. 036 

-1.802 

1.50 

15.225 

10.543 

1 1 .865 

-6.996 

5. 688 

-2.427 

1 . 40 

1 9. 158 

8.72 1 

13.893 

-6 . 440 

6.845 

-2.743 

1.30 

23. 977 

11.927 

16.053 

-8.810 

8. 312 

-3. 8.3 5 

1.20 

26.683 

16.379 

18.762 

-1 1 . 497 

9.615 

-4.943 

1.10 

28.707 

19.113 

20.287 

-12.79 1 

10. 330 

-5.818 

1 . 00 

31.184 

17.434 

23.252 

-11.414 

11. 335 

-5. 220 

0. 90 

31.951 

21.482 

24.280 

-14.246 

1 1. 762 

-6.212 

0. 80 

32. 1 99 

20.588 

25.043 

-13.577 

11.361 

-6 . 040 

0.70 

3 1. 592 

20. 573 

25.258 

-12.980 

10.069 

-4 . 070 

0.60 

3 1.202 

19.287 

23.820 

-12.538 

8. 897 

-3.188 

0.50 

30.065 

19.652 

23.688 

-13. 062 

7.678 

-4.034 

0.40 

28.760 

16.533 

21 . 046 

-11.204 

5. 587 

-3.269 

0.30 

27.672 

16. 946 

20.323 

-10.866 

3.714 

-1.1 48 

0. 20 

26 . 840 

14.287 

17.352 

-9. 170 

2. 803 

-1 .605 

0. 15 

26.070 

14.054 

17.045 

-8.561 

2.587 

-2. 100 

0. 10 

26.757 

11.789 

16.431 

-7.190 

2. 6 00 

-2.952 

0. 05 

28.105 

8.112 

16.755 

-5.260 

2. 497 

-2.923 

0.05 

28. 1 C 5 


16.755 


2. 497 


0.04 

28.609 


17.198 


2. 389 


0. 03 

20.967 


17.434 


1. 72 8 
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TABLE 33 


CONSTANT PRESS UK E FLOW ALONG CONCAVF WALL 


STATION 3 


X (IN) = 

7 0.00 


R(TN) = 

-17. 830 

y (IN) 

0**2 

V*<=2 

W**2 

uv 

II Vi 

V W 

1.60 

0.305 

0.354 

0.2 38 

-0.041 

c. 00 1 

0.050 

1.60 

0.466 

0.649 

0. 4 52 

-0.159 

0 .0 06 

0.031 

1 . uo 

0.929 

0.977 

0.934 

-0.169 

-0.008 

0.015 

1.30 

2.799 

2.646 

1 .704 

-1 . 579 

-0.092 

0.140 

1.20 

5.862 

5.391 

2.930 

-3.740 

-0. 376 

0.585 

1.10 

1 1.722 

7. 308 

6.091 

-6.212 

-0.406 

0.408 

1 . 00 

16.769 

8.896 

10.578 

-7 . 9/3 6 

-1.216 

0.459 

0.90 

23.398 

11.293 

12.566 

-10.232 

-1.182 

0.423 

0.80 

27.340 

15.746 

16.376 

-13.797 

-1. 621 

0 . 969 

0.70 

29.905 

15. 233 

18. 886 

-13.877 

-1. 548 

0.138 

0.60 

30.405 

15.230 

19.516 

-1 3.477 

-2.494 

1.651 

0 . 50 

29.866 

18. 552 

20.790 

-14.716 

-2. 106 

0 . 205 

0.40 

20.639 

17.947 

19.940 

-13.669 

-2. 637 

1 . 0b9 

0. 30 

25.624 

17.314 

20.202 

-11.628 

-2. 177 

0 . 0 85 

0.20 

25.337 

15.047 

18.767 

-9.610 

-3. 352 

1 .679 

0. 15 

24.552 

13.944 

19.534 

-6.951 

-3. 425 

0.597 

0 . 1 0 

25. 957 

12.342 

19.071 

-7. 839 

-3.604 

-0 . 2 86 

0.05 

27.037 

12.032 

22.593 

-7.017 

-3. 620 

- 0 . 2 5 0 

0.05 

27.037 


22.593 


-3. 62 0 


0.04 

28. 165 


19.969 


-3.757 


0.0 3 

28.711 


20.146 


-3.400 
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TAbLC 34 


STATION 4 

r (in) 

4. 65 

4.4 5 
4. 25 

4.05 

3.05 

3.65 
3.4 5 
3. 25 
*.05 
2. 85 

2.65 
2.45 
2.20 
2 . 00 
1. 80 
1.60 
1 . 40 
1.20 
} . 00 
0 . 00 
0.60 
0.40 
0 . 20 
0 . 10 
0. 05 
0.05 
0. 04 
0. 03 


CONSTANT PH CSS 0 1 0 I LOW ALONG CONCA VK WALL 



X (IN) 

U**2 

V**2 

2.53 0 

3.40 2 

3.700 

5.255 

5.550 

7.656 

8.600 

11.419 

13.1 30 

15. 324 

16.650 

1 6 • b 5 6 

1 8.380* 

17. 139 

1 0.070 

15. 454 

1 0 .100 

17. 336 

1 8.39 0 

19.719 

16.870 

20.656 

1 6.000 

21.573 

15.214 

20.458 

14.493 

21 . 28 3 

14.847 

21 . 885 

15.611 

17, 91 M 

16 . 076 

19. 986 

16. 162 

19.973 

15.853 

19.91 6 

1 5 . Q 4 0 

18.906 

16.546 

16.874 

19.124 

13.810 

24 . 427 

9.155 

29.209 

7.386 

30,481 

7.798 

30.481 


29.747 


29.977 



96 . no 


W * * 2 

rJV 

3.550 

- 1 .065 

4.659 

- 1 . 93 9 

6.050 

-3.149 

8 . 2 50 

-5.390 

10.530 

-8. 387 

12. 100 

-9.587 

13.250 

-9.655 

14.63 0 

-7.436 

16.250 

-7.644 

18.370 

-7.690 

20.709 

-6.504 

22 . 650 

- 6 . 04 9 

21 .952 

-5.46 1 

21 . 629 

-7.362 

22.323 

-9.611 

22.609 

-3.569 

22 . 880 

- 1 0. 158 

2 1.301 

-10.993 

23.669 

-11. 027 

24.485 

-10.777 

24.583 

-9.765 

21 . 845 

-8.845 

20.899 

-7.552 

19. 581 

- 6 . 4 9 6 

1 9. 05'> 

-7 . 130 

19. 059 


1 9.90 3 


19.251 



6 (IN) 

= -21.870 

n v 

vw 

* . 300 

-1.713 

( , 66\ : 

-2. 34 * 

1. 27 0 

- 3 . H r > 

2.190 

-4 . 398 

2. 8 70 

-6.013 

3. 07 0 

-5.48 > 

2.850 

-5. 1 26 

2. 30 0 

-3 . 304 

1 . 07 0 

-2.74? 

0. 20 0 

-0.773 

0. 400 

-1.247 

1. 090 

-0.63 3 

2. 329 

-1 . 939 

2.675 

-2.097 

3. 2 19 

-3.442 

3. 367 

-4.615 

3. 557 

-5.1 35 

3.317 

-4 . 902 

3. 086 

-4.233 

3.40 4 

-4.202 

4.169 

-4. 7 7 u 

4.932 

l 

UJ 

• 

■D 

5.8 06 

-3.712 

4 . «7 3 

-2.. 863 

4. 280 

- 4 . 10 9 

4. 280 


4.244 


4. 322 
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TABLK i r j 


CONSTANT L’H ESS Ur>. E FLOW ALONO CONCAVE WALL 


STATIC) ?I 5 


X(IN) - 

9 6.00 


r (IN) - 

-21 . h75 

Y (IN) 

0**2 

V ** 2 

W**2 

UV 

OW 

VK 

3.00 

2.113 

3.525 

4.795 

-1.195 

'■.418 

-0. 297 

2.80 

2.405 

3.524 

5.117 

-1.144 

0. 350 

0. 1 90 

2.60 

2.56 3 

4.091 

5. 985 

-1.437 

0 . 195 

0.512 

2. 40 

3.176 

4. 594 

6.094 

-1.629 

0.215 

0.63 M 

2.20 

4.275 

6. 59 3 

6.84b 

-2.927 

r . 13 5 

0.779 

2.00 

5.6 96 

7. 581 

7.741 

-3.709 

-0.074 

0.941 

1.80 

8.852 

9.756 

8.05o 

-3.572 

-0 . 5 b 2 

1.2 69 

1 .60 

11.180 

12.720 

1 0. 948 

-7.960 

-0.623 

1.241 

1 . 40 

1 3.699 

14.683 

12.510 

-9.442 

-0.467 

0.74 2 

1 .20 

15.007 

15. 249 

14.629 

-9. 90 4 

- 0 . 275 

0.4 62 

1 . 00 

15.777 

16. 130 

15. 549 

-10.268 

0.033 

-0. 1 54 

0.8 0 

15.939 

17.774 

16.125 

-10.6 19 

0.093 

-0. 31 R 

0.60 

15.871 

15. 940 

16.896 

-8.924 

-C. 354 

-0.091 

0. 40 

16.536 

12.503 

17.886 

-7.145 

-0. 938 

-0.243 

0.30 

18.417 

14 . 046 

19.109 

-8.366 

-1.955 

1.141 

0.20 

2 1 . 998 

13. 590 

20.270 

-8.578 

-1 . 989 

0.3 62 

0. IS 

24.501 

9.4 45 

22.179 

-7 . 045 

-2. 881 

0.826 

0 . 1 0 

26.907 

12.731 

21.805 

-9.51 1 

-2. 31 3 

-0.965 

0.05 

28.437 

5.068 

23.724 

-5.634 

-1. 845 

-1 . 584 

0.05 

28.437 


23.729 


-1 . 845 


0.04 

28. 769 


2 3.128 


-1. 654 


0.03 

28.963 


23. 100 


-1. 679 
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